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Today, in any major war, an enemy 
would be protected by an almost 
impenetrable curtain of radar and 
guided missiles. But there is a gap 
beneath the curtain; a vulnerable level. 
The NA.39 is designed and 
constructed to strike at this vulnerable 
level and that is why the Ministry 

of Defence has stated— 


* “In the low-level strike role, 
the NA.39 is ahead of 
any other aircraft in the world ” 


blackoury | BLACKBURN AIRCRAFP LIMITED, BROUGH, YORKSHIRE 
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One of the largest manufacturers of motive power units 
in the world, Bristol Siddeley Engines Limited produce 
the Gamma. A liquid propellent rocket engine, the 
Gamma delivers 16,400-Ib thrust (7,438 kg) at sea level 
rising to 19,000 lb (8,618 kg) outside the earth’s atmos- 
phere for a total engine bay weight of under 700 Ib. 

The Bristol Siddeley Gamma has four gimbal-mounted 
combustion chambers which are hydraulically actuated 
for vehicle guidance. Each combustion chamber is fed 
with propellents by its own turbopump unit and the 
four units are joined at the centre by a common 
manifold. The Gamma burns hydrogen peroxide (HTP) 
with kerosene and uses silver-plated nickel gauze as a 
catalyst to decompose the HTP into oxygen and super- 


heated steam. 
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End view showing combustion 
chambers which can be 


ENGINEERING ADVANCE 
BRISTOL SIDDELEY 


in- 


Gamma powers Black Knight 


The Bristol Siddeley Gamma powers the Saunders-Roe 
Black Knight, Britain’s highly successful space research 
vehicle. The Gamma has proved itself to be exceptionally 
reliable. In fact, in all firings to date Black Knight has 
never failed to start, and has reached a height of over 
500 miles above the Woomera rocket range in 
Australia. 


Since Bristol Siddeley’s rocket division began work in 
1946 it has developed a wide range of rocket components. 
By combining these components in single or multi- 
chamber layouts, thrust requirements from 500 Ib up 
to very high figures, can be met. 


mm BRISTOL SIDDELEY ENGINES LIMITED 
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...looking ahead, 
Blackheath are using 


the most up-to-date 


equipment for the 


continued production of 


SUPERIOR DROPFORGINGS 


THE BLACKHEATH STAMPING CO. LTD. 


BLACKHEATH BIRMINGHAM Tel: Blackheath 1448/9 


CABIN AIR 
CONDITIONING CONTROL 


ANTI-ICING CONTROLS 


FUEL TANK 
PRESSURISATION 


d AN AS VALYV. 
Anti-icing and de-icing equipment manufactured by 


Teddington Aircraft Controls include components used on windscreen 
de-icing, wing and tail de-icing, engine intake cowl anti-icing, ELECTRO-MAGNETIC VALVES 

fuel filter ice warning, bomb bay heating and 

propeller and spinner de-icing systems in many aircraft at present on 
military service or airline operation. These systems 

make use of hot air or gas, or electrical heater mat, and the controls 
are either electrically or pneumatically operated. PRESSURE SWITCHES 


TIME SWITCHES 


PRESSURE REGULATORS 


POTENTIOMETERS, 


TEDDINGTON AIRCRAFT CONTROLS LTD an 


AEGD. TRADE MARK PRESSURE TRANSDUCERS 
MERTHYR TYDFIL, SOUTH WALES Telephone: Merthyr Tydfil 3261 


LONDON OFFICE: COLNBROOK BY-PASS, WEST DRAYTON, MIDDLESEX Telephone: Colnbrook 2202-3-4 TAC 162 
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Hawker 
Siddeley 
Group 


IN THE FOREFRONT 
OF WORLD AVIATION 
AND INDUSTRY 
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LEADING THE WAY 


TO THE STARS 


Steels developed by Firth Brown for the Aircraft 
Industry were assisting world famous Engine and Airframe 
designers long before the first Schneider race. Unceasing 
research and advanced metallurgical control has enabled 
Firth Brown to produce a range of steels to meet the 


most rigorous requirements of Aeronautics and 


Astronautics alike. 


SPECIAL ALLOY STEELS 


THOS FIRTH & JOHN BROWN LIMITED . SHEFFIELD 
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TRAVEL TIME ALMOST HALVED 707 


Now Flying Westwards 
LONDON + NEW YORK 
SAN FRANCISCO+ HONOLULU - 
NEW ZEALAND « SYDNEY 


*Quick and easy connections to Auckland from Fiji 


Also Flying Eastwards 
LONDON MIDDLE EAST 
PAKISTAN + INDIA THAILAND 


MALAYA + SYDNEY 
\ 


S 
ROUND-THE-WORLD AIRLINE 


WITH B.0.A.C, T.E.A.L ANDO S.A.A 


LZ 


Tickets and helpful advice from all appointed Travel Agents or Qantas, corner of Piccadilly and Old Bond Street, London, W.1. (Mayfair 9200) 
or any office of B.0.A.C. 


JOURNAL OF THE ROYAL AERONAUTICAL SOCIETY 7 (ADVERTISEMENTS NOVEMBER 1959 


| ON | 

DON To | 

|| S$ A 

N FRANC | 

1 
wi thout h SCO Il 

OF tye «ging py | 


prizewinner in 1909 * 


In 1909 Elektron Magnesium Alloys 
were introduced to the world at an 
international exhibition in Germany and 
were awarded the First Prize in a 
competition for light metals. 


still the world’s lightest 
structural metal in 1959 


The past 50 years have more than justified this initial success of Elektron. 
Today the latest Elektron alloys, developed by MEL to meet the exacting 
requirements of industries as diverse as aircraft and textiles, motor transport and 
nuclear engineering, are still the finest magnesium alloys available anywhere in 


the world. 


In this Golden Jubilee year of Elektron, MEL are proud of the part they have 
played in pioneering the use of magnesium in Great Britain. As sole producers 

of the metal, with unique research, development, prototype and fabrication 
facilities, MEL enterprise and technological skill will continue to make available 
to all industries the advantages of magnesium—the world’s lightest structural metal. 


Magnesium Elektron Ltd. 


Clifton Junction, Manchester 
London Office: 5 Charles II Street, St. James's, SWI 
Magnesium Elektron, Inc., New York 20, USA 
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To make a ground-to-air guided weapon system mobile is not simply a matter of 


jacking it up and running some wheels under it. It needs an entirely different — and 


more challenging—approach. We know—for we were chosen by the British Govern- 
ment to concentrate on this very problem. And five years of design and development 


work have gone into the solution. THE ANSWER IS 


THE ‘ENGLISH ELECTRIC: THUNDERBIRD 


—a complete mobile and air transportable weapon system — fully de- 
veloped—now in service with the British Army. 


Its mobility and ease of redeployment give the flexibility that is so essential to meet the nation’s 
ever-changing defence requirements .. . AND THE ECONOMY AS WELL 


ENGLISH ELECTRIC AVIATION LIMITED - GUIDED WEAPONS DIVISION: LUTON + STEVENAGE « WOOMERA 


A MEMBER OF THE ENGLISH ELECTRIC AVIATION 
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“Maintaining direct radio-link with 
their refuelling vehicles out on the 


Se the Shell and BP Aviation 
rvices movement control room 
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AVIATION SERVICES 


Shell-Mex and B.P. Ltd. Shell-Mex House, Strand, W.C.2, 


Registered Users of Trade Marks. Distributors in the 
United Kingdom for the Shell, BP and Eagle Groups t 


4 OUT OF 5 AIRCRAFT * 
USING LONDON AIRPORT 
USE SHELL OR BP 
AVIATION FUELS 


‘Careful planning and teamwork is the order of the day," 
says Leslie Hughes, operations room Movement Controller for 
Shell and BP Aviation Services. ‘‘Close co-ordination between 
ourselves and London Airport's central control tower means 
that we know well in advance just where and when an aircraft 
is landing. By consulting our Aircraft Parking Area Plan we 
can direct the fueller carrying the appropriate grade of 
aviation fuel to meet the aircraft as it taxies on te the apron. 

It is this organization that enables us to cut down to a 
minimum the turn-round time of the aircraft we fuel.” 


* 25 Leading Airlines use the Shell and BP Aviation Services at London Airport. 
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Promised V810 economy 
more than proved in practice! 
The economical performance promised by this 
graph has been more than fulfilled in the first six 
months’ operation of Continental Airlines’ V810 
fleet. 


VICKERS 


AIRCRAFT-MILE COSTS 
II cents below forecast! 


—That’s proof of 


Economy for you! 


After 6 months’ V8I0 operation 
CONTINENTAL AIRLINES REPORT 


Aircraft mile cost 71. 23 cents for an average sector of 500 miles 
33 
— II CENTS BELOW FORECAST ! 


Outstanding V810 economy and passenger appeal brought the airline 
$1.3 MILLIONS in Viscount profit between May and Decem- 
ber 1958 (Average fleet of nine aircraft for the period). 


The Viscount V810 is still without a challenger in its class as 
a modern, medium-sized short/medium haul airliner. Its econ- 
omy, flexibility, 365 m.p.h. cruise speed, and background of 
Vickers engineering and Rolls-Royce turbo-prop experience in 
over 370 Viscounts in service all over the world, make it unique 
as a profit-earner in this field. 


Low fare/High profit note 

The V810 in a 70-seat layout has a seat-mile cost of only 1.34 
cents at 300 miles, and a proved 1.01 cents at 500 miles! With 
only 49 passengers, this aircraft would clear 10%, profit at fares 15% 
down on current U.S.A. coach fares! 


FOUR ROLLS-ROYCE DART TURBO-PROP ENGINES 


e@ Fastest and most economical medium-sized airliner in the world 


VICKERS-ARMSTRONGS (AIRCRAFT) LIMITED WEYBRIDGE SURREY 
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For Superlative Accuracy .... 
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& Cross British made products are manu- 


factured by an outstanding process used 
for the hardening and tempering of steel 
coils and rings. Covered by many patents, 
this manufacturing method has enabled 
components to be made from wire and 
with superlative accuracy. When top 
standards are essential you can rely on 
Cross precision products. 


CROSS MANUFACTURING COMPANY (1938) LIMITED 
BATH - SOMERSET - ENGLAND ~- Phone: Combe Down 2355/8 Telegrams: ‘CIRCLE’ BATH 


This photograph shows a few examples of high-quality precision springs, 

made by Riley. Some were produced to customers own specifications, others 
were the result of recommendations by the Riley design and research department. 
If you have any kind of spring problem, Riley will find the answer to it. 


(The diagram shows information required when ordering this type of spring.) 
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There are no finer springs than Springs by 


ROBERT RILEY LTD., MILKSTONE SPRING WORKS, ROCHDALE, Tel: ROCHDALE 2237 (Slines) Grams: *RILOSPRING’ ROCHDALE TELEX 63-151 
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Whether it’s at one end, as in the quietly efficient Westland Wessex . . . 
whether it’s at both ends, as in the outstanding Bristol 192 . . . whether it’s 
upright, horizontal or in any position in between, the Napier Gazelle free 
turbine engine will always satisfy the exacting demands of helicopter 
operation. Designed specifically for helicopter application, the Napier 
Gazelle, with its “‘any angle” installation, relieves the designer of many 
accepted limitations and gives him more opportunity to make proper 


use of fuselage space. 


NAPIER GAZELLE FREE TURBINE ENGINE 


D. NAPIER & SON LIMITED, LONDON, W.3. A MEMBER OF THE ENGLISH ELECTRIC AVIATION GROUP Seo 
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INTERNATIONAL SERIES ON 
AERONAUTICAL SCIENCES AND SPACE FLIGHT 


Joint Chairmen 


DR. THEODORE VON KARMAN ___ DR. HUGH L. DRYDEN 


This series forms a complete source of study and reference covering the whole field of aero- 
nautical sciences and space flight. Each book covers a specific aspect of the subject, and serves 
as an authoritative review for the well-prepared beginner and a reliable source book for the 
experienced worker. 

Together with the works published by Pergamon Press on behalf of the Advisory Group 
for Aeronautical Research and Development (AGARD), North Atlantic Treaty Organization, the 
volumes in Aeronautical Sciences and Space Flight meet much of the increasing demand for up- 
to-date and authoritative literature prepared by those who fully understand the needs of the 
student and the experienced worker. 


Board of Specialist Editors Title of Division in the Series 

R. L. Bisplinghoff, W. S. Hemp SOLID AND STRUCTURAL MECHANICS 

R. T. Jones, W. P. Jones AERODYNAMICS 

A. D. Baxter, E. R. Sharp PROPULSION SYSTEMS INCLUDING FUELS 
R. j. Lees, E. Rechtin AVIONICS 

G. Melvill Jones, W. R. Lovelace II AVIATION AND SPACE MEDICINE 

M. C. Wilson FLIGHT-TESTING 

M. Alperin, S. F. Singer, M. Stern ASTRONAUTICS 

W. Jj. Harris, Jr., C. Zwikker MATERIALS, SCIENCE AND ENGINEERING 


Individual Editors appointed for each 
Symposium SYMPOSIA 


Some of the titles in preparation: 


RAREFIED GAS DYNAMICS ELASTIC STABILITY OF PLATES AND 
Proceedings of the First International Sym- SHELLS 

posium held at Nice by H. L. Co 

Edited by F. M. Devienne 

theory, as well as experimental aspects of the 

subject, of interest to all concerned with | FLAME PROPAGATION 


aviation and aerodynamics. £6 net ($17.50) by D. M. Simon 


AGARD FLIGHT TEST MANUAL 
Edited by C. D. Perkins, E. J. Durbin and E. Seckel 
(Four Volumes) 

An authoritative, comprehensive up-to-date handbook on the analysis and 
testing of aircraft performance, stability, control and instrumentation. 
Published for the Advisory Group for Aeronautical Research and Develop- 
ment, North Atlantic Treaty Organization, the AGARD FLIGHT TEST 
MANUAL will comprise four loose-leaf volumes, and to keep it completely 

up-to-date annual supplements will be available from 1960 onwards. 
£30 net ($100.00) for 4 vols. 
Subscription to Supplements: £5 5s. per annum 


PERGAMON PRESS London New York Paris Los Angeles 


4G 5 Fitzroy Square, London W.1. 122 East 55th Street, New York 22, N.Y. 
24 Rue des Ecoles, Paris V°. 
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B.O.A.C. H.Q., London Airport. 
Middlesex. 


LUTON 
President: F. H. Potzicutt, F.R.Ae:S. 
Chairman: H. J. Brooxs, A.F.R.Ae.S. 
Hon. Secretary: L. A. 
A.F.R.Ae.S., Hunting Aircraft 
Ltd., Luton Airport. 


MANCHESTER 
President: Sir Roy H. Dosson, C.B.E., 
J.P., Hon.F.R.Ae.S. 
Chairman: C. E. Fievpinc, O.B.E., 
F.R.Ae.S. 


Hon, Secretary: J. A. E. WATERFALL, 
56 Manor Avenue, Asbton-on- 
Mersey, Cheshire. 


MERTHYR TYDFIL 

President: C. S. GARDNER. 

Chairman: J. INGRAM. 

Hon. Secretary: Dr. G. O. FORESTER, 
Chief of Research, Teddington 
Aircraft Controls Ltd., Cefn Coed, 
nr. Merthyr Tydfil, South Wales 


PRESTON 
President: Sir GEORGE H. NELSON. 
Chairman: E. Lovetess, A.F.R.Ae.S. 
Hon, Secretary: J. F. Foss, 
A.F.R.Ae.S., English Electric 
Aviation Ltd., Warton Aerodrome, 
nr. Preston. 


READING 
President: Sir FREDERICK HANDLEY 
Pace, C.B.E., Hon.F.R.Ae.S. 
Chairman: R. J. Fenner, A.F.R.Ae.S. 

Hon. Secretary: N. KITCHING, 
A.R.Ae.S., Handley Page (Reading) 
Ltd., The Aerodrome, Woodley, 
Reading. 


SOUTHAMPTON 
Chairman: Professor E. J. RicHARDS 
O.B.E., F.R.Ae.S. 
Hon. Secretary: F. C. KIRK PATRICK, 
A.R.Ae.S., Folland Aircraft Ltd., 
Hamble, Southampton. 


SOUTHEND 
Chairman: T. C. BERRETT, 
A.R.Ae.S. 
Hon, Secretary: K. A. Cote, A.R.Ae.S., 
c/o Aviation Traders (Eng.) Ltd., 
Southend Airport, Essex. 


SWINDON 
Chairman: S. R. HuGuHes, A.F.R.Ae.S. 
Hon. Secretary: D. RicG, A.R.Ae.S., 
Vickers-Armstrongs (Aircraft) Ltd., 
Supermarine Works, South Mars- 
ton, Swindon. 


WEYBRIDGE 
President: Sir GEORGE EDWARDS, 
C.B.E., F.R.Ae.S. 


Chairman: H. GARDNER, F.R.Ae.S. 
Hon. Secretary: E. G. BARBER, 


A.F.R.Ae.S., Vickers-Armstrongs 
(Aircraft) Ltd., Weybridge. 
YEOVIL 


President: E. C. WHEELDON. 

Chairmen: D. L. HOLtis-WILLIAMS, 
F.R.Ae.S. and Dr. E. W. STL, 
F.R.Ae.S. 

Hon. . Secretary: L. A. LANSDOWN, 
A.F.R.Ae.S., Westland Aircraft 
Ltd., Yeovil. 
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NOTICES 


CHRISTMAS HOLIDAY 
The Office and Library of the Society will be closed 
from Thursday afternoon 24th December until 9 a.m. on 
Tuesday 29th December. 


THE LIBRARY 

During rebuilding at the Offices of the Society the 
Library will be operating under some difficulty. Members 
and other users are therefore asked to make allowances 
should the service fall slightly short of its usual standard 
for some months. 

The Library will be closed on 8th, 9th and 10th 
December for the Associate Fellowship Examinations. 


PRIZES IN AERONAUTICS 

The Society’s Prizes in Aeronautics, awarded annually 
to the best student in the Aeronautical Department at 
Universities and Colleges organising courses in aeronautics, 
have been awarded for 1959 as follows 

Bristol University—Joint Award: J. W. Britton and 
D. R. Philpott. 

Cambridge University—C. G. B. Mitchell of St. John’s 
College. 

Glasgow University—Kenneth M. Nicoll 

Imperial College, London—B. N. Tomlinson and M 
Lacey. 

Queen's University, Belfast—Kenneth J. Holden. 

Northampton College of Advanced Technology, 
London—G. N. Sage. 

Queen Mary College, London—G. J. Shorey. 

Southampton University—R. D. Ford. 


HENLOW CADET PRIZE 
The Henlow Cadet Prize for 1959 has been awarded to 
Pilot Officer D. R. West. The Prize is three years’ sub- 
scription as a Graduate of the Society, with entrance fee, 
and a book. 


ANNUAL SUBSCRIPTIONS 


Members are reminded that their annual subscriptions 
become due on Ist January 1960. The rates are: 


Home Abroad 

Fellows 8 0 0 7 00 

Fellows (over 65) 5 0 0 400 

Associate Fellows 6 0 0 5 0 0 

Associate Fellows (over 65) 40 0 3 0 0 

* Associates 5 0 0 
Associates (over 65) 3; @ g 3 00 

Graduates (aged under 26) 210 O 210 O 
Graduates (aged 26 and over) 3 0 O 3 0 0 
Students (aged under 21) 110 O 110 O 
Students (aged 21 and over) 20 0 2 0 0 
Companions 40 0 4 0 0 
Companions (over 65) i 3 0 0 
Companions (No Journal) 210 O 210 O 


*Any Associate elected before Ist October 1947 may, 
if he wishes, elect not to receive the JOURNAL, and in this 
case his subscriptions will be reduced by £1 10s. Od. to 
£3 10s. Od. 


It will avoid delay and confusion if members, when 
sending remittances for subscriptions, will state their names 
clearly and give their addresses and grades of membership. 
Remittances should be made payable to the Royal Aero- 
nautical Society. 


LECTURE THEATRE 

Work on the construction of the Lecture Theatre was 
started at the beginning of October and good progress is 
being made. 

The complete list of contributions received since the 
October JoURNAL went to press will be published in the 
December issue 

Donations are still needed. 


SYMPOSIUM ON FLIGHT SAFETY 

Admission to the All Day Discussion on Flight Safety 
to be held at the Institution of Mechanical Engineers on 
11th December 1959, will be by special ticket only for 
members and non-members. Application for tickets should 
be made to the Secretary of the Society and until 4th 
December will be issued to members only; thereafter any 
seats available will be issued to non-members in strict 
rotation 

The Introductory papers will be: The Accident 
Investigator’s Viewpoint by Group Captain P. G. Tweedie 
(M. of A.); The Pilot's View by Captain E. C. Miles 
(B.O.A.C.); The Operator’s View by Captain J. W. G. 
James (B.E.A.); and The Constructor’s View by R. H. 
Warde (B.A.C.). 

The Discussion will begin at 9.45 a.m. 


HONOURS AWARDED TO MEMBERS 

Sir Arnold Hall, Immediate Past President of the 
Society has been awarded the Jr A. G. von Baumhauer 
Medal for the period 1954-1959 

This is awarded periodically to the person who, in the 
opinion of the Council and advisers of the Royal Nether- 
lands Aero Club, has made the most valuable contribution 
in either a scientific or practical sense to flight safety. 


Sir George Edwards, Past President of the Society has 
been awarded the Daniel Guggenheim Medal for 1959 for 
a lifetime devoted to the design of military and commercial 
aircraft culminating in the successful introduction into 
worldwide commercial service of the first turbine-powered 
propeller-driven aircraft. 


de Havilland Aircraft Co. Comet. The Elmer A. 
Sperry Award for 1959 was presented to the de Havilland 
Aircraft Co. on the 7th October in New York during the 
7th Anglo-American Aeronautical Conference. The award 
is made annually for a distinguished contribution which 
through application, proved in actual service, has advanced 
the art of transportation, whether by land, sea or air. 
de Havilland have successfully developed the world’s first 
jet-powered air liner, the de Havilland Comet. 

Specifically named in the citation are Sir Geoffrey 
de Havilland (Honorary Fellow}, Mr. Charles C. Walker 
(Honorary Fellow), and the late Major Frank B. Halford. 


MESSAGE FROM A PAST PRESIDENT OF THE SOCIETY 
Sir Arthur Gouge, Past President of the Society has sent 
a note to say that he is retiring and would like all his old 
friends to know. He wishes to take this opportunity of 
thanking them for the help and co-operation given to him 
in the past. 


THE SHUTTLEWORTH COLLECTION 


Members may be interested to know that the Shuttle- 
worth Collection is open from Monday to Friday by 
appointment, from 9 a.m. to noon and 2.30 p.m. to 4.30 p.m. 
and on some Saturdays from 9 a.m. to 11.30 a.m. 

Applications for visits, enclosing a stamped addressed 
envelope, should be sent to L. A. Jackson, Esq., Old 
Warden Aerodrome, Biggleswade, Bedfordshire. 
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Summaries of Lectures 


Fifteenth British Commonwealth Lecture 
to be given on 19th November 1959 


SOME RECENT PROGRESS IN AIR SURVEY WITH 
SPECIAL REFERENCE TO NEWLY DEVELOPED 
TERRITORIES 


by W. P. SMITH, M.B.E., B.A., F.R.LC.S. 
(Director, Fairey Air Surveys Ltd.) 


The paper considers briefly the state of air survey 
immediately after the Second World War and then deals 
in some detail with important aspects of technical progress. 

These include the development in camera and lens 
design leading to improvement in photographic quality, in 
navigational methods, and to a lesser extent, in aircraft 
available. The paper then reviews progress in the appli- 
cation of aerial photography to problems of geology, 
agriculture and forestry in newly developing territories. 
This is followed by an appreciation of the value of mapping 
by air survey methods in connection with projects overseas 
and deals particularly with hydroelectric programmes in 
Africa, the survey aspects of land tenure problems, and 
large scale contouring for mining and other purposes. 
These lead to a consideration of progress made in the use 
of aerial triangulation and modern plotting machines and 
to the desirability of reducing the need for ground control. 
The paper discusses airborne electronic methods such as 
Shoran, Hiran and the Airborne Profile Recorder, all of 
which may introduce savings in time and cost when map- 
ping remote areas. More recent systems such as the 
Tellurometer and Aerodist are touched upon. 

The third section of the paper is concerned with geo- 
physical surveys used in mineral exploration and includes 
reference to recent developments in surveys for radio 
activity, magnetic anomalies, and conductivity. The recently 
developed two-aircraft system for electro-magnetic surveys 
is described and the paper ends with an attempt to forecast 
some of the directions in which progress may be made in 
the near future. This section includes references to new 
“super-wide-angle” aerial cameras, to analytical methods 
of stereo-restitution, to automatic map plotting, mapping 
by radar photography and finally to gravimetric survey 
from the air. 


Third Lanchester Memorial Lecture 
to be given on 26th November 1959 


SOME DEVELOPMENTS IN BOUNDARY LAYER 
RESEARCH IN THE LAST THIRTY YEARS 


by Professor Dr. H. SCHLICHTING 
(Technical University of Braunschweig) 


Boundary layer theory is the cornerstone of our know- 
ledge of the flow of air and other fluids of small viscosity 
under circumstances of interest in many engineering appli- 
cations, especially also in aeronautics. Many complex 
problems in aerodynamics, as for instance the problem of 
skin friction, which was theoretically attacked very early 
by F. W. Lanchester, have been clarified by studying the 
flow within the boundary layer and its effects on the 
general flow around the body. 

Research work on boundary layers, as started by 
Prandtl in 1904, was for the first twenty years—up to 
Prandtl’s Wilbur Wright Memorial Lecture to the Royal 
Aeronautical Society in 1927—almost entirely restricted to 
Prandtl’s Institute at Géttingen. But since about 1930 
boundary layer theory has been generally accepted, and in 
the past thirty years there has been an almost exponential 
rise of the number of contributors to its further develop- 
ment. 
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The author tries to trace certain lines, along which this 
important branch of modern fluid dynamics has developed 
in the last thirty years. In this connection the following 
topics are treated to some extent: 


(1) Transition from laminar to turbulent flow. 
(11) Boundary layer control for high lift and low drag. 
(III) Aerodynamic heating at high speeds (high Mach 
numbers). 
(IV) Boundary layer effects on swept wings and on 
rotating bodies. 


(1) The theoretical investigations of the problem of 
transition start from Reynolds’ and Lord Rayleigh’s 
hypothesis of the instability of laminar flow. After 
many unsuccessful attempts Tollmien in 1930 
finally succeeded in calculating the critical Reynolds 
number for the boundary layer on a flat plate. More 
than ten years later Tollmien’s stability theory 
was completely confirmed by very careful experi- 
ments of Dryden and his co-workers. 

After many wind tunnel experiments the investiga- 
tions of boundary layer control for high lift of 
aerofoils led to the construction of two aeroplanes 
with boundary suction at the Aerodynamische 
Versuchsanstalt Géttingen in 1938. Later on also a 
considerable reduction of the skin friction of aero- 
foils was obtained with the advent of the laminar 
flow aerofoil. 

In flow at high Mach numbers the velocity boundary 
layer is accompanied by a thermal boundary layer 
which is caused by frictional heating. The large 
increase in the temperature of a solid surface in a 
high speed stream which can be calculated from 
boundary layer theory only, poses a serious prob- 
lem to aeronautical engineers (“thermal barrier”). 
The aerodynamic characteristics of swept wings and 
delta wings are largely governed by the behaviour 
of their boundary layer. Some of the draw-backs 
of such wing plan forms can be remedied by 
boundary layer control, as for instance by a 
“boundary layer fence.” For turbomachines the 
influence of the centrifugal forces on the boundary 
layer plays an important role for their aerodynamic 
coefficients. 


DIARY 


LONDON 


13th November 
AGRICULTURAL AVIATION Group MeetTING—General Dis- 
cussion on the International Conference at Cranfield. 
Library, 4 Hamilton Place. 7 p.m. 

17th November 
Lecrure—Magnetogasdynamics. Dr. J. A. Shercliff. 
Library, 4 Hamilton Place. 7 p.m. 

19th November 
FIFTEENTH BRITISH COMMONWEALTH LECTURE—Some 
Recent Progress in Air Survey with Special Reference to 
Newly Developed Territories. W. P. Smith. The Institu- 
tion of Mechanical Engineers, 1 Birdcage Walk, S.W.1! 
6 p.m. (Tea at 5.30 p.m.) See summary. 

24th November 
Lecture—Heat and Mass Transfer in Aeronautical Engin- 
eering. Professor D. B. Spalding. Library, 4 Hamilton 
Place. 7 p.m. 

26th November 
THIRD LANCHESTER MEeEmorRIAL LeEcTURE—-Some Develop- 
ments in Boundary Layer Research in the last Thirty Years. 
Professor Dr. H. Schlichting. Church House, Westminster, 
S.W.1. 6 p.m. (Tea at 5.30 p.m.) See summary. 

27th November 
MAN PoweRED ArRCRAFT GRouP LecTuRE—Man as an 
Aero Engine. Dr. D. R. Wilkie. Library, 4 Hamilton Place. 
7 p.m. 
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Ist December 
LecrurE—Transonic Flow over Swept Wings. E. W. E. 


Rogers and Dr. I. M. Hall. 4 Hamilton Place. 7 p.m. 
3rd December 
MAIN LecTURE AT GLOUCESTER AND CHELTENHAM BRANCH 
Some Servo-Mechanisms for Supersonic Aircraft. S. G. 
Glaze. St. Mary’s College, Cheltenham. 7.30 p.m. 
8th December 
LecTurE—Some Practical Aspects of Compressor Blade 
Vibration. Dr. E. K. Armstrong and R. E. Stevenson. 
Library, 4 Hamilton Place. 7 p.m. 
December 
SYMPOSIUM ON FLIGHT SaFETY—AIll-day Discussion. Institu- 
tion of Mechanical Engineers, Birdcage Walk. S.W.1. 
945 a.m. See Special Notice. 
AGRICULTURAL AVIATION Group LECTURE—Aerial Opera- 
tions Planning in the United Kingdom: Results and 
Difficulties. Peter King. Library, 4 Hamilton Place. 7 p.m. 
17th December 
The Engineering Advantages of Rear Engines in Pure Jet 
Transport. E. S. Allwright. Institution of Mechanical 
Engineers, Birdcage Walk, S.W.1. 6 p.m. (Tea at 5.30 p.m.). 


GRADUATES’ AND STUDENTS’ SECTION 

20th November 
Winter Dance. 4 Hamilton Place. 

2nd December 
Aerodynamics of Racing Cars. F. A. Costin. Library, 4 
Hamilton Place. 7.30 p.m. 

BRANCHES 

16th November 
Henlow—Flight Development of a Modern Prototype 
Aircraft. C. F. Bethwaite. The Assembly Hall, R.A.F. 
Technical College, Henlow. 7.45 p.m. 

17th November 
Christchurch——Lecture by Maurice Brennan. Kings Arms 
Hotel. 7.30 p.m. 
Gloucester— Prospects for Space Flight in the Common- 
wealth. J. E. Allen. Wheatstone Hall, Brunswick Road. 
7.30 p.m. 

18th November 
Coventry—General Aircraft Electrical Systems. G. G 
Wakefield. Herbert Art Gallery. 7.30 p.m 
Preston—Aircraft Production Techniques. L. G. Burnard. 
R.A.F.A. Hall, Preston. 7.30 p.m 
Southampton—-Supersonic Propulsion. Dr. R. R. Jamison 
Engineering Lecture Theatre, Univ. of Southampton. 8 p.m 
Weybridge—-Flight Testing at Supersonic Speeds. Wing 
Cdr. R. P. Beamont. Apprentice Training School, Vickers- 
Armstrongs (Aircraft) Ltd. 6.10 p.m. 

20th November 
Luton— Annual Dinner. Leicester Arms Hotel 

21st November 
Hatfield Visit to N.P.1 

24th November 
Belfast-Film Night. Lecture Hall LG8&, David Keir 
Building, Queens University. 7 p.m 

25th November 
Bristol—The Training of a Production Engineer for the 
Aircraft Industry Professor J. Loxham. Main Lecture 
Theatre, Faculty of Engineering, University of Bristol. 
7 p.m. 
London Airpert— Air Traffic Control. L. Woodruffe. B.E.A. 
Viking Centre Cinema, London Airport. 6 p.m. 

28th November 
Haltca—Annual General Meeting and Dinner. Branch 
H.Q., R.A.F. Halton. 3 p.m. 

30th November 
Henlow—Film Night. The Assembly Hall, R.A.F. Tech- 
nical College, Henlow. 7.45 p.m. 

Ist December 
Glasgow—Astronautical Navigation. P. H. Tanner. Royal 
College of Science and Technology. 7.15 p.m 
Luton—Junior Paper/Discussion Paper Evening. Napier 
Senior Staff Canteen, Luton Airport. 6.15 p.m. 
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2nd December 
Bristol— Junior Members Paper Competition. Filton 
House. 6 p.m 
Preston—-The Armstrong Whitworth Argosy. E. D. Keen. 
St. John’s Ambulance Hall, Preston. 7.30 p.m. 
Swindon—Man Powered Aircraft. B. S. Shenstone. The 
College, Victoria Road. 7.30 p.m. 
3rd December 
Gloucester and Cheltenham--Main Society Le&cTURE— 
Some Servo-Mechanisms for Supersonic Aircraft. S$. G. 
Glaze. St. Mary’s College, Cheltenham. 7.30 p.m. 
Isle of Wight—Nuclear Radiation Applications in Industry. 
E. B. Bell. Clubhouse, Saunders-Roe Sports and Social 
Club, Church Path, East Cowes. 6 p.m. 
7th December 
Boscombe Down—Blind Landing. W. J. Charnley. Lecture 
Hall, A. & A.E.E. 5.30 p.m 
Derby—Mechanical Design Aspects of the R.A. 29. F. 
Morley. Rolls-Royce Welfare Hall, Nightingale Road. 
6.15 p.m. 
8th December 
Luton—Young People’s Lecture. 
College of Technology. 
9th December 
Belfast——The Bloodhound. D. J. Farrar. Lecture Hall LG8, 


David Keir Building, Queen’s University. 7 p.m. 
Deck Landing of Naval Jet Aircraft. Lt. Cdr. 


A. W. Bedford. Luton 


Brough 

D. J. Whitehead. Royal Station Hotel Writing Room. 
7.30 p.m, 

Chester—Experiences of an Airline Pilot. Capt. O. P. 


Grosvenor Museum. 7.30 p.m. 


Jones. Lecture Theatre, 
Stanton- 


Christchurch—Saunders-Roe Hovercraft. D. 
Jones. King’s Arms Hotel. 7.30 p.m. 
Hatfield—Wind Tunnels. J. A. Kirk. de Havilland Restau- 
rant. 6.15 p.m. 
Southampton—Experiences of the Trans-Antarctic Expedi- 
tion. Wing Cdr. J. H. Lewis. Engineering Lecture Theatre, 
University of Southampton. 8 p.m. 

Weybridge—-R. K. Pierson Memorial Lecture. H. H. 
Gardner. Apprentice Training School, Vickers-Armstrongs 
(Aircraft) Ltd, 6.10 p.m. 


10th December 
Cheltenham— The Development of Automatic Landing for 


Civil Airline Operations. K. Fernside. St. Mary's College. 


7.30 p.m 

December 
Belfast—Annual Dance. Midland Hotel, Belfast. 

14th December 
Henlow Presidential Address. R. Hills. The Assembly 
Hall, R.A.F. Technical College, Henlow. 7.45 p.m. 

15th December 
Bristol—-Junior Committee Lecture. The Supersonic Air- 
liner. R. G. Thorne. Filton House. 6 p.m. 

16th December 
Coventry—Film Evening. Herbert Art Gallery. 7.30 p.m. 
Preston—The Next Generation of Civil Aircraft. G. H. 
Lee. Queen’s Hotel, Lytham. 7.45 p.m. 

17th December 
Isle of Wight—Annual General Meeting and Film Show. 
Clubhouse, Saunders-Roe Sports and Social Club, Church 


Path, East Cowes. 6 p.m 


JOURNAL BINDING 

Permanent Binding 
1959 Volume (including packing and postage 

in the United Kingdom) 5s. Od. 
Previous Volumes (including packing and 

postage in the United Kingdom) .. aia eS 

Journals, with a note of the name and address of the 
sender, should be sent direct to The Lewes Press, Friars 
Walk, Lewes, Sussex, and the remittance to the Secretary 
at the Offices of the Society. 

Members are asked to be certain that the address to 
which they want their Journals sent is the same on their 
letters to the Lewes Press and to the Society. 
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Copy oF LETTER FROM CHIEF INSPECTOR OF TAXES BRANCH, 
INLAND REVENUE, DATED 17TH OcTOBER 1958 
Reference: C.1./SUB/ 306 


Dear Sir, 

I have to inform you that the Commissioners of Inland 
Revenue have approved The Royal Aeronautical Society 
for the purposes of Section 16, Finance Act, 1958, and that 
the whole of the annual subscription paid by a member 
who qualifies for relief under that Section will be allow- 
able as a deduction from his emoluments assessable to 
income tax under Schedule E of the whole of his annual 
change in the circumstances of the society should occur in 
the future you are requested to notify this office. 

I should be glad if you would inform your members 
as soon as possible of the approval of the society. The 
circumstances and manner in which they may make claims 
to income tax relief are described in the following para- 
graphs, the substance of which you may care to pass on 
to your members. 

Commencing with the year to Sth April 1959, a member 
who is an office holder or employee is entitled to a deduc- 
tion from the amount of his emoluments assessable to 
income tax under Schedule E of the whole of his annual 
subscription to the society provided that— 

(a) the subscription is defrayed out of the emolu- 

ments of the office or employment, and, 

(b) the activities of the society so far as they are 

directed to all or any of the following objects— 

(i) the advancement or spreading of knowledge 
(whether generally or among persons belong- 
ing to the same or similar professions or 
occupying the same or similar positions); 
the maintenance or improvement of stand- 
ards of conduct and competence among 
the members of any profession; 
the indemnification or protection of mem- 
bers of any profession against claims in 
respect of liabilities incurred by them in the 
exercise of their profession; 

are relevant to the office or employment, that is 
to say, the performance of the duties of the office 
or employment is directly affected by the know- 
ledge concerned or involves the exercise of the 
profession concerned. 

A member of the society who is entitled to the relief 
should apply to his tax office as soon as possible after 31st 
October 1958, for form P358 on which to make a claim 
for adjustment of his pay as you earn coding. 


Yours faithfully, 
(signed) T. DUNSMORE, 
Senior Principal Inspector of Taxes. 


ELECTIONS 
The following is a list of elections and transfers of 
membership of the Society :— 


Associate Fellows 
Philip Bruce Atkins 
(from Graduate) 
Greville Gordon Beale 
(from Graduate) 
Geoffrey Bennett 
Angus Lochiel Cameron 
Frank Jones 
Associates 
Kevin Henry Knowles 
Graduates 
William Eric Bodley 
James Dunstan Campbell 
Yat Yung Chan 


Students 
Saeed Ahmed Khan 
Choudhry 
Guy Nils Goodman 
Christopher Martin Leech 


Milton Lalas 
(from Graduate) 
John David Roberts 
(from Graduate) 
Owen James Stevenson 
(from Graduate) 


Robert Wallace Grant 
Alan Douglas Sussex 
Heinz Siegfried Weller 


George Brian Meadley 
Thomas James Triggs 
Donald John Woodman 
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NEWS OF MEMBERS 


Roy ATKIN (Graduate) formerly in the Aerodynamics 
Department, Blackburn Aircraft, Brough, is now in the 
Advanced Projects Group, Hawker Siddeley Aviation, 
Kingston-on-Thames. 

F. J. Brapsury (Associate Fellow) formerly a Research 
and Development Engineer, Automotive Products Co. Ltd., 
is now Chief Engineer, Lockheed Precision Products Ltd. 

H. Brown (Graduate) formerly with Blackburn Air- 
craft Ltd., is now doing Postgraduate Research Work at 
Bristol University. 

L. J. Coates (Associate) formerly Development Engin- 
eer, Hymatic Engineering Co., Redditch, has been 
appointed Experimental Officer, National Engineering 
Laboratory, Mechanisms and Metrology Division. 

G. E. Davies (Graduate) formerly an Aerodynamicist 
on Stability and Control, de Havilland Aircraft Company 
Ltd., is now an Aerodynamicist, Light Aircraft Division, 
Short Brothers and Harland Ltd., Belfast. 

H. J. Davies (Associate Fellow) of the Mathematics 
Department, University of Southampton, is now a Research 
Associate, James Forrestal Research Center on the Hover- 
craft Project, Princetown, New Jersey. 

JoHN E. C. Dickison (Associate) formerly at the 
Ministry of Supply, St. Giles Court has been posted to 
Guided Weapons Dept., Ballistic Missile Division, R.A.E. 

J. DUNHAM (Graduate) formerly a Technical Assistant 
Rolls-Royce Ltd., Derby, has been elected a Junior 
Research Fellow, Pembroke College, Cambridge, and is 
working in the Engineering Laboratory of the University. 

A. L. R. Dztwinski (Graduate) formerly Petroleum 
Laboratory Technician, Technical College, Henlow, is now 
an engineer, London office, de Havilland Propellers Ltd. 

K. FEARNSIDE (Associate Fellow) formerly Technical 
Executive (Research), Smiths Aircraft Instruments Ltd., 
has been appointed Director of Research. 

Wing Cmdr. A. E. Foster (Associate Fellow) formerly 
Group Engineer Officer, Air Headquarters, Malta, is now 
at R.A.F. Honington as Officer Commanding Technical 
Wing. (Correction to October Notices.) 

D. J. Foster (Graduate) formerly an Aerodynamicist, 
Avro Aircraft Ltd., Malton, is now an Engineer, Perform- 
ance Office, Atomic Power Division, English Electric Co. 
Ltd., Whetstone. 

Fit. Lt. NoRMAN GARDNER (Associate Fellow) formerly 
at the R.A.F. Technical College, Henlow, is now Officer 
Commanding Armament Squadron at Tengah, Singapore. 

I. A. B. Gaunt (Associate Fellow) of the General 
Electric Company Ltd., has been appointed Project 
Engineer and Deputy Manager in Tokyo in connection 
with the construction of the Nuclear Power Station for 
the Japan Atomic Power Company. 

A. Netw Gossy (Graduate) formerly with Hunting 
Engineering Ltd., at Luton Airport is now a lecturer in 
the Department of Aeronautical and Mechanical Design 
at Hatfield Technical College. 

L. H. GouLp (Associate Fellow) formerly a Section 
Head, Stress Department, Avro Aircraft Ltd., Malton, is 
now with the Republic Aviation Corp., Farningdale, New 
York as a Specialist Stress Engineer, Research and 
Development Division. 

J. C. GREEN (Graduate) formerly Technical Assistant, 
Stability and Control Section, Vickers-Armstrongs (Air- 
craft) Ltd., Winchester, is now an Engineer III, United 
Kingdom Atomic Energy Authority (Production Group), 
Capenhurst Works. 

G. R. HAMPSHIRE (Associate) has left Aden Airways, 
Aden, and joined Air Links, Gatwick, as Chief Engineer. 

J. C. HARDING (Associate) formerly Chief Flight Engin- 
eer, Airwork Ltd., in the Air Transport Division, is now 
Flight Test Engineer, Handley Page Ltd., Radlett. 

J. F. HARRIMAN (Associate Fellow) formerly with Erics- 
son Telephones Ltd., is now Development Manager, 
Raleigh Industries Ltd., Nottingham. 
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Some Factors Affecting the Choice of Materials 


HEATH 


(Assistant Chie) Structural Engineer, A. V. 


Roe and Co. Ltd.) 


Introduction 

For the purpose of this paper, certain fundamen- 
tal parameters are assumed to be outside the control 
of the structures designer. It is, however, realised 
that the designer is often able to show that small 
changes in one or more of these parameters may have 
far-reaching effects. 

These basic parameters include the aircraft size, 
weight, operating speed and altitude, aerodynamic 
characteristics and load factors. Fig. 1 shows how 
these combine together in various ways to give the 
loads acting on the structure, the temperature environ- 
ment and the stiffness required to prevent aeroelastic 
disturbances. 

At this stage the remaining specified parameter is 
introduced, that of service life. For a particular 
material this parameter, coupled with the loads and 
temperature conditions, establishes the stress levels 
which can be tolerated if undesirable fatigue and 
creep phenomena are to be avoided. 

The structures designer is now in a position to 
review the various materials of which the properties 
are known. It is perhaps not unfitting at this juncture 
to point out that even in 1959 when the design and 
construction of a trans-Atlantic supersonic air liner is 
being seriously considered, tremendous gaps exist in 
our knowledge of material properties at elevated 
temperatures. 

Two other considerations must not be overlooked. 
Firstly, any material which is considered must be 
available at the time and in the form (sheet, extru- 


sions, forgings, and so on) which the design 
necessitates. Secondly, the workshops must be 
prepared to handle the material. This latter point 
will also influence the designer in his choice of the 
type of structure, which may involve the workshops 
in fabrication processes requiring heavy capital 
expenditure on new plant. 

This paper considers how structural efficiency, 
stiffness, creep and fatigue may affect the choice of 
materials which are suitable for structures operating 
at temperatures up to 500°C. The main,concern is 
with currently available materials, i.e. aluminium 
alloys, titanium alloys and steels. Brief reference is 
made to the problems of flight at temperatures up to 
1,000°C. 
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SPECIFIC TENSILE STRENGTH ~INCHES x 10° 


50 400 500 
TEMPERATURE 

FiGcure 2. Specific tensile strength of sheet materials 

(10 hour soaking time). 


° 


2. Structural Efficiency 
2.1. GENERAL 

Structural efficiency may be defined as a strength/ 
weight or stress/density ratio having the units of 
length. For materials it is often quoted as the ratio 
of the ultimate tensile strength to the density and for 
fabricated components (especially compression 
panels) it is usually expressed as the load at failure 
divided by the weight per unit length or width. 

The definition is extended slightly in this paper to 
embrace specific stiffness, i.e. the ratio of the elastic 
modulus to the density. 


2.2. SPECIFIC STRENGTH 

Most structural designers are familiar with curves 
of the type shown in Fig. 2, which depict the reduction 
of specific ultimate tensile strength (U.T.S./density) 
with increasing temperature. Usually several materials 
are represented and for a true comparison each 
material must have been “soaked” for the same period 
prior to testing. 

Curves of this type suffer from several disadvan- 
tages, and because these have not always been fully 
appreciated the curves have often been misused. It 
is thus worthwhile to consider the various shortcom- 
ings of Fig. 2, so that its pitfalls may be avoided. 

Very little of an aircraft structure is designed with 
tensile strength as the sole requirement. More than 
half of the structure has stability in either compression 
or shear as its criterion, and the requirements of stiff- 
ness, fatigue and creep play such an important part in 
the design of the remaining structure that tensile 
strength of itself is of little importance. There may 
be some use for these curves in the design of the short 
life missiles, where creep and fatigue are unimportant, 
but even in these cases the high speeds involved may 
well make stiffness an over-riding criterion. The use 
of curves similar to Fig. 2 without reference to other 
relevant data can therefore be misleading. 


2.3. SPECIFIC STIFFNESS 
For those components which are designed 
primarily by aeroelastic considerations, knowledge of 
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SPECIFIC STIFFNESS INCHES x 10” 
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Ficure 3. Specific stiffness of sheet materials 
(10 hour soaking time). 


the specific stiffness of the various materials can be of 
use. Fig. 3 shows a typical plot of this parameter, 
represented by the ratio of Young’s modulus to 
density, plotted against temperature. It is of interest 
to note that the curves fall in a somewhat. different 
pattern to those of Fig. 2. 


2.4. SPECIFIC STABILITY 

Since so much of the airframe is designed to be 
stable under compression and shear loads, it is con- 
venient to consider a parameter which may be used 
to show the relative merits of various materials in 
withstanding loads of this type. Stability is essentially 
a stiffness problem, and at first it may be thought that 
the specific stiffness values of Fig. 3 are sufficient 
indication of how the various materials will behave 
when built up into, say, skin-stringer panels. 

It is not difficult to show, however, that the 
structural efficiency of an element of a compression 
structure may be written as 


where P critical load per unit width of element 
W weight per unit area of element 


C constant, depending on edge conditions, 
etc. 


b leading dimension (length or width) of 

element 

plasticity reduction factor 

Young’s modulus 

material density 

second moment of area per unit width of 

element 

effective thickness of element, i.e. thick- 

ness of single sheet having same weight. 
The structural efficiency thus depends on _ three 


variables : — 
(a) The structural index (P/ 5), which is a function 
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FiGure 4. Stability parameter of lightly loaded plates 
(10 hour soaking time). 


of the applied loading and the size of the ele- 
ment only. 

(b) The material property parameter (yE)'/,, 
which depends not only on the shape of the 
stress/strain curve but also on the stress level, 
since » varies with stress. The optimum stress 
is also a function of the structural index. 

(c) The non-dimensional shape parameter /*//, 
which depends on the type of structural 
element considered (e.g. skin-stringer panel, 
honeycomb sandwich panel, and so on). The 
optimum value of this parameter is also a 
function of the structural index. 

There is thus only one independent variable, which 
is the structural index. For low values of the index, 
the optimum stress is low and thus »=1. For a given 
type of construction operating at low values of the 
structural index, the parameter E'/,p is thus a guide 
to the efficiency of each material. A curve showing 
typical values of this “specific stability” plotted 
against temperature is given in Fig. 4. The order of 
merit of the various materials is quite different to that 
of the specific stiffness curves of Fig. 3. 

For high values of the structural index, the effici- 
ency is limited by the maximum allowable stress, 
which is taken as the compressive yield or 0:2 per cent 
proof stress, 

The maximum obtainable efficiency is thus o,/p, 
and a plot of this parameter against temperature is 
given in Fig. 5. 

The two “specific stability” parameters discussed 
above are valid only at the low and high ends of the 
structural index scale respectively. The intermediate 
values are often those on which interest is focused, 
and these are considered in the next Section. 


2.5. TYPE OF STRUCTURE 
For compression structures in particular, it is well 
known that different types of constructions exploit the 
material properties to different degrees, depending on 
the value of the structural index. There thus arises 
the condition where a material exhibiting high values 
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of the “specific stability” parameters may yet 
give a very low efficiency if it is made into the 
wrong type of structure. This is particularly true at 
low to medium values of the structural index (values 
found, for example, in low aspect ratio wings) where 
the optimum stress is lower than the yield stress. 

Consider, for example, Fig. 6 which shows the 
relative merits of two similar types of construction, 
honeycomb sandwich and corrugation sandwich, in 
three different materials after short-time soaking at 
250°C. It is obvious that over a very wide range of 
values of the structural index, titanium honeycomb 
sandwich panels are by far the most efficient. But if 
this type of construction involves difficult workshop 
processes (as indeed it does) which preclude its use, 
what is the next most efficient construction? For 
certain values of the structural index it is certainly not 
titanium corrugation sandwich; in fact an aluminium 
corrugation sandwich may be slightly more efficient 
and a steel honeycomb sandwich considerably more 
efficient than the “obvious” alternative. 

It is therefore concluded that the choice of the 
material and the type of structure must be made 
simultaneously, since these two results are mutually 
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FiGure 6. Structural efficiency of sandwich construction at 
250°C. 
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dependent. They are also dependent on the configu- 
ration of the aircraft; in other words the choice of 
structure and material for one aircraft does not 
necessarily apply to another aircraft even though both 
may be operating at the same temperature. 


3. The Effects of Time on Structural Design 
3.1. INTRODUCTION OF TIME-DEPENDENT EFFECTS 

When considering long-life structures, such as 
supersonic transport aircraft which will be expected to 
have a service life of up to 30,000 hours, allowance 
must be made for at least three time-dependent effects. 
These are fatigue, corrosion and creep. 

The latter effect is new to the airframe designer, 
and information relevant to the materials, temperature 
and times of interest to him is scanty. Although 
fatigue and corrosion are familiar problems, little is 
known about the effects of high temperature on these 
phenomena, and even less about their interaction with 
creep. 


3.2. FATIGUE 

Fatigue damage depends as much on detail design 
and working stress level as it does on the basic 
material. Any assessment of the relative merits of the 
resistance of test coupons to fatigue must therefore 
assume that the same standard of detail design can be 
achieved in practice, no matter what type of construc- 
tion is contemplated. 

Test results from fatigue tests on simple specimens 
may be used to give a provisional comparison of the 
S-N curves of various materials. If these curves are 
also plotted as “specific fatigue strength” by dividing 
the alternating stress by the appropriate density (see 
Fig. 7) some idea can be gained of the relative 
efficiencies of the materials. 

These simple tests will presumably be carried out 
under rapid fluctuations of stress at the required tem- 
perature. Long term tests should be performed in 
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Figure 7. Specific fatigue strength. 
Alternating stress = mean stress. 
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“real time” so that the effect of soaking time can be 
evaluated. In order to include the usual factor for 
scatter, periods of up to 100,000 hours (12 years) would 
be required for transport aircraft. 

Within this fantastic time scale more useful results 
could be obtained by testing portions of actual struc- 
tures under realistic conditions including transient 
heating and cooling, take-off, flight and landing loads 
and periods of rest. A supersonic transport aircraft 
may well be accumulating flying hours at almost the 
same rate as its fatigue test specimens—a very good 
reason for using “fail-safe” structures for such aircraft. 


3.3. CREEP 

A criterion for evaluating the creep resistance of a 
particular material is difficult to define, since creep 
“damage” can take several forms. These include 
creep rupture, excessive creep deformation (which may 
lead to control and aeroelastic problems), and creep 
buckling of compression members. 

Creep is dependent on both time and temperature 
and it is desirable to link these two together into a 
single parameter to enable direct comparisons of test 
results to be made. Larson and Miller’ have proposed 
a relationship between tensile rupture time and tem- 
perature for a given constant stress in the form 


T (a, + a, log t)= constant 
where T absolute temperature 
a,,a, constants 
t time to rupture. 


Recent tests by the N.A.S.A. have shown that a 
similar relationship holds good for the creep buckling 
of flat plates in compression. 

A plot of rupture or buckling stress against the 
Larson-Miller parameter thus forms one simple guide 
to the creep resistance of materials. Such a curve is 
shown in Fig. 8, using specific stress for the vertical 
scale. 

Similar curves may also be drawn for creep in 
bearing and shear so that the deformation of typical 
joints may be studied. Fig. 9 shows such a joint, 


T: TEMPERATURE : °K 
t: TIME HOURS 


SPECIFIC CREEP STRENGTH INCHES » 
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Ficure 8. Specific creep strength of sheet materials 
(one per cent creep deformation). 
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T (21+ Logt) x 16° 
T = temperature °K. t=time in hours 
FiGuRE 9. Creep strength of riveted joints. Matl: 2024-T3 
(8 per cent rivet deformation) 


together with a curve of stress plotted against time 
required to produce either rupture or 8 per cent (of 
rivet diameter) deformation. This type of creep 
deformation can give rise to serious structural 
troubles, such as looseness of connections and exces- 
sive displacements of components. 

The various types of stress—tension, bearing and 
shear—which are present in such a joint must be care- 
fully balanced so that all produce similar rates of 
creep. 

The Larson-Miller relationship holds good only if 
no metallurgical transformations (e.g. precipitation) 
occur during the soaking period. Both extrapolation 
and interpolation may be misleading if such changes 
are likely. There is no reliable substitute for creep 
tests at the temperature and for the time under 
consideration. 

Like fatigue, the full effects of creep can be assessed 
only by reference to the load spectrum for the 
particular aircraft in question. For transport aircraft 
operating at around a Mach number of 2, studies 
show that creep may well entail a restriction on the 
use of aluminium alloys, since short bursts of over- 
speeding and excess g can produce large creep 
deformations. 

Unlike fatigue, however, creep “damage” is 
partially recoverable when the applied loads are re- 
moved. The extent of this recovery is a matter for 
controversy, and at this time it is not proposed to 
discuss this phenomenon. 


CHOICE 


OF MATERIALS 


3.4. INTERACTION OF CREEP AND FATIGUE 

Evidence regarding the effects of fatigue loading 
on creep are conflicting; some results suggest a 
deterioration in creep properties due to the superim- 
position of an alternating stress while others suggest 
an improvement. It is clear that, at present, these two 
time-dependent phenomena must considered 
separately. 


4. New Materials 

As the design temperature increases the well-known 
alloys of aluminium, titanium and steel gradually 
disappear from the field of usefulness, and as tempera- 
tures of the order of 1,000°C are approached new 
materials must be found. 

Molybdenum shows considerable promise but it 
oxidises easily and the oxide is volatile, so that the 
metal burns away. Chromium plating appears to pro- 
vide a reasonable surface protection against this 
catastrophe. At 1,000°C, the 0-5 per cent titanium 
alloy has a specific stiffness similar to that of the 
aluminium alloys at room temperature; its “specific 
stability” values are considerably lower, however, and 
creep is likely to be the design criterion for quite short 
soaking times. 

Other possible materials include niobium and the 
cermets and carbides. So little information is available 
on these materials that comment is valueless; one can 
only bear them in mind as “possibles” until research 
reveals their true place in the scheme of things. 


5. Conclusion 

Simple curves relating the specific strength, stiffness 
and stability of materials to temperature and soaking 
time have a limited usefulness. The information they 
present is often conflicting and the designer may be 
misled if undue emphasis is placed on one particular 
aspect. 

The final choice of material is strongly linked to the 
type of structure, and both of these may be considerably 
influenced by the availability of the material and the 
workshop processes entailed. 

Much more research is needed into the time-depen- 
dent effects of creep and fatigue. Over-simplified 
laboratory tests may give optimistic results when com- 
pared with realistic programmed tests on aircraft 
components, and the choice of any particular material 
will be ultimately justified only by full operational 
flying. 
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Design of Transparencies 


by 
J. S. PRZEMIENIECKI, Ph.D., B.Sc., A.F.LA.S., A.F.R.Ae-.S. 
(Bristol Aircraft Limited) 


SumMary :— A study of engineering problems encountered in the design of transparencies for 
supersonic aircraft is presented with a special emphasis on the effects of kinetic heating. The 
paper is written primarily with optical transparencies in mind but many of the general 
principles discussed are equally applicable to the design of radomes for elevated temperature 
applications. 

The first parts of the paper consider briefly the methods of calculating temperature 
distribution in transparent panels subjected to kinetic heating (or cooling) and this leads 
naturally to discussion of the effects of the non-uniform temperature distribution, differential 
expansion and the edge restraints preventing thermal deformations. Methods of calculating 
thermal stresses are suggested through the concept of the “ self-equilibrating ” temperature 
distribution; the main attention, however, being concentrated on the derivation of simple 
engineering formulae rather than on the mathematical rigour of the solution. 

The choice of design cases for combined pressure and thermal loading is considered and 
it is shown that, in general, deceleration manoeuvres combined with outward pressure con- 
stitute the design criteria for transparencies. The effects of various parameters on the 
magnitude of thermal stresses are analysed; these include the effects of altitude, temperature. 
deceleration and material properties. It is shown that, for given material and panel geometry, 
the thermal stresses depend mainly on the following parameters: total temperature change, 
acceleration (or deceleration), and the Biot number based on the panel thickness. 

Practical design aspects are also considered, such as choice of materials, strength, panel 
geometry, mounting of panels and the fail-safe design. Finally, testing of transparencies 
under pressure and thermal loading is briefly discussed, together with possible methods of 
strain measurements, and a tentative proving procedure for transparencies is suggested. 


1. Introduction meteoritic dust and adequate protection from radiation 
hazards. 

The paper is written primarily with optical trans- 
parencies in mind—with a special emphasis on the 
design of flat panels—but many of the general 
principles discussed are equally applicable to the design 
of curved panels or radomes. 

A number of controversial points will most certainly 
be raised and if the material presented here is sufficient 
to spark professional argument and lead eventually to 
a clearer understanding of the design problems of 
transparencies for elevated temperatures, the main 
purpose of this paper will be realised. 


In the design of supersonic aircraft serious con- 
sideration must be given to the elevated temperature 
effects on the transparent materials for window panels 
and windscreens. The temperature limitation of the 
present organic plastics imposes severe restrictions on 
the use of these materials above 80°C. For higher 
temperatures, glass or fused silica are the alternative 
materials that may be used; their use, however, creates 
many engineering problems, largely unexplored, par- 
ticularly those concerned with thermal effects due to 
aerodynamic heating and cooling and the safe working 
stress in the material. 

As its principal aim, this paper attempts to present 
a study of engineering problems encountered in the 
design of transparencies for supersonic aircraft with a x,y, 
special emphasis on the effects resulting from kinetic 


NOTATION 


N 


rectangular co-ordinates (z is measured from 
the middle surface) 


heating. Some of the problems discussed are naturally . — a 
common to all types of structures for high temperatures, L 
others are peculiar only to transparencies treated as time 
structural components. E Young’s modulus 
The design of transparencies for supersonic aircraft 2 coefficient of thermal expansion 
poses a difficult and challenging problem to the struc- ¥ Poisson’s ratio 
tural engineer. The material selected must satisfy k thermal conductivity 
many requirements, some of which are conflicting: it p density 
must be capable of withstanding pressures specific heat 
arising from aerodynamic loading and cabin pressurisa- : = diffusivity =k / pc 
tion, Tapid aerodynamic heating or cooling, heat from 
de-misting and de-icing systems, and finally impact and p. density of air 
erosion, especially from rain and hail. For transparen- V velocity 
cies designed for space travel additional requirements @¢ Stefan-Boltzmann constant in the radiation 


may be added: resistance to erosion caused by law 


4 
§ 
a 


J. S. PRZEMIENIECKI 


D=Ed*/12(1 

f., f, stresses in the x- and y-directions respectively 
s tensile strength (or modulus of rupture) 
w panel deflection in the z-direction 
pressure 

T., adiabatic wall temperature 

T, surface temperature 

T, ambient temperature 

T, initial temperature 

T, steady state temperature 

T temperature distribution through the panel 
thickness = T (z, 1) 


T,, mean temperature, defined by equation (14) 
AT temperature gradient, defined by equation 
(15) 


T*  self-equilibrating temperature distribution, 
defined by equations (12) and (13) 
heat input 
h convective heat transfer coefficient 
B Biot number = /d/k 
M Mach number 
Re Reynolds number 
St Stanton number 
R, boundary resistance (analogue circuit) 
applied voltage proportional to (ana- 
logue circuit) 
F_ Aijry stress function 
w, the jth root of equation (11) 
}-, B, parameters defined by equation (19) 
@ parameter defined by equation (23) 
Yn=mz=b/2a 


Laplace operator— 
Cx 


v4 


biharmonic operator 
CX 


2. Design Requirements 

The main function of optical transparencies is to 
provide an adequate field of vision to enable the pilot 
to control the aircraft and while it is possible to 
imagine an aircraft flying optically blind, some restric- 
ted field of vision must always be provided for take-off 
and landing or for emergency. 

In the early days of aviation no optical transparen- 
cies were needed except perhaps for a pair of goggles 
worn by the pilot who usually had an excellent, 
unobstructed view, in practically all directions. As 


Windscreen structure of a large subsonic aircraft 


FiGuURE 1. 


KINETIC HEATING—DESIGN 


OF TRANSPARENCIES 


aircraft speeds increased and pressurised fuselages were 
introduced, the need for an enclosed cockpit became 
essential and this led to the present form of windscreen 
assemblies and cabin windows. A typical windscreen 
assembly of a large subsonic aircraft is shown in Fig. 
1; here large quadrilateral transparency panels are used 
and the field of view is generally good. On smaller 
aircraft the bubble type canopy became a characteristic 
feature, as illustrated in Fig. 2. For supersonic aircraft, 
operating at elevated temperatures, transparencies con- 
sisting of relatively small flat panels supported in a 
metal framework must be used. Thus some restriction 
of the field of vision must inevitably be expected. All 
panels would be duplicated and have bevelled edge type 
mounting to give a flush external finish and to avoid 
unnecessary edge restraints. Naturally it is also 
possible to use curved panels, and it has already been 
tried in one British supersonic aircraft, but the inherent 
difficulties in the manufacture, and the high thermal 
stresses arising from the edge restraints on the curved 
boundaries, make this type of construction rather an 
unattractive proposition 

The requirements for passenger windows are 
altogether different. Windows are provided mainly for 
psychological reasons and it is therefore not unreason- 
able to assume that future airlines and passengers could 
be educated to accept supersonic air liners without 
any windows in the passenger compartments. This 
would undoubtedly remove, at the same time, many 
difficult problems in the design of pressurised fuselages 
with numerous cut-outs, and a greater structural safety 
and economy would be ensured. Of course, a splendid 
view from a passenger seat of a supersonic air liner 
flying at a Mach number of 2 or 3 may present an 
additional attraction for the potential passengers; but 
surely, this small asset could be easily sacrificed for 
the sake of safety and economy of the structure. A 
closed circuit television camera could provide a reason- 
able view to reassure any passengers suffering from 
claustrophobia; or they could use slower aircraft with 
conventional windows. 

Public reactions to windowless air liners could be 
tested now so that any decisions affecting window trans- 
parencies for future generations of supersonic air liners 


FiGuRE 2. Typical bubble type canopy. 
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would be formulated accordingly. However this could 
only be done with the co-operation of a scheduled 
airline willing to undertake this particular investigation. 


3. Aerodynamic Heating or Cooling of 
Transparencies 

The temperature of a transparency can be affected 
by any of the following: 

(i) Change in ambient temperature (i.e. change in 

altitude). 

(ii) Kinetic heating or cooling. 

(iii) Solar radiation. 
(iv) De-icing or de-misting methods involving heat. 

(v) Cabin heating or cooling. 

It appears, however, that if a supersonic aircraft is 
designed to operate at Mach numbers greater than two 
only the first two causes are of primary importance in 
the calculation of temperature distribution for design 
purposes. Attention must be paid also to the proper 
design of de-icing and de-misting systems to avoid any 
likelihood of thermal conditions which are worse than 
those arising from kinetic heating and changes in 
altitude. 

To determine the effects of kinetic heating and rapid 
changes in altitude on a transparent panel, a knowledge 
of the temperature distribution as a function of space 
and time is required for the panel itself and for the 
surrounding structure. Since the thermal conductivi- 
ties of materials proposed for elevated temperature 
transparencies are relatively low, changes in speed and, 
or altitude produce transient temperature gradients 
through the thickness of a panel. Therefore all trans- 
parent panels may be referred to as “thermally thick” 
panels for which the heat flow through the thickness 
must be considered. 


3.1. BASIC EQUATIONS FOR HEAT TRANSFER 

If an aircraft is accelerated to a high speed, then 
all external surfaces are heated aerodynamically and 
the heat input Q per unit surface area is given by 


(1) 


where the first term on the right hand side represents 
forced convection heating (kinetic heating) and the 
second term is the radiation heat loss from the surface 
at the temperature T,,. 

The so-called adiabatic wall temperature 7... is a 
temperature which would be attained at a surface in 
equilibrium without any radiation losses; it depends 
mainly on the local ambient temperature 7,. the Mach 
number M, and the type of boundary layer. It can be 
calculated simply from: 


T.,.=T, 1 +0167 . (2) 
for laminar boundary layer 


for turbulent boundary layer 


The convective heat transfer coefficient fA will 
be discussed in greater detail in Section 3.2. 

The relative importance of the radiation term in- 
creases with altitude and temperature. However, since 
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the critical design cases for thermal loading on trans- 
parencies usually represent deceleration manoeuvres at 
low altitudes, as will be shown later, the radiation term 
can be neglected, at least as a first approximation. 

The heat input density due to kinetic heating (or 
cooling) of a panel is approximately uniform and there- 
fore conduction of heat in the panel can be assumed to 
proceed mainly in the direction normal to the surface, 
say the z-direction, and thus the problem reduces to 
finding the solution of the heat conduction equation in 
one dimension, i.e., 


0 oT oT 


If, however, thermal conductivity is independent of tem- 
perature, the above equation reduces to the well-known 
diffusion type equation 


10T 
(9) 
ct 


where «=k/pce is the thermal diffusivity. 

To obtain a complete solution to equation (5) 
boundary conditions are required at the outer and inner 
surface. The boundary condition at the outer surface, 
which can be either heated or cooled depending on 
whether 7, is greater or less than T,, can be expressed 
as 


oT 


MOT... z=d/2 . (6) 


where for convenience the z-direction is measured 
from the middle surface towards the outer surface. In 
this boundary condition the rate of heat absorption (or 
loss) by the material at the outer surface is equated to 
the rate of heat gain (or loss) by forced convection. 
For the second boundary condition, one can assume 
that either the surface is insulated, in which case 


oT 

—<=0; z=-d/2 . . 
ez 

or, if the panel is in perfect thermal contact with 

another panel, 


where suffixes | and 2 refer to the two panels in contact. 

This latter condition may be used if the transparent 
panel is of laminated type, such as a laminated glass 
panel with an interlayer of vinyl or other similar 
material. Sometimes it may be necessary to provide 
cooling of the inner surface. The boundary condition 
would then express the balance between the heat loss 
from the surface and the heat extracted by the cooling 
air. In the general discussion that will follow, however, 
only the boundary condition of zero heat transfer, i.e. 
equation (7), will be used; this forms a suitable design 
assumption and leads generally to conservative designs 
(see Section 7.1.). 
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Figure 3. Variation of local heat transfer coefficients for 
laminar and turbulent boundary layers in incompressible flow 


3.2. CONVECTIVE HEAT TRANSFER COEFFICIENT 
The convective heat transfer coefficient /: is deter- 
mined from the following equation 


(9) 


where Sr is the Stanton number, c, is the specific heat 
of air at constant pressure, p, is the air density and V is 
the velocity of air; all these quantities being evaluated 
for the aerodynamic conditions just outside the boun- 
dary layer. The product c,p,V for any given flight 
manoeuvre is a known function of time since it depends 
only on the variation of altitude and speed and c, can 
be taken as a constant. The Stanton number Sr varies 
not only with the Reynolds and Mach numbers, which 
are known functions of time, but also with the ratio of 
external surface temperature to ambient temperature, 
(T./T,). This latter dependence is more pronounced 
if the boundary layer is turbulent; it also introduces 
non-linearity into the boundary condition at the exter- 
nal surface. 

A convenient way of determining the heat transfer 
coefficient is to calculate first its incompressible flow 
value A, which depends only on the Reynolds number 
Re and then to calculate its compressible flow value h 
from the variation of h/h, with the Mach number M 
and the temperature ratio T,/T,."' 

The variation of the incompressible heat transfer 
coefficients has been plotted in Fig. 3, which shows that 
h, is a decreasing function of Re. For laminar boundary 
layers it would be sufficiently accurate to use /, instead 
of the compressible value /; this would give values of h 
on the high side, but it should not be more than 10 per 
cent too high for Mach numbers up to 4. For turbulent 
boundary layers compressible flow formulae must be 
used since / decreases appreciably with an increase in 
the Mach number or the temperature ratio 7,/7,, as 
shown in Fig. 4. A perusal of Fig. 4 indicates clearly 
that only for small variations in M and T,/T,, the 
assumption of a constant heat transfer coefficient may 
be justified. In general, however, variable heat transfer 
coefficients should be used for calculations on 
transparencies. 

Since it may not always be possible to maintain 
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FiGure 4. Variation of local heat transfer coefficient for 
turbulent boundary layer in compressible flow (h,=incom- 
pressible value). 


laminar boundary layers over transparencies, turbulent 
flow formulae should be used for design purposes. It 
should be mentioned, however, that the main difficulty 
in calculating h, (see Fig. 3) is the determination of the 
effective Reynolds number Re. The equation used for 
producing Fig. 3 is essentially a flat plate formula for 
which Re is based on the distance from the transition 
point which may be particularly difficult to determine 
on transparencies. To overcome this difficulty the 
effective transition point for design purposes could be 
assumed to be one to two feet forward from the edge 
of the panel; this would probably ensure that the 
estimate for the heat transfer coefficient is conservative. 
Experimental data in this field would be extremely 
valuable. 


3.3. METHODS OF SOLUTION OF THE HEAT 
CONDUCTION EQUATION 


The following methods can be used for calculating 
transient temperature distributions in thermally thick 
panels : 

(i) Analytical methods 

(ji) Numerical methods. 

(iii) Semi-numerical methods. 

(iv) Analogues 
Each of these methods will be discussed briefly and 
possible methods of solutions will be suggested for 
calculating transient temperatures in transparencies. 

Several objections may be made to the use of the 
analytical methods. Firstly, solutions may not be suffi- 
ciently realistic due to the various simplifications which 
are necessary to get any solution at all; secondly, the 
analytical expressions are usually extremely complica- 
ted for practical applications; and thirdly, the analytical 
solutions usually converge poorly near t=0 unless 
special asymptotic solutions are used for small values 
of time. 

The arbitrary variation of A and T,, with time and 
also the dependence of / on temperature, rules out any 
possibility of a successful analytical solution for the 
temperature distribution through the thickness of a 
transparency panel. It should, however, be emphasised 
that in many cases sufficient engineering accuracy may 
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be achieved with constant average values of the heat 
transfer coefficient A for which analytical solutions are 
available. In particular, these methods could be used 
if only a study of the effects of various parameters is 
required in the preliminary design rather than an 
accurate quantitative study. 

The temperature distribution through a thermally 
thick panel for suddenly increased (or decreased) adia- 
batic wall temperature from 7, to 7, and for constant 
heat transfer coefficient is given by: 


T -—T, sin COS w ( kt 
(w, + sin w; COS w’,) 
(10) 


where w, is the jth positive root of the transcendental 
equation 


w, tan w;=hAd/k. (11) 


The above formula for temperature can also be applied 
to the case of arbitrary variation of 7, by using the 
Duhamel’s integral. 

The numerical methods are based essentially on 
step by step integration of the equation of heat conduc- 
tion. These methods have the additional advantage 
over analytical methods that the variation of thermal 
properties of the material with temperature and the 
temperature dependence of the heat transfer coefficient 
can be included. Similarly any radiation losses and 
other non-linear effects, e.g. free convection, can be 
included in the solution. Numerical methods applied 
to heat transfer problems are, in general, time-consum- 
ing and require relatively small intervals of time. 
However, they can be programmed for automatic com- 
putation and several such programmes for digital 
computers are already in existence. 

Semi-numerical methods are basically repetitive 
applications of the exact analytical solution to the 
response to triangular variation of T,, over a relatively 
large time interval. This is possible because of the 
linear character of equations (5), (6) and (7). In these 
methods, the heat transfer coefficient A is usually 
approximated as a step function. A semi-numerical 
method has been described by Hill’ and the design 
charts for the application of this method to thermal 
stress problems in thermally thick plates have been pre- 
pared at Bristol”®’. The semi-numerical methods are 
very rapid and are suitable for computation on desk 
machines; in general, they are preferable to numerical 
methods *’. 

There are three main types of analogue computers 
which can be used for calculating temperature distribu- 
tions in transparencies: 

(i) Resistance-capacitance network.” 

(ii) General purpose analogue computer (differen- 

tial analyser). 
(iii) Pure resistance network. 
To discuss the relative merits of various types of 
computers would be beyond the scope of the present 
paper; only the basic features of each type will be 
indicated. 
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The basic principles of the resistance-capacitance 
analogue and the differential analyser are similar; the 
main difference is that in the resistance-capacitance 
network actual resistances and capacitances are used in 
each space element while in the differential analyser 
amplifier units are used for that purpose. On the one 
hand the accuracy of the differential analyser is gener- 
ally better than that of the simple resistance-capacitance 
network. On the other hand the resistance-capacitance 
analogue is considerably cheaper to construct. Both 
types use finite differences in space only; time variations 
are continuous and are therefore suitable for automatic 
plotting. 

The pure resistance network has the advantage of 
somewhat greater flexibility and versatility since it can 
handle temperature-dependent coefficients. Its major 
disadvantage arises from the fact that this type of 
analogue uses finite differences not only in space but, 
also, in time, so that unless automatic procedures are 
used to set up the starting conditions for each time 
increment the solution of a problem of any size is a 
lengthy process. 

It appears, however, that the resistance-capacitance 
analogue has certain advantages. In particular, it 
allows continuous plotting of temperature-time varia- 
tion and the resulting thermal stress at any point 
through the thickness of the panel. An analogue com- 
puter of the resistance-capacitance type has been built 
and successfully operated at Bristol Aircraft Limited 
for routine calculations of temperature distributions in 
one- or two-dimensional structural components subjec- 
ted to aerodynamic heating or cooling. The particular 
application of this analogue to transparent panels will 
be discussed in Section 3.4. 


3.4. THE BRISTOL RESISTANCE-CAPACITANCE ANALOGUE 

The resistance-capacitance analogue computer is 
based on the fundamental similarity of the conductive 
heat flow within a homogeneous body to that of the 
charge in a non-inductive electric circuit. The ideal 
analogue with distributed capacitance must be replaced 
by a “lumped” system and in that sense the method is 
similar to the solution of the differential equation of 
heat flow by finite differences except that, unlike any 
numerical method based on finite differences, this 
analogue gives a continuous variation of temperature 
with time. 


z=d/2 m=NUMBER OF ELEMENTS z=-d/2 
HEATED (OR COOLED) INSULATED 
SURFACE SURFACE 
Ry R R | 
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Ficure 5. Analogue circuit for kinetic heating (or cooling) 
of thermally thick panels e.g. transparencies. 
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FiGure 6. Block diagram for kinetic heating analogue 


A typical circuit used at Bristol for solving equation 
(5) with the boundary conditions (6) and (7) is shown 
in Fig. 5. Heat storage effects are represented by 
electrical condensers C, while the electrical resistances 
R simulate the thermal resistance of the panel. The 
input voltage V,,, is made proportional to the absolute 
value of T,,—7, while R, is made inversely propor- 
tional to the convective heat transfer coefficient /. By 
suitable choice of the scale factors this analogue circuit 
can be operated at any convenient voltages and time 
intervals. 

The variable voltage input V,, and the boundary 
resistance R, are generated in an electro-mechanical 
control unit operated by a programmed paper tape. A 
block diagram showing the basic components of the 
Bristol analogue, connected up for temperature 
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Figure 7. Typical variation of heat transfer coefficient h 
with time for accelerated flight. 
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investigation in transparencies, is shown in Fig. 6. The 
correct variation of R,0cl/h (Re, M, T,/T,) to be pro- 
grammed on tape can only be obtained by successive 
approximations, since the variation of 7,/7, with time 
is unknown. This can be achieved by estimating 
approximately variation R, with time, operating the 
analogue while recording the voltage V, (proportional 
to the surface temperature 7,) and then re-evaluating h 
using the variation of 7, obtained on the analogue. In 
general, one such iteration would be quite sufficient for 
engineering purposes. However, to avoid the difficulty 
of the iterative procedure, particularly where a large 
number of different flight programmes is to be investi- 
gated, the following simplified approach has been used. 

Consider a typical variation of A with time for 
either accelerated or decelerated flight at a constant 
altitude by taking (1) T,=T7,,, and (2) T,=T;; these two 
variations are shown in Figs. 7 and 8. The correct 
variation of h must lie between the curves (1) and (2) 
(dotted curve). Thus for the successive approximation 
one could start from either of the two curves. 
Alternatively two programmed tapes could be prepared: 
one for T,=T,, and the other for 7,=7;. The analogue 
computer is then operated using the two tapes 
in turn and all the required voltages are recorded for 
the two variations in 4. Mathematically this means that 
upper and lower limits for the temperature distribution 
are obtained by this procedure; these limits in general 
will differ by not more than 10 per cent which is quite 
sufficient for the preliminary calculation. If necessary, 
a second approximation could be used for a few 
selected cases. 

A single channel recorder was used on the Bristol 
analogue, but more recorders could be connected if 
required. Up to a maximum of 20 separate outputs can 
be taken from the analogue network via cathode 
followers to four summing amplifiers which in turn may 
be coupled or used separately. The final output voltage 
from the summing amplifiers can be either the tempera- 
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Figure 8. Typical variation of heat transfer coefficient h 
with time for decelerated flight. 
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(b) TEMPERATURE GRADIENT AT 


(c) SELF-EQUILIBRATING TEMPERATURE T* 


Ficure 9. Thermal deformations of a flat rectangular panel 
due to T=Tm + AT (z/d)+T™*, all edges unrestrained. 


ture-time variation at a selected point or the summation 


20 

x a,T, where —I<a,<+1 and T, is the tempera- 
nel 

ture of the mth point. The latter variation is particu- 
larly important since it can be made proportional to 
the thermal siress in a transparency. 


4. Thermal Deformations; 
“Self-equilibrating” Temperatures 
The temperature distribution 7(z,?) through the 
thickness of a panel can be separated into linear and 
“self-equilibrating” temperature distributions, namely 


where T,, (f) is the mean temperature, AT is the 
linear temperature gradient and 7*(z,?) is the self- 
equilibrating temperature distribution for which 
d/2 di2 
-d/2 
Thus only the first two terms on the right hand side of 
equation (12) contribute to the overall deformation of 


10 


Ficure 10. Variation of mean temperature 7, with non- 
dimensional time N = xt/d? (B=Biot number). 
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the panel whereas the 7* distribution, being self-equili- 
brating, produces zero resultant deformations. 7, and 
AT terms produce thermal stresses only in the presence 
of edge restraints while stresses due to T*, away from 
the edges, are not affected by the edge restraints and 
they occur even if all edges are completely free. 

Using now equations (12) and (13) it follows that 
the two linear temperature components are given by 
da/2 
Tdz . (14) 
d 


~da/2 


T,, 


d/2 
12 


« A 
and T(t) 


Alternatively, ST could be obtained by equating the 
first moment of the temperature 7 (z.1) to the first 
moment of the AT (1) z/d distribution about the middle 
plane. 

Figure 9 shows the overall thermal deformations on 
a rectangular panel with free edges; it also jllustrates 
well how the thermal stresses arise in a panel as the 
result of edge restraint from the surrounding structure. 
The mean temperature 7,, gives overall expansion in 
the plane of the panel, the temperature gradient AT 
bends the panel into a spherical surface, while the self- 
equilibrating temperature distribution T* produces no 
overall deformation. (However, some small “self-equi- 
librating” deformations are present near the edges.) 

Variation of T,,, and the self-equilibrating 
temperature components at both surfaces with the 
non-dimensional time parameter N = «t/d* and the Biot 
number B= hd/k is shown in Figs. 10 to 13 for a sudden 
change in the adiabatic wall temperature (7, —7,). The 
variation of the mean temperature 7,, is of an 
exponential type approaching asymptotically the equi- 
librium value while both AT and T* components 
increase to some maximum value and then gradually 
decrease asymptotically to zero. Thus the effects 
caused by AT and 7* are of transient nature. The 
range of values of B for typical cases is approximately 
1 to 10. 
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Ficure 11. Variation of temperature gradient AT with non- 
dimensional time N =x«t/d?. 
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Figure 12. Variation of  self-equilibrating temperature 
T* (d/2. 1) with non-dimensional time N = «t/d 


5. General Equations for Stresses and 
Deformations 


The large deflection theory of plates subjected to 
arbitrary pressure and thermal loading indicates that, 
in general, there is coupling between the membrane 
stresses in the middle plane of the plate and the bending 
stresses. The basic differential equations for the 
general case are as follows * ’ 


Oxey Ox" dy 
(16) 
Ed 
ox” CXCV CXCV ox” oy 
D= Ed’ /12 (1 (18) 


Where the x and y co-ordinates refer to any set of 
mutually perpendicular axes in the plane of the plate, 
as shown in Fig. 14. The resulting stresses are then 
determined from the solutions for the Airy stress 
function F and the transverse deflection, » 

In any design calculations, it is usually sufficiently 
accurate to assume that both the pressure and thermal 
loading are constant over the whole panel, so that p, 
T,, and AT then become functions of time only while 
T* is a function of z and /. 

If one or two pairs of edges are restrained in the 
plane of the panel then a strong coupling between the 
membrane stresses and bending stresses may exist. This 
has been analysed extensively by Williams*:*, who 
solved the problem of the simultaneous loading of 
pressure and heat (7,, and AT) on an infinitely long 
strip with the lengthwise edges either clamped or pin- 
ended. Parkes'’” analysed the case of a flat panel, with 
initial eccentricities and with all edges restrained against 
expansion subjected to uniform temperature increase 
a). 

If there are no restraints in the plane of the panel, 
and the deflections are small, the effects due to the 
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Figure 13. Variation of  self-equilibrating temperature 
T* (—d/2, t) with non-dimensional time N = «t/d?. 


pressure p and temperature gradient AT can be con- 
sidered separately. The self-equilibrating temperature 
distribution produces only self-equilibrating stresses 
which can also be considered separately and then 
superimposed on the other stresses. 


5.1. PRESSURE STRESSES 

The pressure loading on windscreen panels is 
affected by two components: 

(i) cabin pressurisation (internal pressure), and 

(ii) local aerodynamic pressures (external pressure). 
The first component is a function of altitude only while 
the second component is affected by the shape of the 
windscreen and canopy, incidence, sideslip, speed, Mach 
number and altitude. In pressurised aircraft, however, 
the internal pressure component is likely to be the 
major one, (usually about 5 to 7 Ib./in.*) so that 
although the external pressure loading is far from being 
uniform, it can be treated as constant when combined 
with the higher internal pressure. 

On forward facing panels, the maximum outward 
pressure will tend to occur at high incidences, especially 
in the transonic region, and at maximum cabin pressure, 
while the maximum inward pressure will occur at low 
or negative incidences at high speeds and cabin 
unpressurised. Sideslip can also affect external pres- 
sure loading, especially in V windscreens. 

On side and upper panels, the maximum outward 
pressure is similar to that in the forward facing panels 
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Ficure 14. Co-ordinate system for rectangular panels. 
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Pressure stresses at the centre of a simply- 
supported rectangular panel. 


Figure 15. 


but sideslip may be more important. Inward pressures 
will generally be small, even if the cabin is unpress- 
urised. 

The solution for calculating stresses in a rectangular 
panel under uniform pressure with all edges simply- 
supported is well known. The maximum stress occurs 
at the centre of the panel and may be expressed as '' 


where 8, and 8, are functions of the aspect ratio a/b; 
these functions are plotted in Fig. 15. 


5.2. THERMAL STRESSES 
5.2.1. Mean Temperature T,,. 

The mean temperature 7,, induces stresses only if 
one or two pairs of edges are restrained against expan- 
sion in the plane of the panel. In such cases 
compressive or tensile stresses are induced if the panel 
is respectively heated or cooled more rapidly than its 
surrounding structure (mounting frame). Therefore, it 
is advisable in designing transparency mounting to 
avoid any restraint on the expansion of the panel. This 
also almost certainly excludes any design in which the 
transparent panel is bolted firmly to the structure, such 
as is the case with the present laminated glass panels 
where vinyl layer with light alloy insert round the 
edges is bolted to the main structure. 

Thermal stresses due to 7,, are proportional to 
2ET,,,./(1—¥) and this could have a dangerously high 
value''*’—in fact over 24 times higher than the stresses 
due to AT and T* distributions. 


5.2.2. Temperature Gradient AT. 

The effect of AT is to produce a uniform curvature 
in the panel and therefore in the absence of any edge 
restraints the panel bends into a spherical surface. 
Because of the restraints along the edges of a rectangu- 
lar panel, which prevent the edges from leaving the 
plane of the supports, the temperature gradient AT 
produces certain reactions along the boundary and 
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bending stresses in the panel. For circular panels, with 
simply-supported edges, these edge reactions and 
stresses are zero since the edges of the panel remain in 
the plane of the support under the action of AT. 

For a rectangular panel clamped at four edges, the 
bending moments due to AT are constant over the 
entire panel, and twisting moments and transverse shear 
forces are equal to zero. The thermal stresses are 
simply given 


-2E (3) (20) 


which means that the maximum thermal stress occur- 
ring at the surface (z= +d/2) is +2EAT/2(1- ») 
For a rectangular panel simply-supported at four 
edges, thermal stresses have been solved by Maul- 
betsch"*’; the thermal stresses depend on the aspect 
ratio of the panel (a/b) but are independent of the 
Poisson’s ratio. These stresses are given by '*' 


* 
d 
(22) 
d 
sin (mzx/a) cosh (nzy/a))mcosh y,, (23) 


The function 7 can vary between 0 and «/4 and thus 
the maximum value of the stress will be 


f. or fy= +42EAT 


The variation of the thermal stresses f, and f, at the 
centre of a rectangular panel (x =(a/2), y=0) with the 
aspect ratio a/b is shown in Fig. 16; the maximum 
stress, +42EAT is only realised as a/b tends to ~, 
and the minimum thermal stress is only obtained on 
square panels (a/b= 1°. 

For any aspect ratio a/b=0° the maximum thermal 


a£(z/d)aT 


(z/d)AT 


0 
+0 20 2°5 3-0 35 
ASPECT RATIO, (a/b) 


Ficure 16. Thermal stresses at the centre of a simply- 
supported rectangular panel due to temperature gradient 7 
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stresses are always reached at the boundary and are 
given by 


” 


The analytical solutions obtained by Maulbetsch’* 
and Timoshenko''®) indicate that the concentrated 
reactions and the twisting moments at the four corners 
tend to infinity. In ductile materials this would 
probably lead to some local yielding resulting in the 
redistribution of the stresses near the corners. Maul- 
betsch''” suggested a possible solution on these lines 
and recently this problem has also been discussed by 
Bijlaard''® in connection with the analysis of a sandwich 
panel with bolt attachments subjected to a temperature 
gradient. With glass panels which exhibit apparently 
perfect elasticity this mathematical infinity in the 
solution for the stresses may give cause for concern 
A possible practical solution to avoid these stress con- 
centrations is to use generous corner radii. 

Other polygonal shapes of panels can also be 
analysed; for example thermal stresses in triangular 
panels have been discussed in Ref. 14. 


5.2.3. Self-Equilibrating Temperature Distribution T* 
The temperature distribution 7* produces a self- 
equilibrating stress distribution which is given by ’’ 


7/2 2 
1/2 7/2 
— (27) 


For design calculations it is sufficiently accurate to 
assume that the stress distribution given by equation 
(27) is uniform over the entire panel. Only at the edges 
of the panel do the stresses vanish but it may be 
expected that they reach their full values within a short 
distance from the edges. 


5.2.4. Order of Magnitude of Thermal Stresses 

Before discussing the combined pressure and thermal 
stresses it may be appropriate to indicate the magnitude 
of the maximum thermal stresses. It has already been 
mentioned that the probable range of the Biot numbers 
is between | and 10 and therefore attention will be 
confined to these two extreme values. Assuming 
Poisson’s ratio to be v=0-24 and using the results of 
Sections 5.2.1. and 5.2.3. and Figs. 10 to 12, the 
magnitudes of the thermal stresses due to various tem- 
perature components can be estimated. These are 
shown in Table I for a rectangular panel with all edges 
simply-supported. 

The table applies, strictly speaking, only to the case 
of a sudden increase in the adiabatic wall temperature 
(7, —T,) and a constant heat transfer coefficient h. It 
nevertheless provides a useful guide to the magnitude 
of the thermal stresses in transparencies; it also shows 
that it is important to avoid edge restraint in the plane 
of the panel—particularly on outer panels—as otherwise 
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TABLE | 
IHERMAL STRESSES IN TRANSPARENCIES 
( 0:24, glass) 


4 pproximate value of thermal stress 
Temperature 


E(T,—T)) 
Component 0 
B=1 | B=10 
rs 1:32 1-32 
Al + 0°16 
T* 
0-12 0:50 


external surface 


thermal stresses of the order of 1-32z2E (7, —7,) can be 
induced in the panel. 


5.3. COMBINED STRESSES 

Assuming that the thermal expansion in the plane 
of the panel is not prevented, the total stresses will 
consist of stresses due to 7*, AT and the pressure 
loading p which may be either inward or outward. The 
peak values of T7* and AT do not coincide and there- 
fore separate stress time variations must be added to 
obtain the maximum thermal stress and the total 
thermal stress is then combined with the pressure stress. 

The design criterion for transparencies is based on 

the maximum tensile stress. Various design cases are 
quoted in Table II showing the signs of the stresses. 

From this table it can be seen immediately that the 

maximum tensile stress may occur in any of the follow- 
ing cases: 

(i) Deceleration; thermal loading combined with 
outward pressure; stress at z=d/2 (external 
surface). 

(ii) Deceleration; thermal loading combined with 
inward pressure; stress at z= —d/2 (internal 
surface). 


(iii) Acceleration; thermal loading combined with 
inward pressure; stress at z= —d/2 (internal 
surface). 

It has been found, however, that case (i), i.e. thermal 
loading due to deceleration of the aircraft and outward 
pressure, usually gives the highest tensile stress, and 
may therefore be used as a design criterion. 


TABLE Il 
SIGNS OF MAXIMUM SURFACE STRESSES 
ve sion 
ve Compression 
The mal Stresses 
{cceleration Deceleration 
Sircss (Heating (Cooling) 
Caused Outer Inner Outer Inner 
by Surface Surface Surface Surface 
1/2 1/2 z=d/2 z=—dj2 
+ + 
Al 
Pressure Stresses 
p 
Outward 
P 4 + 
Inward 


+ 
4 
i 
: 
i 
Le 
al 
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Figure 17. Typical variation of total stress with panel thick- 
ness; rectangular glass panel. 


Since thermal stresses due to T* and AT increase 
with increasing thickness of the panel (see Figs. 11 and 
12) and pressure stresses decrease with increasing thick- 
ness, there will be some optimum thickness for which 
the combined stress becomes a minimum. A typical 
plot showing the variation of the combined stress with 
thickness for a deceleration case is shown in Fig. 17, 
here the optimum thickness giving the minimum stress 
is about 0-5 in. This figure also illustrates the signifi- 
cant feature of thermal stresses, that an increase in the 
amount of structural material does not necessarily lead 
to a reduction in the magnitude of thermal stress. 


6. Temperature Stresses in Laminated 
Panels 

In laminated glass panels, which are often used in 
windscreen assemblies, the middle layer is usually made 
from polyvinyl butyral plastic which has different 
thermal properties from those of the outer glass layers. 
When these layers experience changes in temperature, 
the middle layer expands or contracts at a different rate 
from the outer layers. Temperature stresses are there- 
fore set up in the layers; in particular large stresses 
may exist near the edges. 

A similar problem has been analysed by Aleck '’ 
who studied the distribution of thermal stresses in a 
rectangular beam, clamped along one side, using the 
principle of least work. The solution obtained could 
not be applied directly to the analysis of laminated glass 
panels subjected to temperature changes, because of 
the assumption of rigid clamping and the essentially 
two-dimensional character of the problem considered. 

Further work by Durelli and Tsao’: *'’, who used 
a photoelastic technique, on a three-ply laminated panel 
provided very useful information on the distribution 
of stresses and their magnitude near the edges. The 
experimental results were again only applicable to two- 
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dimensional distribution. The correlation of the 
photoelastic results with the temperature required to 
produce failure in glass at the interface and with the 
location of the point of failure near the edge, seems to 
indicate that two-dimensional thermo-photoelastic 
technique gives a satisfactory approximation. It may 
be interesting to mention that the temperature change 
in these tests was obtained by lowering the temperature 
to — 40°C. 

Additional experimental and theoretical studies in 
this field, particularly of the three-dimensional stress 
field at the corner of three-ply laminated panels, would 
be extremely valuable. 


7. Effects of Design Parameters 
7.1. STEADY FLIGHT 

In steady flight, when equilibrium conditions are 
achieved, there will be no thermal stresses unless any 
differential expansions of the transparency material and 
the surrounding structure are prevented. The only 
problem is that of finding a material with sufficient 
strength at the equilibrium temperature to withstand 
the pressure loads. 

If the duration of the high speed flight is not 
sufficient to produce a uniform temperature distribution 
through the thickness of a transparency, then the effec- 
tive temperature drop due to any _ subsequent 
deceleration will be somewhat smaller, resulting in a 
decrease of the maximum thermal stress; the same ts 
true if the internal surface temperature is kept lower by 
the supply of cooling air. 


7.2. SPEED AND ALTITUDE 

Thermal stresses are proportional to the total 
change in the adiabatic wall temperature which is a 
function of the Mach number and altitude. Variations 
in the Mach number and altitude have an additional 
effect on thermal stresses through the variation in the 
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Ficure 18. Effect of altitude on thermal stress; analogue 
solution. 
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ger 
Figure 19. Maximum thermal stress f, in simply-supported 


rectangular panels; aspect ratio a/h=2, Poisson's ratio 0:24. 


convective heat transfer coefficient / or the variation in 
the Biot number, B=/hd/k. 

To illustrate the effect of altitude on thermal stresses 
in transparencies, a hypothetical case has _ been 
considered in which the aircraft was decelerated with a 
constant deceleration of 2:5g starting at various 
altitudes and flying either level or in a dive. The panel 
analysed was j in. thick, of aspect ratio of about two, 
and with all edges simply-supported; the materials 
considered were glass and fused silica. The results of 
this investigation are shown in Fig. 18; the solutions 
were obtained on the Bristol resistance-capacitance 
analogue using variable adiabatic wall temperature and 
variable heat transfer coefficient. 

Figure 18 illustrates well how thermal stresses could 
be appreciably reduced by decelerating at high altitudes; 
this procedure should be possible on a supersonic civil 
transport aircraft, but not necessarily on a fighter air- 
craft for which no rigid flight programme could be 


defined. 


7.3. ACCELERATION OR DECELERATION 

The variation of the total thermal stress due to AT 
and 7T* at the external surface with acceleration (or 
deceleration) will be illustrated for the case of linearly 
varying adiabatic wall temperature from 7,,=0 to 
(T,—T,) in time WN, (non-dimensional), thereafter 
remaining constant; the heat transfer coefficient is 
assumed to be constant. In Figs. 19 and 20 carpet plots 
of the maximum total thermal stresses f, and f, are 
presented in non-dimensional form for a series of values 
of N, and the Biot number B=/id/k; the stresses 
represented by these plots are for a simply-supported 
rectangular panel of aspect ratio a/b=2 and Poisson’s 
ratio v=0°24 (glass). For fused silica, which has »v 
0:17, the corresponding curves are only slightly lower. 
It should be noted that the stresses are negative for 
positive values of 7,—7,; thus accelerations induce 
compressive (negative) stresses and conversely decelera- 
tions (negative value of 7, —7,) induce tensile (positive) 
stresses. For square panels, f, and f, will both be equal 
to the mean of f, and f, for the rectangular panel. This 
is apparent from the symmetry of equations (21) and 
(22). 
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Ficure 20. Maximum thermal stress f, in simply-supported 
2, Poisson's ratio v=0°24, 


rectangular panels; aspect ratio a/b 


For a given material, panel thickness and flight plan, 
approximate relations between N, =«t,/d* and average 
acceleration (or deceleration) and between B=hd/k and 
the Mach number and altitude can be obtained. This 
enables the non-dimensional results of Figs. 19 and 20 
to be expressed as the variation of stress with altitude, 
Mach number and acceleration (or deceleration). 

It should, however, be emphasised that there may 
be large variations in the Biot number during the 
acceleration (or deceleration) period (see Figs. 7 and 8) 
and therefore care should be exercised in interpreting 
the results of Figs. 19 and 20 which are based on 
constant values of B. If the variation of the Biot num- 
ber B is not too great, average values of B can be used 
and they would give sufficient engineering accuracy. 
For larger variations more exact methods should be 
used (e.g. Ref. 3), perhaps for a few selected cases. 

For convenience, the variation of the Biot number 
per unit thickness (B/d) with the Mach number and 
altitude is shown in Fig. 21 for Reynolds number based 
on a distance L=3 ft. and 7,=T,,. (Note that the 
latter assumption over-estimates the heat transfer co- 
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FiGureE 21. Variation of Biot number with altitude and Mach 
number (Reynolds number is based here on L=3 ft., 7, is 
assumed to be equal to T,,,). 
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efficient, and hence B, for deceleration manoeuvres, as 
shown in Fig. 8.) A rather interesting feature is 
that for M>2 the Biot number is practically constant 
at any given altitude. 


8. Choice of Materials 
To facilitate the choice of materials, a thermal shock 
resistance parameter may be used as a_ possible 
criterion for comparing different materials; it is given 

by the expression 
Ks 


aE (28) 


where s is the tensile strength at elevated temperatures 
(modulus of rupture). The material with a high value 
of this parameter will generally have betier thermal and 
strength characteristics. 

The low values of the product zE increase the 
thermal shock resistance parameter. For various types 
of glasses, however, the Young’s modulus E is approxi- 
mately constant but the coefficient of expansion z can 
vary appreciably. The use of glass with low thermal 
expansion to improve thermal shock resistance is not 
always desirable, if the glass is to be used in the 
toughened form. Glass is toughened by heating it to a 
temperature near its softening point and then rapidly 
chilling both surfaces to produce locked-up strains, 
compressive in the surfaces and tensile in the centre 
plane. Thus any resultant tensile stresses in the surface 
due to the thermal and pressure loading zre smaller and 
consequently the effective strength of glass is increased. 
The amount of toughening depends on the thermal 
expansion and it may not be feasible to increase the 
strength of the low expansion glasses by toughening to 
the extent possible with glasses of higher thermal ex- 
pansion. For comparison, coefficients of thermal 
expansion for various transparent materials and a few 
typical structural materials are given in Fig. 22; it can 
be seen that fused silica has the lowest coefficient of 
expansion (0°54 x 10° in./in./°C), followed by borosili- 
cate glass (about 10-° in./in./°C). 
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FiGure 23. Thermal conductivity of various materials 


High values of thermal diffusivity (or conductivity ) 
always decrease thermal stresses. The thermal conduc- 
tivities of various materials are presented in Fig. 23, 
here again fused silica appears to be a very attractive 
material because of its high thermal conductivity. Its 
thermal conductivity is 0-0032 c.g.s. units, while that of 
soda-lime glass is about 0-0024. The thermal conduc- 
tivity of plastic materials is about one order of 
magnitude smaller. 

An additional advantage of fused silica is its low 
value of the Poisson’s ratio; a value of v=0°17 has been 
quoted but further experimental evidence would be 
required to confirm this. The Poisson’s ratio of 
ordinary soda-lime glass is about 0-24. 

Other properties which must also be considered are: 
light transmission, ease of manufacture, creep resist- 
ance, and resistance to erosion at elevated temperatures. 


8.1. ORGANIC TRANSPARENCIES (PLASTICS) 

Figure 24 shows the variation of tensile strength 
with temperature for several typical transparent 
plastics. Tensile strength is generally high at room 
temperature but decreases rapidly with temperature 
imposing severe limitations on the use of these 
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Ficure 24. Tensile strength of various plastic sheets 
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materials at elevated temperatures. The present tem- 
perature limit may be taken as around 80 to 100°C but 
further development may possibly extend this by about 
sustained loading, fatigue and high rates of creep at 
elevated temperatures would always make plastic trans- 
parencies a somewhat doubtful proposition. Alterna- 
tive materials for elevated temperatures are inorganic 
transparencies and since these offer much greater scope, 
they will be discussed at greater length. 


8.2. INORGANIC TRANSPARENCIES (GLASS OR 
FUSED SILICA) 

The strength of glass varies widely with its chemical 
composition, but even more so with the surface con- 
dition. rate of loading, degree of toughening and the 
presence of moisture at the surface. 

The low strength of glass is usually explained by 
Griffith’s concept of surface flaws which act as 
stress raisers so that the stress at the base of the crack 
is much greater than the average stress on the surface. 
The size of the cracks is affected mainly by the moisture 
and this is primarily responsible for the delayed failure 
of glass (static fatigue). This theory is generally sup- 
ported by experiments and a mechanism of dissolution 
has been proposed’** in which alkali ion self-diffusion 
controls the initial steps of water corrosion leading to 
breakdown of the glass structure; this in turn permits 
any existing flaws to grow to critical dimensions and 
bring about spontaneous crack propagation and failure. 
The rate of this reaction is influenced strongly by 


temperature. 
Figure 25 shows the experimental results of delayed 
failure tests on soda-lime glass rods**’. From —50°C 


to 150°C the maximum tensile stress decreases with the 
increase in temperature and the time to failure. Above 
150°C the failing tensile stress increases slightly with 
temperature and is practically independent of the time 
of loading, except perhaps for very short time loading 
when there is a slight increase in strength. In a 
similar study Preston**’ found that the minimum 
strength occurred at about 200°C. This peculiar be- 
haviour is generally true of all glasses including fused 
silica. 

The increase in strength at higher temperatures is 
attributed to the decreasing effects of moisture which is 
driven off from the surface. Furthermore, high tem- 
peratures produce sufficient movement of atoms to heal 
the surface flaws and this effect probably outweighs the 
influence of thermal motion, which in itself lowers the 
strength of the glass. Treatment of the surface with 
5 per cent hydrofluoric acid may also tend to heal the 
flaws and increases the strength; this effect, however, 
is only temporary. 

All strength calculations should be based on the 
delayed strength of glass at a temperature corresponding 
to the minimum failing stress (about 150-200°C). 
Sometimes in the absence of any other information, 
short time strength is used with a factor of 2} in 
addition to any factor to cover scatter in strength 
between specimens. 
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FiGurE 25. Delayed failure curves for soda-lime glass. 


8.2.1. Annealed Glass 

The mechanical and optical properties of annealed 
plate glass are sufficiently well-known not to be des- 
cribed in detail. Coefficient of thermal expansion is 
about 8-7x10-° in,/in./°C and this results in high 
thermal stresses. The accepted maximum allowable 
working (tensile) stress is about 800 to 1,200 Ib./in.’. 
In general, the strengih of annealed glass is too low 
to consider it for transparencies 


8.2.2. Toughened Glass 

The toughening process (sometimes called temper- 
ing), increases the tensile strength of glass. The 
generally accepted allowable working in stress is about 
4,500 lIb./in.* for ordinary soda-lime glass in the 
toughened form. 

Borosilicate glasses, of which the most familiar type 
is Pyrex, have a lower coefficient of thermal expansion 
compared with the soda-lime glass (see Fig. 22). As a 
result of this they offer good resistance to kinetic 
heating, but have lower strength in tempered form 
because the lower coefficient of expansion precludes a 
high degree of toughening. Whether this will result in 
a useful material depends on the relative changes in 
both thermal and allowable stresses and further investi- 
gation is obviously needed. The optical properties of 
borosilicate glasses are usually inferior to those of the 
soda-lime glasses 

Aluminosilicate glasses have been suggested as a 
likely material but too little is known yet to evaluate 
their potentialities. 


8.2.3. Fused Silica 

The most attractive material for high temperature 
transparencies appears to be fused silica. It has excep- 
tionally low coefficient of thermal expansion (about 7, 
that of soda-lime glass) and its thermal conductivity is 
about 1:5 times greater than that of ordinary glass; 
both these factors tend to reduce thermal stresses 
appreciably These important factors are easily 
demonstrated by the frequent use of silica ware in 
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Ficure 26. E for glass at elevated temperatures. 


chemical laboratories, in conditions which demand 
heating to red heat followed by immersion in cold 
water. 

The mechanical strength is less well established for 
sheet material and further investigation would be 
required. The working stress would probably be 
slightly higher than for annealed glass, but bearing in 
mind that the thermal stresses would generally be about 
1/20 of those in a similar soda-lime glass panel 
(combined effect of low coefficient of expansion and 
high thermal conductivity), it becomes obvious that 
fused silica is a material of considerable value. This 
is also demonstrated in Fig. 18 where a comparison of 
thermal stresses in soda-lime glass and fused silica is 
made; the thermal stresses in fused silica are reduced 
to such an extent that they are small, even compared 
with the lower allowable stress, and it should therefore 
be possible to ensure adequate safety with much more 
certainty. 

Unlike other types of materials, the Young's 
modulus of fused silica increases slightly with tempera- 
ture (see Fig. 26). 


9. Practical Design Considerations 
9.1. MOUNTING OF PANELS 

Present day designs of windscreen panels use 
laminated glass and the methods of mounting glass 
panels are dependent on the plastic interlayer being 
extended beyond the edges and providing a flange 
which can be bolted to the structure (see Fig. 27). This 
feature will be absent in high temperature applications, 
and alternative methods will be required without the 
necessity for drilling the transparent panels. Bevelled 
edge type mounting may be necessary, as shown also 
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FiGuRE 27. Types of edge mounting on transparencies. 


in Fig. 27, and this may also have the additional 
advantage of flush external surfaces. Furthermore, 
edgewise restraints would not be excessive. 

The frames for glass panels should be made as stiff 
as possible, to avoid straining of the glass due to 
structural deformations under inertia and aerodynamic 
loading. 

9.2. FAIL-SAFE DESIGN 

The fail-safe philosophy can be adopted in the 
design of transparencies**'; an obvious solution is the 
use of a double panel with one panel unloaded but 
ready to take up the full pressure if the other fails. For 
example, the interspace between the two panels can be 
vented to cabin pressure and thus normally only the 
outer panel is subjected to full pressure and thermal 
loading. In an alternative design, the interspace could 
be vented to the air outside and the inner panel would 
then be subjected to cabin pressurisation, with the 
outer panel taking the thermal loading. 

The fail-safe design has already been demonstrated 
in tests on double panel transparencies. In Fig. 28 a 
double panel window is shown afier a test in which the 


FiGure 28. Fail-safe test on a window panel. 
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Membrane stresses in vinyl interlayer 
to pressure loading after initial failure of glass layers 


Figure 29 reacting 


inner panel was loaded under 8:3 Ib./in.* and then 
fractured instantaneously by punching out a scribed 
disc and thus subjecting the outer panel to the full 
pressure suddenly. There was no failure or damage to 
the outer panel, which continued to hold pressure satis- 
factorily. The window panels were made of Perspex 
and were intended for low temperature application. 
This nevertheless shows the usefulness of the principle. 

Naturally designs incorporating duplication as 
safety measure tend to be heavier than the single panel 
designs but in view of the uncertainty about the 
strength, edge restraint, thermal and pressure loading, 
it may still be better to duplicate the panels. 

The use of a vinyl interlayer in laminated glass 
panels, designed to take the membrane stresses, has the 
additional advantage in that, should both sheets of glass 
fail, the vinyl is still capable of carrying the pressure 
loading. This is manifestly demonstrated in Fig. 29 
showing a test on an initially flat rectangular panel in 
which the vinyl is holding a pressure of eight times the 
working pressure after both glass sheets had failed. 

Laminated glass for inner panels could be used if 
cooling air could be supplied to the interspace to keep 
the temperature of the vinyl at a reasonable level. Such 
a scheme for a double panel window is shown in 
Fig. 30. 


10. Acceptance Tests 
10.1. TENTATIVE PROVING PROCEDURI 

Realistic tests on transparencies will require artifi- 
cial cooling of the external surface of a panel. The 
artificial cooling can be obtained by discharging a fluid 
coolant and controlling the rate of discharge, or alterna- 
tively the flow rate could be maintained constant while 
the rate of cooling experienced by the surface could be 
controlled by radiation from a heating source 

The factors required on the transparencies are so 
high relative to the factors on the rest of the structure 
that the transparencies in representative mountings 
could not achieve their factored loads before failures 
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occurred elsewhere. The use of unrepresentative 
mountings, to achieve the required factored loading, is 
a dangerous practice since such mountings would 
seriously modify the load reactions. As an example of 
what is involved, a case has been known in which a 
typical panel mounted in a representative structure had 
an ultimate strength of less than half of that achieved in 
a rigid test rig. 

Although factored pressure loads could be applied 
with non-representative mounting, factoring of thermal 
transient stresses is not easily achieved, and therefore 
a different proving procedure must be employed. 
Assuming that a suitable technique for measuring panel 
stresses can be developed the proving technique might 
be as follows: 

(i) Determine the stresses due to unfactored pres- 

sure loading. 

(i) Determine the maximum thermal stresses under 
representative temperature drop and represen- 
tative rate of cooling. (If no suitable experi- 
mental technique can be developed, theoretical 
thermal stresses would be used for this stage). 

(iii) If possible the maximum stresses should also 
be determined for the combined loading (i) and 
(ii), to check any coupling effects. 

(iv) Add the stresses from (i) and (ii), multiply by 
the appropriate factor and then determine the 
reserve factor by comparing the resulting 
stresses with the maximum allowable stress 
representing a statistically “safe” stress divided 
by the factor of, say, 5 or 6. 

This procedure depends on the development of a 


successful method of strain measurements in trans- 
parencies. It would be outside the scope of this paper 
to discuss various possible experimental methods. 


There is, however, a new photoelastic technique’, 
developed initially in France, which could be applied 
to transparencies. This technique is usually known by 
the trade name “Photo Stress”, and the necessary 
equipment is already being manufactured by an 
American company. 


10.2. “PHOTO STRESS” TECHNIQUI 

“Photo Stress” is a photoelastic technique in which 
the actual full scale structure to be analysed is coated 
with a special clear plastic that becomes doubly refrac- 
tive when stressed. Coating can be obtained by means 
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of liquid plastic or a plastic sheet; the latter can be 
either clear or metallised on one surface. The material 
for coating allows measurements to be taken in the 
range —55°C to + 135°C. 

The method is based on observation and measure- 
ment of coloured fringes by a special electronic photo- 
stress meter. The minimum measured strain, claimed 
by the manufacturer, would be more than adequate for 
the application of this method to transparency testing. 
The measurements of thermal strain would probably 
have to be taken on the internal surface since the plastic 
coating would affect the heat transfer rates unless its 
effect could be compensated by an increased rate of 
heating (or cooling). 

At the time of writing this paper, there has been no 
known actual practical experience with the “ Photo 
Stress ” method in this country. It is hoped, however, 
that this experimental technique will be thoroughly 
explored and, if successful, it will eventually be 
applied to tests on complete transparency assemblies 


11. Conclusions 

In conclusion, I would like to say that I am aware 
that only the fringes of this topic have been touched. 
However, | hope that I have been able to present at 
least a summary of the most important design problems 
and to suggest a few tentative solutions for the practical 
design of transparencies for elevated temperatures. The 
presentation was mainly from the point of view of a 
structural engineer, and I am very conscious that there 
are a great number of problems which have not been 
discussed, particularly those concerning the strength 
properties and the possible improvements in the 
existing glass and plastic materials—these are essen- 
tially problems for glass and plastics technologists. 

I would like also to emphasise that high temperature 
transparencies still present a great number of challeng- 
ing problems and in this new field there is undoubtedly 
a clear need for further research and development. 
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Some Practical Aspects ot Kinetic 
Heating Calculations 
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1. Introduction 

This paper is concerned primarily with some of the 
practical difficulties encountered in connection with the 
prediction of kinetic heating temperatures. Attention 
will be concentrated upon methods for estimating tem- 
peratures and heat transfer rates for practical aircraft 
designed to fly at Mach numbers up to about five. 

One factor common to all kinetic heating cuicula- 
tions is the variation of temperature through the 
thickness of the boundary layer, with consequent 
variation of viscosity. At Mach numbers above about 
3, these temperature variations also lead to considerable 
variations of other properties of air—which are 
commonly assumed to remain constant—even in 
classical compressible flow aerodynamics. These fac- 
tors complicate the aerodynamic equations 

Over certain regions of an aircraft (notably in the 
air intake ducts) it happens that the Jocal Mach number 
outside the boundary layer may differ considerably from 
the free stream value. Since the basic heat transfer 
relationships are derived for flat plate conditions it 
follows that allowance must be made for differences 
between local and free stream conditions. 

To estimate transient temperatures for the aircraft 
structure (and, in the limit, equilibrium temperatures) 
it is necessary to evaluate the overall effect of all the 
various modes of heat transfer which affect the struc- 
ture, as well as the heat storage capacity of the 
structure. For practical estimates it is also necessary 
to deal with any prescribed variation with time of 
boundary layer temperature and heat transfer 
coefficient. This calls for numerical methods of esti- 
mation, which must be very adaptable, yet not too 
laborious. 

To show the principles of some methods which have 
been used with success, the aerodynamic heat transfer 
relationships are first examined, with the object of 
separating the effects of skin surface temperature from 
those effects which depend on prescribed sortie condi- 
tions. The equation governing transient heating for 
shell structures is then shown. From this, a step-by- 
step method of solution is derived. Finally, some 
typical results are shown. 


NOTATION 
a_ speed of sound (ft./sec.) 
C; coefficient of skin friction 
C, specific heat of air at constant pressure 
(Fig. 1) (C.H.U./slug °K.) 
C’, temperature-mean specific heat (see text). 
(C.H.U./slug °K.) 


Hawker Aircraft Ltd.) 


E emissivity coefficient for Stefan-Boltzmann 
law 
Fy Reynolds number correction factor, equa- 
tion (14) 
h heat capacity of skin material (C.H.U./ft® 
K.) 
i enthalpy, or total heat, of air (C.H.U./slug 
K.) 
k conductivity of air (C.H.U./ft. sec. °K.) 
ky see equation (4) (equivalent to Stanton 
number) 
M Mach number=u/a 
Pr Prandtl number =(u C,/k) 
p Static pressure, absolute (Ib. / ft.) 
q_ heat flux (C.H.U./ft.*sec.) 
R_ gas constant for air (=3,089 ft. Ib./slug °K) 
R. Reynolds number 
r temperature or enthalpy recovery factor 
(equation (2)) 
temperature (° Kelvin) 
t time (seconds) 
local velocity (ft./sec.) 
y ratio of specific heats 
€ pressure recovery factor 
(Pryor) local: (pyoy) ambient 
u viscosity (slug/ft. sec.) 
p density (slug/ft.*) 
Stefan-Boltzmann constant 


(=2°78 x 10-'* C.H.U./ft.? sec. (°K)*) 


Suffixes 

central mass 

forcing 

see above 

total, or reservoir condition 

wall 

ambient 

local 

evaluated at effective enthalpy condition 
quantity N at mid-point of time interval 


x 


= 
10 


2. Development 
2.1. THERMODYNAMIC PROPERTIES OF AIR 
In normal compressible flow aerodynamics it is 
assumed that the two specific heats of air are both con- 
stant, so that the ratio (y) also is constant. It is also 
assumed that the Prandtl number (Pr) is constant. 
Figure | shows the actual variation with temperature 
of these properties, together with the variation of the 
viscosity («). This shows that for temperatures above, 
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Ficure 1. Thermal properties of air (from N.B.S. Circular 564). 


say, 600°K (corresponding with about M,=3) the 
specific heat at constant pressure (C,) rises significantly, 
while y decreases. The Prandtl number (or its square 
root, as plotted) varies only a little, except at low 
temperatures. 

It is worth noting at this point that the total heat 
(or enthalpy, i) of air (Fig. 2) is given by the integral of 
C, with respect to temperature 7. The mean value of 
C, over the temperature range zero to temperature 7, 
will be denoted by C’,, while the value of C, at this 
temperature will be denoted simply by C,,. 

Eckert’? and Monaghan®” have shown that to deal 
with the variation of specific heat it is preferable to 
express the heat transfer equations in terms of enthalpy, 
rather than temperature. The energy equation (for 
example) becomes: 


instead of 


where suffix “1” denotes conditions just outside the 
boundary layer. Also we assume an enthalpy recovery 
factor (r), as showa:— 
ip —i, 
p= 
lyor 


instead of the temperature recovery factor ry: 


Here, suffix “F” denotes the virtual “forcing” 
enthalpy or temperature. This forcing enthalpy is pri- 
marily a function of the aerodynamic flow and only to 
a minor extent is it influenced by wall conditions. It 
is felt that the forcing enthalpy is best regarded as a 
virtual enthalpy in the boundary layer, and not some 
property to be assigned to the wall. For this reason 
the name chosen for this quantity avoids the implication 
of the word “wall”, as involved in such a term as 
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ENTHAL PY, i (C.H.U,/SLUG) 
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000 4200 
TEMPERATURE 
Ficure 2. Enthalpy of air at atmospheric pressures. 


* adiabatic wall ” enthalpy, or implied by suffixes such 
as “AW” or “WO”. 

Clearly, when kinetic heating is taking place. there 
are large variations of temperature through the thick- 
ness of the boundary layer; consequently the properties 
of the air vary considerably through the layer— 
particularly the viscosity. Now Monaghan has argued 
(after Eckert) that if the properties of the fluid are 
evaluated at a certain weighted average temperature or 
enthalpy, then most of the available experimental 
results can be correlated quite well. Here, this average 
enthalpy will be called the “effective enthalpy”, and 
conditions evaluated at this enthalpy will be denoted by 
an asterisk or star. For present purposes the effective 
enthalpy is given by Monaghan’s “ mean enthalpy” 
formula : — 

i* =i, + 0°54 (iy —i,)+ 0°16 (ip —i,) 


which can be written: 


where i*,=0-16ip+0°30i,, which depends only on 
ambient and local conditions, and i*w=0°54 iy, which 
depends on the wall temperature only. It is not sugges- 
ted that the coefficients given in equation (3) are 
necessarily the best possible choice, but they are 
sufficiently plausible for present purposes. 

The heat transfer coefficient then becomes: 


k*y= 
p* u, (ig —iw) 
compared with 


kn 


Uy C, (Ty 


2.2. HEAT TRANSFER IN TERMS OF KNOWN DATA 

In the general case the enthalpies or temperatures 
outside the boundary layer are not known directly. It 
is instructive, therefore, to consider how the heat trans- 
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fer relationships may be expressed in terms of known 
quantities. 


2.3. SORTIES 

Let us consider what is known of the boundary 
layer in advance. Fig. 3 shows a typical sortie for an 
aircraft operated as a research vehicle or as an 
interceptor. 

This high-speed sortie starts with the aircraft already 
at a high altitude after a long cruise at subsonic speed, 
during which the structure has time to soak to a low 
temperature. The aircraft then accelerates rapidly and 
climbs to the specified altitude, remaining at the 
specified altitude and speed until the required distance 
has been covered. In this type of sortie the rates of 
heating and therefore the thermal stresses are at a 
maximum, since the acceleration is rapid and the initial 
structural temperatures are low. This type of sortie 
will be considered more closely later. 

Note that when the sortie has been specified, both 
the aircraft Mach number (M,) and the altitude are 
known at any time. From the altitude any of the am- 
bient conditions such as pressure (p,) or temperature 
(T,) can be found. 


2.4. LOCAL CONDITIONS 

As local velocities over certain parts of an aircraft 
may differ very much from free stream, it is necessary 
to find relationships between local and ambient 
conditions. To find local conditions over the various 
regions of the aircraft in terms of ambient conditions 
it is convenient to deal with the local Mach numbers 
(M,) outside the boundary layer. 

Figure 4 shows typical variations. It can be seen 
that for many regions M, is not very far from M,, 
although the differences are certainly significant. For 
the intake ducts, however, the local Mach number is 
very different from the aircraft Mach number. 

Now consider how the basic heat transfer equations 
may be expressed in terms of the known ambient con- 
ditions (suffix “0"), in conjunction with the local Mach 
number which is presumed known. 

The energy equation may be written: 


igor =i +4. Uy? =i, +4 
or 
— i, = 4a,” M,? 
where a,=speed of sound in free stream. 


Remembering that the temperature-average specific 
heat, C’, is defined by i=C’,T the energy equation may 
be written: 


2.5. LOCAL ENTHALPY 


Therefore i, ox 
y¥-1 G 
{1+ k (6) 
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FiGurE 3. Typical sortie (long subsonic cruise, followed by 
acceleration, climb and high-speed run). 


Note that the term in square brackets is a function 
only of the temperature or enthalpy indicated by the 
suffix. In ordinary compressible flow aerodynamics, of 
course, this term is taken to remain constant at re 4 
or for y=1-4, at a value of 0-2. Fig. 5 shows the actual 
variation of the term with local temperature. It can be 
seen that at 1,300°K (which corresponds roughly with 
total temperature at a flight Mach number of 5) the 
term has dropped about 12 per cent below the value 
usually assumed. For temperatures less than about 
600°K, corresponding with M, less than 3-0, the error 
due to taking this term equal to 0:2 is essentially less 
than 3 per cent. For larger temperatures, of course, it 
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FREE STREAM MACH NUMBER (M,) 
Ficure 4. Typical local Mach numbers in terms of free stream 
Mach number. 
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Ficure 6. Local enthalpy in terms 
of M, and M.. 
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LOCAL TEMPERATURE T,(*%) 
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The enthalpy recovery factor r* 
varies only a little, for flight 
conditions, being roughly 0°84 
for laminar flow and about 0°89 
for turbulent boundary layers. 
A rather more detailed con- 


sideration is given in_ the 
appendix. 


2.7. EFFECTIVE ENTHALPY 


It was pointed out earlier 
that part of the effective enthalpy 


for the boundary layer was in- 


dependent of the wall tempera- 


ture. This part (i*,) also can be 


expressed in similar fashion to 


FREE STREAM MACH NUMBER, Mog 


is wiser to take note of the variation, as has been done 
in Figs. 6, 7 and 8. 


2.6. FORCING ENTHALPY 

Now from equations (2), (5) and (6) it is possible to 
express the boundary layer “forcing enthalpy” (or 
“adiabatic wall enthalpy” as it is often called) in terms 
of the total enthalpy and the local Mach number as 
shown in equation (7): 


2 C, 


TEMPERATURE 
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4000 4000 000 4000 11900 “000 18,000 
ENTHALPY (CH U/SLUG) 
Ficure 5. Coefficient of M?. 


that shown for the forcing 
enthalpy, thus: 


x (0:46+0-16 [751 


=], M?). (8) 


Evidently for the special condition when the wall tem- 
perature is equal to the forcing temperature (zero heat 
transfer) the whole of the effective enthalpy can be 
expressed similarly : 


TOT 


« M,*). (9) 


It is interesting to note the similarity of equations (6) 
to (9). 

Now it is obvious that both the total and the local 
enthalpies can vary over very wide ranges, but it is 
convenient that a wide range of possible altitudes can 
be linked to a few values of ambient temperature. For 
brevity, the graphs shown are based on the standard 
stratospheric temperature. It is possible, without 
serious error, to use the enthalpy graphs for other 
atmospheric temperatures by applying the appropriate 
factor to all ordinates. 

For given ambient temperature, of course, the total 
enthalpy depends only on the free stream Mach number. 
The remaining terms of equations (6) to (9) then depend 
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Ficure 7. Boundary layer “forcing” | 
enthalpy (“adiabatic wall enthalpy”). ¥ 
| | | 
on both M, and the local Mach —* | 3 | | 
enthalpy. Any of these equa- u s* i 4 
tions, therefore, can be plotted 
graphically as a function of M, 3 8 tity 
and M, only. 3 Yo 
Figure 6 shows a graph of 3 eg ; ied 
local enthalpy in terms of M, 
and M,/M,, through equation j 
(6). The curve for M, =0 gives 8 GY _| | 
the total enthalpy. It is evident 
that local enthalpy varies over a 
wide range. Fig. 7 shows a 
similar plot of the forcing 
enthalpy, showing that i, varies | 
relatively little with local Mach + | ' 
number. 
A “carpet” plot could, | 
of course, be used to expand 0 2 3 ‘ $ 
the curves for forcing FREE STREAM MACH MQUMMER {Ma 
enthalpy. Such a presentation 
is shown in Fig. 8 for the boundary layer dependent number for a local point would have to be based on 
part of the effective enthalpy. the boundary layer thickness at the point. 
It has now been shown how many of the required “3 
boundary layer enthalpy values can be found from the Taking R.* ; 
sortie data, without reference to the wall temperature. 
The heat transfer rate equations, however, contain and using p* a wae il 
terms which depend on the wall temperature, as will a RT* 
be seen next. 
3. Heat Transfer Coefficient 9 
The heat transfer relationships will be developed ° . 
first for the more important case of turbulent boundary M,= LOCAL MACH 
layer flow. The laminar flow equations will be shown ‘ on ed 
briefly at a later stage. These developments start on 8 


the same lines as used by Monaghan. 


To= 216 66 °K 


So 
3.1. TURBULENT BOUNDARY LAYER == 
Cy* =0-288 (log,, R.*)?** ) 4 
whence ky*=0-176(log,, R.*)** (11) - 
296, 


3.2. REYNOLDS NUMBER 


1000 
\¥ 

\ 
A 
\ 


How can the effective Reynolds number be found in a — 
terms of known data? Let us assume that the Mach 
number outside the boundary layer is practically con- 
stant from the leading edge of the body to the location Ficure 8. Effective enthalpy i,* (the part independent of wall 
considered. If this were not the case, the Reynolds temperature). 
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A difficulty enters here since it is not easy to relate 
p, to M, by the treatment given for enthalpy changes. 
Naturally, in some cases the value of p, will be obvious, 
while if M, is known the pressure can always be found 
from an enthalpy-entropy chart. It is interesting, how- 
ever, to examine the relationship between ambient and 
local Reynolds numbers on an approximate basis, 
taking C, and y constant. 

To make allowance for non-isentropic compression 
from free stream conditions, the ratio of local reser- 
voir pressure to free stream reservoir pressure may be 
denoted by the local pressure recovery factor, «,. 
This factor is usually known (e.g. for air intake ducts) 
as a function of M,,. 

For constant specific heats, the effective Reynolds 
number can now be expressed thus: 


Be Mo 


T* (13) 


= (Fy R.,.) 


where 
or for y=1-4:—_. 1+ M,’\* 


Going back to the heat transfer coefficient equation 
(11), this becomes: 


ku* =0-176| log,, (Fu R...) +10g,, as) 
which can evidently be plotted as a function of the two 
arguments contained in the parentheses. The first of 
these arguments can be defined completely from the 
sortie data, but the term involving (7* «*) depends on 
the wall temperature. We now proceed to see how this 
coefficient fits into the calculation of heat transfer rates. 


3.3. HEAT TRANSFER RATE 
Equation (4) was: 


q=p* u, ky (ip —iw) 


but proceeding as above this can be re-written in 
the form: 


q= Po M, R Fy T* (ip — iw) (16) 


This equation is instructive, since it shows that the 
rate of heating is proportional to the atmospheric 
pressure, to Fy and to (ky*/7T*). Now an increase in 
Fy, also leads to a slight decrease in (ky* / T*), as shown 
in Fig. 9, but the explicit factor is much more powerful. 

In Fig. 9, the parameter quoted for the various 
curves is i*/i,. Strictly, the proper parameter is 
(T,u,)/(T*4*), but for any practical atmospheric tem- 
perature the error due to using i*/i, as parameter is 
negligible. For calculations involving various ambient 
temperatures it is more accurate to use graphs of k,*, 
instead of (k,*/T*), and make separate allowance for 
T*. It is felt that this procedure is simpler and quicker 
than that outlined in the present Royal Aeronautical 
Society data sheets. 
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Ficure 9. Heat transfer coefficient. 


Figure 10 shows some typical variations of Fy. 
The large values of Fy for the air intake are noteworthy; 
these show why intakes heat rapidly. Note that if 
M,=M.,, but ¢,<1, as for a blunt-nosed flat plate at 
high supersonic speed, then Fy is less than unity—so 
that heat transfer is reduced, and transition is delayed. 


3.4. LAMINAR BOUNDARY LAYER 
The expressions appropriate to a laminar boundary 
layer can be developed in a similar way to that given 
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for turbulent flow. The starting point is the 
familiar equation: 


ky* =4C,* (Pr*)-? (17) 
Taking 
0-332 
we get ky* = 


(Pr*y 
Using equation (13) this becomes: 


0-332 
(Fuk) 


Again, this equation can be plotted in terms of the 
expressions separated in the denominator, as shown in 
Fig. 9. 

It has now been shown how the basic equations of 
boundary layer heat transfer can be expressed in terms 
of known data, with the effects of wall temperature 
isolated, as far as necessary. It is now possible to con- 
sider how these equations may be applied to the 
calculation of transient structure temperatures. 


ky* 


4. Transient Temperatures in a Shell 
Figure 1l(a) shows a simple element of the surface 
of an aircraft structure. The basic equation involved 
is simply that of continuity, or conservation of heat. In 
literal terms this may be written for an arbitrary time 
interval : 
Rise in amount of heat stored 
heat input—heat output . , (20) 


The literal form of this equation is emphasised, for 
in a practical structure there may be a large number 
of heat input and output terms. This equation, of 
course, is exact for any time interval. The application 
of this equation will be illustrated by a simple example. 

Consider the state of knowledge about the element 
during the finite time interval from time (f) to time 
(t+ 6t), as shown in Fig. 11(b). To keep the principles 
clear, for this example it is assumed that the sole heat 
input to the element is the convective heat transfer 
from the boundary layer, while the only heat outputs 
are radiation from the outer surface and convective or 
conductive heat transfer to an interior mass. The 
diagram indicates that the boundary layer enthalpy 
(iy) is known over the whole of the time interval, but the 
temperatures of the shell (7) and of the interior (7) 
are known only at the beginning of the interval, as a 
result of the computations for the previous time interval. 

If a time interval (6f) is chosen so that the variations 
of temperature with time of the various interacting 
masses are not markedly non-linear, and such that the 
radiating surface temperature at the end of each interval 
is not much more than, say, 10 per cent greater than 
that at the beginning of the interval, the parts of the 
heat transfer equation may be written with good 
accuracy as follows. (The bar over any quantity 
denotes a time average over the interval, which is 
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assumed to be practically equal to the value of the 
quantity at a time halfway through the interval; i.e. at 
(t+ 40f).) 

Rise in heat stored =h,6Ty 
where / is the heat capacity of the shell 


(C.H.U./sq. ft. °K). 


Heat input from boundary layer 
—= p a,M Fy (ky* T*) (i; — ly 4C w) ot 


where C, is evaluated at temperature Tw. 
Heat output due to radiation 
~EoTy* (1+ =—6Tw) ot 
Ty 
Heat output by internal convection 
or conduction = K (Ty + 46Tw) ot 
(21) 
Collecting all the terms involving the shell temperature 
increment we get 
oT 
R QM. Fy (ky* / T*) (ig — iw) K (Ty —Tc) EoT 


h/at- Pad M. Fy) (ku* /T*)+4K + 2EoTy? 
(22) 


This equation looks rather horrible, but in practice it 
is easy enough to evaluate, the numerator and 
denominator being merely sums of simple terms. The 
only terms involved which are not known quite 
accurately at the beginning of the time interval are To, 
C, and (ky*/T*). To estimate these, graphs can be 
plotted as the calculations proceed, and extrapolated 
to the middle of the time interval concerned. This 
procedure amounts to taking a second approximation 
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in advance. No case has yet been encountered where 
a further iteration has been required. Extension of 
this process to any number of interacting masses and 
heat sources is obvious. It is also obvious that varia- 
tions with structure temperature of the surface 
emissivities, the structural heat capacities and the 
convective heat transfer coefficients to the interior can 
all be accommodated readily by using the appropriate 
coefficients at the mid-points of the time intervals. 

A discussion of methods for estimating natural and 
forced-convection heat transfer coefficients in the 
interior would require a separate paper. In the absence 
of a special discussion, it must suffice to state baldly 
that the estimates shown in this paper are based on 
generalisations of the usual dimensionless heat transfer 
data. 

It may be asked whether there is any overall gain in 
time or labour compared with simpler formulae which 
require shorter time intervals. Comparative calcula- 
tions have shown that for given accuracy the present 
method requires only about one-fifth the number of 
time intervals, so that the computation time overall is 
reduced appreciably. 


5. Typical Results 

Let us now look at some typical results of practical 
kinetic heating calculations. Consider one type of rear 
fuselage cross section which has to be considered (Fig. 
12). 

It can be seen that in the cross section shown there 
are a considerable number of interacting components. 
Thus surrounding the upper three-quarters of the 
periphery of the jet pipe is a radiation reflector. Out- 
side this lies the fuselage skin with immediate stringers 
and frames. At this station very little ventilating air 
passes down the fuselage between the jet-pipe and the 
radiation reflector and between the reflector and the 
fuselage skin, though most sections of this type have 
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FiGure 12. Typical cross section of rear fuselage. 


large forced convection airflows. Radiative heat trans- 
fer occurs between jet-pipe and reflector, and between 
reflector and fuselage. Convection heating of the 
ventilating air occurs in both annular passages, and the 
ventilating air warms up as it passes downstream. The 
outer skin, of course, undergoes the usual convective 
and radiative heat transfer processes on the outside. 

The lower part of the fuselage is divided into three 
compartments, inside which natural convection and 
radiation are the major heat transfer mechanisms, with 
various pipes in these compartments forming the major 
masses. The air in these lower compartments acts as 
an intermediate small mass facilitating convective heat 
transfer to the pipes. 

Results of calculations for this cross section are 
shown on Fig. 13. The temperatures attained by the 
various bits of outer skin are practically equal. It is 
notable that the pipes and the undercarriage legs in the 
lower bays barely warm up at all, although the air in 
these compartments gets quite hot. The radiation 
reflectors also get hot, while keeping the outer skin 
temperatures down. 
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FURTHER TYPICAL RESULTS 
Some more typical results are shown 


on Fig. 


14. The most interesting point about these is the 
difference of shape of the various transient tem- 
perature curves. Other significant results include 
the high intake skin temperature, due to the large 


Fy and absence of radiation cooling 


and the 


high final temperature attained by the centre 
fuselage skin. This was because the radiation 
reflector was omitted over the jet engine com- 
pressor, because this was a relatively cool part of 


the engine. 
installed as a result of this calculation. 


In practice, a reflector would be 


Figure 13. Transient temperatures for typical cross 
section of rear fuselage. 


d 

CONVECTION / \ 

a f 

8 
1.0 | AD. 

| 
| | 
| 

| | L, Z | — | 
eof 
\ | 
G 
| 
-50 


C. L. BORE KINETIC HEATING-—-PRACTICAL ASPECTS OF CALCULATIONS 645 


ase FiGure 14. Typical transient skin temperatures. 


FUSELAGE In some calculations, however, it may be desir- 
able to work as accurately as possible, so that the 
known variations should then be taken into account. 

es For laminar boundary layers, it is well known 
that the enthalpy recovery factor may be taken as 


For turbulent boundary layers, it is indicated in 

Ref. 2 that the enthalpy recovery factor tends to 

decrease with increasing Reynolds number. Taking 

— a plausible mean line through the results quoted, a 
reasonable law would appear to be: 


TEMPE RAT 


6 


r*0-898 —0-015 log, , (= ') : (24) 


LZ MINUTE > Making allowance for the fact that the quoted 
| results are based on wind tunnel experiments, and 
“$ ——_ making use of equations (12) and (13), it is possible 
to express (24) as follows: 


his sort of calculation can be repeated all over an 
aeroplane. When all the results are complete, it is 
interesting to see how all the regions fit together. It is 
then possible to see whether it is necessary in some 
localities to take account of conduction along the length This is plotted in Fig. 15, with i*/i, as parameter. 
of the skin. 


*,,* 
r*~0-98—0-015 | log,, (Fy R.,)—log,, (7 (25) 


oro 


REFERENCES 
6. Concluding Remarks 1. Eckert, R. G. (1954). Survey on Heat Transfer at High 
Speeds. W.A.D< Technical Report 54-70. April 1954. 
MONAGHAN, R. J. (1955). On the Behaviour of Boundary 
Layers at Supersonic Speeds. Proceedings of the Fifth 
Anglo-American Aeronautical Conference, Los Angeles, 


An outline has been given of means for predicting 
the transient and equilibrium temperatures in practical . 
aircraft structures subjected to kinetic heating. 


4 


This outline was considered in two parts, the first 1955. Institute of the Aeronautical Sciences, 1955. 
of which comprised a new presentation of the boundary 3. MONAGHAN. R. J. (1955). Formulae and Approximations 
layer heat transfer relationships, which includes pro- for Aerodynamic Heating Rates in High Speed Flight. 
vision for any region of the aircraft where the local Unpublished Ministry of Supply Report, Oct. 1955. 


Mach number differs from the free stream value 
Emphasis was placed on the separation of terms 
dependent on wall temperature from those which HN | 
depend only on the prescribed flight path 7 mo? 
conditions. 

The second part outlined a method for the 
estimation of transient structural temperatures. It 
was shown that this method could be applied easily 
to complicated structures with general modes of 
heat transfer between the various parts. It is simple 
to deal with variations of heat transfer coefficient, 
heat capacity and thermal emissivity. Examples 
have been given of typical results obtained by 
means of this method. 


APPENDIX 
ENTHALPY RECOVERY FACTOR 
In connection with equation (7) it was remarked 
that (for flight conditions) the enthalpy recovery 
factor varied only a little. 


FiGurE 15. Tentative enthalpy recovery factor 
(turbulent boundary layer). “CORRECTED” REYNOLDS’ NUMBER, ER 
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Advantages and Limitations of Models 


A. J. SOBEY, M.A. 


(Structures Department, Royal Aircraft Establishment) 


Summary: The use of models for structural test investigations in the presence of 
kinetic heating effects is examined. The principal features of the complex process to be 


represented are discussed under the classifications external air flow, internal heat transfer, 
elastic response. Of these the second is found to influence most model design, and an 
analysis of a typical structure is included to illustrate the various contributions to 


1. Introduction 

To establish the value of an aircraft structure 
designed to fly at supersonic speeds, certain ground 
tests will have to be undertaken. Requirements for such 
an undertaking have been given by Walker*’’, who 
has considered problems of full scale testing. 

Models are one possible source of information on 
the value of a structure and in view of the extensive use 
of them in aeronautical research in the past, it is 
pertinent to invite consideration of their assets and 
liabilities in structural testing. 

In aeronautics, models have been used more as small 
scale work pieces than as working representations of 
complete systems. As such they have enabled extensive 
experiments to be made, from among the results of 
which general observations have been made by inter- 
polation. In this way theory is verified and information 
relevant to conditions not accurately simulated may be 
determined. 

In the case of a structural test, in the presence of a 
complex system of mechanical forces under a general 
thermal field, the replica test, if it is to be of value, 
must reproduce faithfully the conditions of the 
prototype. The opportunity for sequence tests and 
interpolation is not possible. Accordingly, accurate 
representations of the thermal field are considered as 
obligatory in any satisfactory testing with models. 

The most obvious advantage that a model affords 
is the reduction in the volume of the complex equipment 
required for ground examination, should synthetic 
heating be used. Clearly, loading and heating a structure 
in a wind tunnel or in an open jet is hardly feasible 
and synthetic heating of a full scale structure is accepted. 
For the model a choice of techniques is available, and 
whether wind tunnels or artificial sources of heat are 
used, the equipment required will, of necessity, be very 
specialised. 

However, as will be shown, unless a material for 
modelling with a lower diffusivity may be found, a 
higher intensity of surface heating in the model will be 
needed, even though the actual quantity of heat required 
may be reduced. 

In order to rationalise the modelling requirements 
in the general case, some observations on models in 
general follow. 


NOTATION 
Two-letter symbols in parentheses are non- 
dimensional parametric groups, taking the first two 


internal heat transfer. 


letters of the name with which they are generally 
associated. 

The units used are the most convenient in terms 
of which to express each parameter. The number of 
basic dimensions is four, but additional dimensions are 
introduced for reasons of clarity or when justified in 
the text. 

Small case letters are used, where possible, for scale 
lengths and suffixes m, p refer to model and prototype 
respectively*. Thus /=L,,/L, defining the scale of 
length. Suffixes a, i, r, s stand for airstream, insulant, 
refrigerator and structure respectively. 

The length 1, time T, force F, mass M and heat 
flux Q quantities are total measurements. Temperature 
@ is a relative measurement so that ? implies 

@ relative to datum in model 
6 relative to datum in prototype ~ 
However, the use of the Absolute Scale for temperature 
in radiative or material property variation requires 


6 abs., model 


v= — 
abs., prototype 


V_ velocity of fluid (external or in cavity) LT~' 
a_ velocity of sound in fluid Yio 
p density L-*M 
C, specific heat MQ¢~' 
viscosity of fluid L~'T-'M 
z coefficient of linear expansion 6-? 
K thermal conductivity 
H_=heat transfer coefficient 
Stefan-Boltzmann constant L-*T-'Qoe-* 
{ surface emissivity none 
R_ interface thermal resistance L°TQ ‘6 
g gravitational acceleration sa 
Young’s modulus 
Poisson’s ratio none 
7 stress components 
U, displacements 
Strain components none 
V, body force components Lo 
thickness of insulant 


B, parts of surface (see Appendix) 


2. Review of Models 

Two systems are similar if, in the general case, scale 
lengths exist (as many as the order of the base of 
dimensions describing the physical process) in terms of 


*Except C,. 


we 4 
! 
7 


KINETIC 


A. J. SOBEY 


which all physical items in the one system have values 
determined (i.e. derived) in terms of the basic scales. 
This ideal similarity constitutes a projectivity between 
kindred systems and is essentially geometric in inspira- 
tion. This viewpoint is admirably presented in Duncan’s 
book *’ using the vector concept of base dimension and 
measure formula. Such a system can rarely be achieved 
except in simple models with a limited field of enquiry. 

In fact a physical model in which work units may 
have different scales depending on their origin is not 
unusual. In the thermal structure the energy stored 
elastically is of an entirely different order of magnitude 
to the thermal energy and may be considered 
independent of it, so that thermal effects based on one 
set of scale lengths may coexist with elastic effects based 
on a different system, e.g. as far as the thermal structure 
is concerned, the base of dimensions is four for purely 
thermal effects, plus a fifth appropriate to the elastic 
case, say the force unit. If, further, dynamical effects 
are present, of which the mechanical forces are part, 
the system of scales adopted must be consistent with 
the elastic similarity established. The mixing of units 
is permissible if, and only if, there are grounds for 
dissociating phenomena, and that each class of pheno- 
mena is dimensionally self-consistent. A second instance 
of this occurs in the process of natural convection where 
the mass unit associated with the fluid motion is, in 
practice, independent of the unit of heat transported 
(see e.g. Bosworth’). 
3. Limitations on Experimental Design 

Whenever it is necessary to regard as significant a 
physical parameter, a condition is imposed on the 
choice of scale lengths which is in the nature of an 
experimental constraint. Each such constraint reduces 
the number of such degrees of freedom in the choice of 
scale lengths, and clearly the number of degrees of 
arbitrariness in choice of scale length is the excess of 
base dimension over number of constraints; if this is 
zero a unique model is possible. In particular we note 
that if this excess is negative (e.g. number WN of significant 
parameters (including base dimensions) is greater than 
twice the order of the base dimension M) then the 
existence of a model is fortuitous. 

To facilitate model design we require to keep down 
N. When N=2M and a unique model is possible, it 
is important that at least one parameter shall change, 
otherwise the model is identified with the prototype. 
In particular, if material properties are important, then 
the number of such parameters which may be scaled 
cannot be as great as M, unless the material is changed. 

The importance of material parameters in model 
design is obvious. The choice of material implies the 
acceptance of certain parameters and the appropriate 
scaling conditions which go with them. If the number 
of these is large, difficulties are apparent, and in the 
thermal transfer/elastic response similitude analysis, 
there are numerous such internal parameters which 
dominate similitude (e.g. K, C,, », 

Although five such material properties are significant 
in this problem, they do not determine scale lengths 
entirely since elastic work and thermal work are treated 
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as independent. Further z and E (which, like the non- 
dimensional v, are of no interest in the heat transfer 
problem) only occur in disjunction when strain must be 
scaled as one to one, which is not obligatory in many 
of the practical cases. 

Accordingly we can expect two constraints from 
material properties on the heat transfer process (leaving 
two degrees of freedom in the choice of scales) and a 
supplementary condition based on elastic similarity 
which is taken up in an essentially independent fifth 
scale (say elastic work, or force). Since E, 2, v are all 
significant parameters peculiar to the elastic problem, 
there may be additional constraints on the thermal 
similarity process induced by elastic requirements. For 
example, if v varies with temperature in a random 
manner, and » is important, then a one-to-one tempera- 
ture scale would have to be accepted. (For practical 
materials, such as light alloy or steel, there is no 
evidence of any significant variation of v over the 
working range.) 


4. Thermo-Elastic Similitude 

The similitude conditions appropriate to the 
statically /dynamically loaded structure in the presence 
of thermal effects have been considered by several writers 
including O’Sullivan*’’, Madejski“*’ and Heldenfels"”’. 
The first two writers consider heat transfer within the 
structure by conduction alone, but Heldenfels notes the 
significance of external radiation and joints. 

Many writers have given similitude criteria for 
general heat transfer, including Bosworth®’, Fishenden 
and Saunders*’ and usually from the 
viewpoint of understanding the convective transfer 
process as a substitute for precise quantitative analysis, 
and other writers have considered dimensional analysis 
as applied to elastic phenomena. These include 
Goodier’*, Durelli, Phillips and Tsao’ and by the 
three authors cited who consider the compounding of 
similitude requirements for thermal and elastic effects. 
In addition thermoelastic similarity has been discussed 
by Mindlin and Salvadori"*’ and by Green, Radok and 
Rivlin’”’. 

In the present paper, emphasis on the internal heat 
transfer process is made, so that similarity conditions 
for thermoelastic phenomena are dominated by accept- 
ance of the need for accuracy in heat transference. When 
this is done, the choice of scale lengths is so restricted 
(except in the simplest case) that the elastic process is 
seen as supplementary and scaled using as many of the 
scales derived from thermal arguments as is possible. 
5. Similarity Criteria—Supersonic 

Aerodynamics 

The nature of supersonic flow is well understood 
and the prediction of heat transfer from the supersonic 
boundary layer has been given extensive mathematical 
treatment and experimental verification. (See e.g. 


Taylor**’ and the references given there to the work 
of Monaghan and Eckert.) 

The general survey (on a dimensional basis) shows 
that the non-dimensional heat transfer coefficient (the 


: 
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Stanton number, S*) is a function of the Reynolds 
number, Re, Prandtl number, Pr, and the Mach number, 
Ma, for a given profile 


where St=H/pVC, 
Re=VpL/u 
Pr=pC,/K 
Ma=V/a 


and all the dimensioned parameters refer to the air. 
In a practical attempt to simulate the supersonic flow, 
it will, in general, be necessary to have similar Re, Pr, 
Ma in order to obtain a correct St. If this is attempted 
in a wind tunnel using air as the working fluid, it is 
usually impossible to arrange for Re and Ma both to 
be similar. Moreover, the choice of units appropriate 
to the air and to the structure must be the same, e.g. 
time and temperature scales must be consistent in fluid 
and structure [i,=%,]. 

Because of the well established validity of the 
supersonic flow theory and its mathematical interpre- 
tation, the computed heat transfer may be a more 
satisfactory basis for experimental heating tests than 
attempting general aerodynamic similitude. Unless 
properties of the airstream are significant (as in aero- 
elastic tests), the synthesis of heat transfer can be 
considered, using, for example, lamps. (See Taylor**’.) 


6. Similarity Criteria—Internal Heat 
Transfer 

Thermal conduction, being dependent on material 
properties (thermal conductivity and thermal capacity) 
dominates the modelling in as much as conductive heat 
transfer is always important and the scaling of two 
parameters limits the choice of modelling scales, two 
constraints being imposed. Further, surface effects may 
be important, such as radiation and interface resistance, 
together with transfer media in contact with the surface 
(Fig. 1). 

Accordingly the modelling criteria are developed 
from the central idea of conduction always being present. 
Radiation is added both as an internal distributive 
phenomenon and as an influence on external heat input. 
Joints are treated as an adjunct to conduction. 

If the transfer to and from, and within the conductive 
parts of the structure is simulated, the choice of scale 
lengths is seen to be very limited, if possible, and 
convection, insulation and refrigerative elements are 
simulated as additional phenomena to the main transfer 
process. 

When heat transfer within the structure is by con- 
duction alone, we require the identity in the two systems 
of the 


KT 

pC,L? 

and the Biot number = Bi=HL/K. 
The former requirement governs the correctness of 
the thermal diffusion within the structure and the latter 


ensures the appropriateness of the heat intensity on the 
surface of structure. These conditions follow at once 


Fourier number = Fo = 
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Ficure 1. Schematic representation of heated aircraft structure. 


by noting if K, and (pC,), are to scale consistently 
so that the scale for time, f, is chosen to be 
P (K/ pCy)n/ / pCy), 
and the heat flux, temperature scales are related by 
G=PI 

These conditions are consistent with a surface 

intensity of heating which will scale as K ['? which 

p 
in most practical cases will exceed unity. In order to 
avoid this it is clearly necessary to use a material of 
lower conductivity in the model. To avoid compression 
of the time scale, a lower diffusivity is required. Fig. 2 
indicates the order of diffusivity of some common 
materials and it will be seen that avoiding too severe a 
compression of the time scale requires the use of 
materials such as steel for light alloy, with some non- 
metallic medium as a working material for steel, e.g. if 
the prototype is light alloy and the model is in stainless 
steel, pure conductive requirements are satisfied with 
scales g=1°52/, and the heat transfer 
coefficient scales as 0-0867/~-'. Surface intensity of 
heating scales is 0-0867/- "i. 

This system would be feasible for thermal (analogue) 
investigations were it not that the scaling requirements 
may not be sustained when allowance is made for the 
material property changes with temperature. 

If radiation and conduction are both significant, as 
in an integral construction, evacuated box, we require 
identity of four criteria: namely, two conductive and 
two radiative. 

HL/K, KT/L?, C,, & and o@L/K must be similar. 
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If compatible scale lengths are found, /, gq. 
require 


(K,,/K,) 


from which we deduc 


Kn 
l 


e.g. L.A. prototype, stainless steel model yields 


[=0-09)-*, t=0-141) 


CONCLUSIONS 

If the original material is used, and we are forced 
to use an identical temperature scale, the only model is 
the identity. For the scales given the surface intensity 
of heating will scale as #*. Thus for a reduced intensity 
we require <1 and a reduced total heat input is best 
achieved by using a low conductivity model, e.g. stain- 
less steel for L.A. or even glass for stainless steel. 

If the radiation is a direct gain or loss from an 
external surface* rather than an exchange between 
interior surfaces, we may relax the condition of 
independent satisfaction of 


© and 


to read (6 c@L/K), 


This may be assumed whenever the “‘form factor’ for 
the radiation transfer is not critical upon the greyness 
of the surface to give 
pax Km 9-3. q=(*") {pC 

Kytom K, 

(K,C,)» 

The form factor can be estimated (see e.g. Fishenden 
and Saunders”®’, Bosworth®’, McAdams“”) and _ its 
dependence on greyness determined. 

When the material properties vary with temperature, 
similarity of two systems is ensured if the two systems 
are described by a common non-dimensional boundary 
value problem. The heat conduction uniqueness 
theorem for this case is given in the Appendix. 


In the general conduction case, similarity of tw 
systems (m, p) is assured if 


P=t and d=] 


*When only the radiation from the external surface is 
significant, and synthetic methods of heating the structure are 
used, the radiative constraint on scaling can be avoided, and 
design proceeds as for the purely conductive model since the 
appropriate surface heating technique may be presumed to be 
adjusted to allow for the radiative gain or loss. 
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The boundary conditions impose: 

over B, noconstraint 

over B, a constraint on the impressed heat intensity 

over B, a constraint on the heat transfer coefficient. 
These are 

lti=1 
and 
H,./H,=i-". 


We note that the requirement that 
K.=K,, (pC,).=GC,), 


leads to scale lengths with only one degree of freedom, 
namely, 
(qg,/,t, are 1, FP, 1). 


If one further condition is used to determine / then this 
condition must be between physical properties different 
in the two systems else g=/=t=i=1 and no modelling 
is possible. This eliminates the possibility of further 
scaling of 

(a) absolute items (e.g. radiation), 

(b) further material properties (e.g. 2). 


The conclusion reached is that if the variation of 
material properties is important, the imposition of 
identity in the temperature scale admits of similitude 
only in the conductive case. This is of interest only 
for high solidity structures such as thin solid wings 
where radiation effects are absent. 

Is #=1 obligatory? To avoid this restriction we 
must examine the temperature dependence of significant 
material properties on the absolute temperature. 

The characteristic temperature dependence of 
material properties on temperature is given in Fig. 3* 
from which we note that over quite large and practical 
ranges of temperature, the four properties shown are 
significantly power laws in the absolute temperature. 
This can be used to scale temperature, namely: 


if K is K,@ 


(pC,) is (pC,),O* 
and 
16, (6,,, 8, on Absolute Scale) 
then 
K.=K,. (pC,). = 


and the heat conduction equation in the model is 


6 


Tn 


> 
Cc 
divin (Kn grad, Om) 


00 

which is identical with the equation for the prototype if 


On the boundary we require that the heat inflow shall 
be suitably scaled, i.e. HL/K the same in both systems, 


i.e. H,=1-'H, 


*Prepared from Refs. 16 and 19. 


; 
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Ficure 3. Variation of material properties with absolute 
temperature. 


from which there follow for the basic scale lengths: 
q="'P, Land # arbitrary (two degrees of 
freedom for pure conductive similarity). 

If, in addition, we have radiation effects of import- 
ance so that § has to be the same and also the radiant 
conductive heat units are to scale, then 


In practice the indices A, « are quite small, and the 
problem of complying with the requirements for this 
case is not significantly changed from the isotropic case. 

Practical values for A, « are given in Table I. 

For modelling with a reduced linear scale, we require 
ie. (i.e. surface heat intensity scale > 1). 
If the prototype is operating at — 56°C (i.e. 6,=217°K) 
and the model in the laboratory at 20°C (i.e. @,, =293°K) 
#=1-35 and / is 0-43 for L.A. and 0-46 for steel. These 
scale lengths are unlikely to be useful in practice for the 
aircraft case. Further the time scale is severely shortened 
and likelihood of useful data being obtained diminished. 

For a significant modelling scale, /, we need a large 
change in the emissivity, and the blackening of interior 


TaBLe I 
ae BA. Stainless Steel 
(D.T.D.683) (D.T.D.176) 
"Temperature range 0-200°C 0-500°C 
024 
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faces of a model may lead to a practicable scaling. 
However, the prototype surface conditions may be 
temperature-dependent (McAdams) and the inconstancy 
of { may prove a major limitation. Rough steel surfaces 
have {. of the order of 0-5 so that a further 50 per cent 
reduction in / may be possible. The experimental 
prospects are, however, slim. 


7. Allowance for Resistive Joints 


In all the cases considered (conduction with or 
without radiation, isotropic or non-isotropic materials) 
the effect of joints has been ignored. They are readily 
incorporated into the scaled system when the scale 
lengths are determinate, namely, 


R,,/R,=1 x 

This relation must be independent of temperature and 
will be so, for the cases discussed, if either R is constant, 
or a power law in the Absolute temperature. 

Thus a smaller scale model demands a reduction in 
joint resistance for the same material, which in practice 
will prove very difficult to accommodate. If joints are 
important, they are almost certain to prove obstructive 
in model design. 


8. Influence of Convection 

The general (Dimensional Analysis) criteria require 
the identity in model and prototype of the Nusselt 
number, Nu, the Grashof number, Gr, and the Prandtl 
number, Pr, given by 


Nu=(HL/K) | where the physical proper- 

a ae ties K, «, and so on, refer 
| to circulating (convecting) 
Pr =:(uC,/K) 


fluid in cavity 

If conductive (and/or radiative) effects have been used 
to fix the scales for g, /, t and # in the structure. these 
parameters fix properties of the fluid. 

The Nu criterion ensures that the units of heat flux 
in solid and fluid are consistent. 

Gr and Pr ensure similarity of the flow process. 
It is clear that to satisfy all these criteria a fluid with an 
unusual range of parameters will be required. Attempts 
to simplify these requirements by insisting only that Nu 
and GrPr be identical as suggested by Bosworth” 
will alleviate somewhat the experimental difficulty. 

The physical properties of fluids are not variable 
over the range likely to be demanded and if natural 
convection is known to be important, a mismatching of 
the convective transfer is likely. 


9. Insulants 

In practice, an insulating layer is a region of very 
low conductivity thinly applied so that surface-wise heat 
transfer is negligible compared with normal heat transfer. 
Accordingly the full simulation of thermal diffusion in 
the insulant may be ignored and the thickness of the 
insulant adjusted so that Ht,/K, remains constant in the 
two systems. This should present no practical difficulty 
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unless a poor conductor is used as a model material 
(e.g. glass for steel). 


10. Heat Sinks 

These may be of two kinds: - 

(a) Heat absorption associated with change of state 
of the sink material. This requires the tem- 
perature of the change of state to scale (Abs.) 
as 3. This should not present practical diffi- 
culties. If the heat sink completely changes 
state, complications arise, but this case is not 
of interest. 

(5) Heat is absorbed by a refrigerative element 
according to 


QO, =H, (6, —8) 


This requires H,L,/K, to scale in order to have 
a consistency of the thermal units. Implicitly 
=O 


rm 


11. Similarity Criteria—Elastic Effects 

In the general isotropic case we have P, E, v, L, 
7, 2, 9 involved, requiring 7,,/E, P/EL*, 26 to be 
consistent in both model and prototype. 

In this context F, v and 2 are assumed independent 
of 6. Thus 


E f(p/EL’, a6). 


Under special conditions, we may drop the depend- 
ence of 7,, on v and may expect to change the scale for 
(z#) or strain. To do this we require a mathematical 
examination of the boundary value problems associated 
with (linear) small-strain elasticity. Similarity criteria 
based on this kind of assumption have been given by 
Mindlin and Salvadori“*’, and by Green, Radok and 
Rivlin“”’, but the non-isotropy of the material properties 
has not been investigated. When material properties 
vary with temperature, these requirements imply that if 
the original material is used, = 1, which has been seen 
in the section on heat transfer to be necessary when 
radiation and convection are important. To examine 
#1 requires the detailed elaboration of the linear 
problem of elasticity with its less restrictive criteria. To 
justify the similarity criteria based on a common form 
we require a uniqueness theorem, which is given in the 
Appendix. 

If two elastic systems have the same non-dimensional 
boundary value equations they lead to the same stress 
system. 

Elastic Similarity Conditions in a three-dimensional 
medium with elastic constants dependent on temperature 
can be most readily developed by the following 
argument. 


using the tensor notation of Sokolinkoff@*. 

The compatibility of strain and equilibrium of stress 
subject to certain stress boundary conditions ensure the 
correctness of the stress/strain system. 
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Compatibility requires 
mn + Emn, =Eim, mgt Enjim 
Equilibrium 
ti, 4+ (3) 


where V; are the body force vector components per 
unit volume. 

If we arrange that 

at all points, e being a constant, then for isotropic length 
scaling, if (2) is satisfied in the model it is automatically 
satisfied in the prototype. 

For equation (4) to hold universally 


m 
) 
kh (6 
(24) m ( x0), (7) 
From (5), (6), 
(Tim ale ’ En (E26) 
Scale for stress = « E, ~ (E20), (Eo), (9) 


If F and z are together proportional to some power of 
the absolute temperature, (9) gives a constant scaling 
factor for stress even when temperature is scaled. 


From equation (3) 


IV —1 


giving the scale for body force intensities. 

If we are dealing with two-dimensional stress systems 
in the isotropic case, it is well known that for the 
problem with boundary conditions or stress, v is unim- 
portant for singly connected regions since the Airy 
stress function (stress potential) is independent of »v. 
Only when conditions at the boundaries involve displace- 
ment is v important, or when we have multiply connected 
regions when the displacement analysis is implicit. 

Thus the requirement that v is to be reproduced is 
important in all cases except the special one of plane 
strain or generalised plane stress (when v may be 
arbitrary) for singly connected zones. This has more 
pertinence to components than full scale structures. 
The special conditions on the scales for plane strain 
and generalised plane stress have been given by Green, 
Radok and Rivlin“® who use an advanced treatment 
based on complex stress potentials. Their conclusions 
—that » is important except in the single connected 
case—are that no subset of the similarity criteria holds. 
This follows from connectivity arguments of the kind 
given. 

The scale for stress implies a scale for force, and 
also for work. Units of energy can be mixed in this 
way (heat and work having, in general, different scales) 
provided that there is no significant interaction. This 
has been considered by Hemp“*’, and by Broglio and 
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Santini*’ and found not to be significant for the non- 
thermal shock case. 

If we have dynamical effects, e.g. in the analysis of 
the natural frequencies of systems, we should amend the 
equations of equilibrium to read: 


leading to 


(pU) m/(pU), = (E2), 
But the scale for strain 
whence 
Un = id 
U, 7 


gives scale for displacements and therefore 
Pm | Pp = Emt? / El? 

gives scale for density. This is a condition which will 
be difficult to achieve in practice owing to impositions 
on the time scale from conduction considerations. 

Elastic properties may be examined for the elasto- 
static case, but dynamical effects will be difficult to 
accommodate. 


12. Plasticity 

When the elastic limit is exceeded and time- 
dependent strains are significant, the cumulative strain 
is a function of time, stress, temperature, and, possibly, 
the strain already produced. Madejski®*’ has examined 
an idealised case, but similarity will, in general, require 
an identity of the strain, stress, temperature scales, and 
possibly of time as well (if the strain rate is dependent 
on the cumulative strain as in work hardening). Such 
conditions are not tenable except in idealised com- 
ponent study and the presence of plastic strains is 
regarded as a major limitation on the design of 
experiments. 


13. A Design Study 

To show the likely significance of some of the 
internal heat transfer parameters discussed, calculations 
have been made using the structure of the original 
Hoff-Torda*’ study modified in the following ways: 


(i) By the inclusion of web-flange interface 
resistances. 


(ii) By allowance for radiation heat transfer 
between flange and web. 


(iii) By allowance for free convective heat transfer 
between flange and web. 


The geometry of the system is given in Fig. 4. The 
material is steel and the significant properties are as 
follows. The structure is subjected to an instantaneous 
heating according to 


Q =H Onange) 
Q being the heat intensity 
H the heat transfer coefficient 
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Ficure 4. The Hoff study wing. 


4, is the zero heat transfer wall temperature 
(taken as 300°C). 


Following Hoff, H is taken as 90 CHU/ft.* hr. °C, a 
heat transfer coefficient consistent with Ma=5, Re=10' 
at 40,000 ft. (see Taylor @*’). 

This example has been discussed extensively in the 
literature, in addition to the original reference, by 
Parkes®'’ who added the effect of flange-web solidity 
variation, by Pohle and Oliver®* who gave the exact 
conductive solution, and by Biot’ who used the instal- 
lation as an application of his variational technique. 

In none of the papers cited is the effect of joint 
resistance, internal convection or radiation examined. 

The calculations include the effect of a thermal 
resistance joint at the web-flange junction, the resistance 
of which is varied over a wide range (of 200 to 1) and 
is likely to include all practical constructions. See, e.g. 
Barzelay and Holloway’ and Barzelay, Tong and 
Hollo®’, who give measured values of joint resistances 
for sheet constructions of the stringer-panel type. 
Bernard®) gives an indication of the effect of impact 
pressure on the joint resistance for blocks. 

The radiative transfer is calculated using the classical 
theory (see, e.g. Fishenden and Saunders” or 
Bosworth®’, and as Hoff’) assuming an internal surface 
emissivity of one half. The evidence available on the 
likely value of { is scant, but McAdams"” gives an 
extensive list of known values which bracket the one-half 
region. As this survey is intended only to be com- 
parative, the value chosen is not critical, and the selected 
value is thought to be representative of typical interior 
surfaces. The problem of the dispersal of non-absorbed 
radiation is ignored, and is unlikely to have a significant 
influence on the trends of the calculations. 

To estimate convection effects, the following hypo- 
theses have been used: 


(i) The thermal capacity of the convective fluid 
is Zero. 


(ii) The wall fluid heat transfer coefficient may be 
estimated from flat plate formulae, e.g. a value 
of heat transfer coefficient of H = 1-2 CHU / ft.’ 
hr. °C (McAdams"” and Zemansky*’) so that 
the heat transfer between walls is at the rate of 
(—meant+ Where iS the effective 
transfer temperature of the convective fluid 
and is calculated on the basis of assumption (i). 


Numerical analyses have been conducted using a 
matrix technique suitable for use on an electronic com- 
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puter. Results are presented for a range of resistances 
from 10°* to 10 C hr. per CHU, which 
covers the measured range of values (Figs. 5, 6). A 
likely value of interface resistance for riveted or bolted 
constructions is 5x 10 C ft.* hr. per CHU’, at 
which value a significant transfer is brought about by 
convection and radiation. ; 

This particular example has been selected because of 
its historical significance and the relative solidity of the 
flange and web imply (for the heating rate given) a 
thermal stress which is very high (in the Hoff-Torda case 
of about 80 per cent the maximum possible zF times the 
excess of adiabatic wall temperature over initial tem- 
perature). In this case the effect of the joint is to 
increase the maximum tensile stress from 80 per cent 
to 90 per cent of the maximum possible over the extreme 
range of joint resistance. The resulting alleviation of 
stress by flange-web transfer (convection and radiation) 
is shown in Fig. 7. 

When the flange soak tem- 
perature is raised, the effect of 
radiation and convection is 
intensified and a survey of the 
relative proportions of heat 
received by the web due to pure 
conduction across the joint, to 
radiation, and to convection is 
given in Fig. 8 as a function of 
flange temperature and joint 


THE SHADED 


resistance. TURES THE 


Ficure §. Temperature distributions 
through the study wing as a function 
of joint resistance. 


The evidence of this example 
is that for the jointed structure, 
the simplification of heat transfer 
to a purely conductive process is 
misleading and a proper allow- 
ance for all these effects will 
have to be made in any practical 


testing. 


14. Conclusions 
In the simplest of cases, a 
solid structure, the internal heat 
transfer may be represented 
satisfactorily by a scale model 
with two degrees of freedom. 
Allowance for surface radiation 
(when aerodynamic sources of 
heat are used) reduces the model 
to a one degree of freedom one. Size of model may 
be chosen arbitrarily. In this case the elastic properties 
may be represented also. Change of material is possible 
only when the material properties are constant through- 
out the experiment, i.e. when they are independent of 
temperature. 

When properties vary with temperature it is possible 
to use the same material with an identical temperature 
scale. In this case the conductive representation leaves 
only one degree of freedom and external radiation 
cannot be catered for. For practical materials a shift 
in temperature scale is possible owing to the special 
nature of the temperature dependence of the appropriate 
physical properties. When this is done a model is 
possible, but when conduction and radiation are allowed 
for, a unique scaling is determined. This leads to a 
smaller model only if the absolute temperature is 
increased in scale. This will be the case for stratosphere 
to laboratory scaling. Elastic effects are scaled as before. 
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Figure 6. Effect of radiation and 
convection on the heat transfer 
within the study wing. 


7) ant 


i 
|| 
1 
653 
< 
/ 
stout 
he 
te 
: | 
| 
4 
Tina (anes ) 
ag 
leith 
—aapation 
= 
CONVECTION 
| 
A>” 
= 
° 
BS 


JOURNAL OF THE ROYAL AERONAUTICAL SOCIETY 


CONDUCTION ONLY 


ZERO JOINT RESISTANCE CASE 


CONDUCTION 
PLUS 
RADIATION (t %) 


CONDUCTION 
PLUS 

DIATION 
PLUS 

CONVECTION 


2110°* 10°? 
fog 


Ficure 7. Estimates of maximum thermal stress (tensile) as a 
function of joint resistance. 


If a high solidity structure is built up of components 
whose faces of contact are not on surfaces of zero heat 
transfer then the interface resistance will have to be 
scaled, and this will be very difficult in practice. For a 
reduced linear scale, the joint resistance has to be 
reduced unless a change of material is possible to one 
of greater conductivity. 

For structures of low solidity the conditions of the 
examples given are complicated by internal radiation 
and convection, which are of increasing significance as 
joint resistance increases. If these phenomena are to be 
scaled, the experimental constraints are severe. A unique 
scaling is possible when the internal radiation is correctly 
simulated. To arrange for the convective process to be 
simulated correctly requires the fluid properties to 
change, and the appropriate fluid may not be available. 

Elastic effects are correctly reproduced in the small 
strain case with a change of temperature scale owing 
to the apparent constancy of Poisson's ratio and the 
power law dependence of F with temperature over the 
working temperature range. Effects of interest outside 
the small strain theory (e.g. post buckling behaviour 
and permanent set) will require a one-to-one temperature 
scaling and restrict the field of application to purely 
conductive heat transfer cases, e.g. the solid structure. 


APPENDIX I 
UNIQUENESS THEOREM FOR THERMAL CONDUCTION 
In the non-isotropic case we may no longer use the 
method of dimensional analysis but boundary value 
problem analysis is required. 
The heat conduction equation in the case of non- 
isotropic and/or temperature-dependent properties follows 
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FicureE 8. Comparative study of heat transfer from flange to 
web for various joint resistances. (The ordinate in each graph 
indicates the proportion of the total heat received by the web 
due to the three agencies indicated.) 


at once from the conservation of energy within a volume 
V, as 


06 
div (K grad @) - - ; (1) 
a 


and the appropriate boundary conditions of practical 
interest are: 


over part B, of the boundary, (x, y, z, T) specified 
06 


over part B, of the boundary, 


=@ (x, y, z. T) specified | 2) 


over part B, of the boundary,.- =H’ (6,—6) specified 
On 
(n being a unit normal vector). 


Subject to (2), (1) has a unique solution. For if there are 
two possible solutions 


6, (x, y, z, T) and 6, (x, y, z, T) 
with common 
6, (x, y, z,0) in V 


satisfying (1) and (2), then the difference 


satisfies 


) 


div (K grad 6) = pC 


06 
on =0 on = on B 


an 
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Consider 
Ke as | grad as— | ds 
on 


B,+B,+B 


which is negative semi-definite. 
By the divergence theorem this is 


| (grad 6)°dV 4 ( ab 


Initially | pC, dV =O since 6, (x, y, z, 0) = @, (x, y, z, 9) 


he derivative of {{§{ pC,@*dV cannot therefore be negative 
with respect to time. Hence 4 and grad@ vanish every- 
where 

This uniqueness theorem is independent of the non- 
isotropy of K and pC, with respect to space and tempera- 
ture. Obviously if K is time-dependent the proof fails, 
but this is of no interest. Thus we conclude that if two 
systems satisfy equations like (1) and boundary conditions 
(2) which can be reduced to a common non-dimensional 
form, then they are physically similar irrespective of the 
non-isotropy. 


If T 
then (1) reduces to 


div (K grad @) 


co 


APPENDIX II 


UNIQUENESS THEOREM FOR ELASTIC MATERIAL UNDER 
THERMAL STRESS FIELD (SMALL STRAIN THEORY) 


Equilibrium 
Strain 


Stress strain 


Subject to boundary conditions 


U,=U,, over part B, of the boundary 


‘ (8) 
T,=ryn,;=T,, over part B, of the boundary { 


Then if two solutions (U,),, (U,), of this system exist with 
associated stresses (7;,),, the difference system 


U,=(U)), —(U)). 
(9) 
T= (74), 
with strains ¢, =(¢,,), —(¢y), 
satisfies 
I E 


Uniqueness then follows in the usual way by considering 


LB, 


2 | | 
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since the strain energy is positive definite. 
Over B, U,=0, over B, 7,n;=9; 


hence {\\ U;7,,n,dS=0 and the two stress systems are 
identical. 

This proof shows in addition that if (U,),~(U)),, the 
only difference possible is a displacement vector leading to 
zero stress, i.e. a rigid body displacement. 

Thus a common non-dimensional boundary value 
problem structure ensures similarity of strain and stress, 
and of displacement if the systems are suitably orientated. 

Further the temperature dependence or non-constancy 
of E, » does not affect the proof 
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Summary of the Discussion 


MAUREEN E. MICHAEL, M.A., Grad.R.Ae.S. and JOAN BRUCE, B.Sc., A.F.R.Ae-.S. 


(Technical Department, Royal Aeronautical Society) 


SOME FACTORS AFFECTING THE CHOICE OF MATERIALS 

The first speaker in the discussion following Mr. Heath’s 
paper was from Bristol Aircraft Ltd.; he noted that the 
paper could be divided into three parts, material properties, 
time-dependent effects and materials of the future. Com- 
menting on the first part, he agreed that it would be a good 
idea to throw away existing curves of specific strength, but 
he would like someone to suggest a convenient method of 
presentation for the data on material properties necessary 
to the designer. In the second part of his paper Mr. Heath 
seemed to suggest that structures should be designed on 
the “fail-safe” principle because of the difficulty in 
achieving a full test life with some specimens. Surely the 
need for “fail-safe” structures arose primarily from other 
considerations? In giving guidance on material selection 
Mr. Heath had concentrated on the properties af tempera- 
ture, but it was possible that in future aircraft the plan 
form might be such that at high speeds when heating was 
present the gust loads would be small, and the design 
criteria might well be landing and take-off loads when the 
material would be cool. He agreed with the lecturer that 
the present lack of knowledge on creep was rather fright- 
ening, and put in a plea for more compression creep 
testing, as in some cases compression buckling could be 
important. It could be shown that before creep buckling 
collapse there was a significant amount of lateral deforma- 
tion and this extended over a considerable period of time 
before collapse, so that it would be possible to use 
inspection methods. 

In his reply Mr. Heath said that it would be difficult 
to present data on material properties in a readily usable 
form; curves showing stiffness and tension and compression 
proof stresses would have to be given for each material for 
various lengths of time at any particular temperature. 
With some of the new materials it was remarkably difficult 
to ascertain their properties due to lack of published data, 
but they would probably have to be used for structures 
subject to temperatures above 500°C. Later on in the 
discussion he stressed the point that choice of material 
and choice of type of structure must go together, as one 
could not get a fair picture by looking solely at material 
parameters. The materials mentioned in his paper were 
those for which he had information readily available and 
they were to be regarded only as an iliustrative selection. 


He feared that where some of the new materials were con- 
cerned it might ultimately be the workshops that governed 
the choice. 

He agreed that the primary consideration in “fail-safe” 
design was the same as at room temperature; hig reason 
was merely an extra one for such design. 

Mr. Heath endorsed the plea for more compression 
creep testing and, in reply to another speaker, from A. \ 
Roe and Co. Ltd., suggested that the right proportion 
might be 60 per cent compression testing and 40 per cent 
tension testing. 

Asked to give some indication as to which were the 
most important material properties for the designer to 
consider, he said he could not really help; he felt that the 
answer to such a question could only come from a large 
scale research programme. Although there were many 
establishments in this country carrying out research there 
was little or no co-ordination of their activities. Only by 
setting up a co-ordinating body could they hope to get 
such programmes as the one he envisaged carried out. He 
made a strong plea for such a body to be set up. 

In replying to a question from the Chairman, Mr. 
Heath agreed that in testing large built-up structures sub- 
ject to very severe kinetic heating, it would be impossible 
to separate fatigue and creep testing. 

A speaker from the R.A.E. commented that, unlike the 
engine people who thought they knew all about creep, 
structural engineers did not want to know a great deal about 
the fundamentals of it. They wanted to know what creep 
problems would apply to structural elements used in 
supersonic transports. 

A spokesman from Bristols considered the case of a 
supersonic transport having a pressurised aluminium alloy 
cabin. Taking into account a scatter factor he deduced a 
safe working stress due to pressure. In supersonic aircraft 
with thermal stresses, the working stress would have to be 
even lower. From a structural point of view the problem 
was simplified. The problem became not how to deal with 
these stresses when they arose, but how to stop them 
arising. He suggested that this could be done by cooling or 
by insulation to smooth out the thermal gradients. 

A questioner from the Air Registration Board asked 
how the rapidly changing qualities of thermal structures 


4 
: 
21. 25. 
22. 26. 
3 
27. 
24. 28. 


KINETIC HEATING-—-SUMMARY OF DISCUSSION 


would affect the repair schemes for future supersonic 
transports. Mr. Heath replied that at the present time the 
only answer he could give was to put on a patch and 
hope for the best 

In reply to a questioner from Bristol Aircraft on the 
need for consideration of similitude when dealing with the 
specific strength curves, Mr. Heath agreed, and repeated 
again the point he had made in his paper, that the con- 
sideration of one facet was not enough. A completely 
misleading answer could be obtained if only one factor 
were taken into consideration 

In correspondence after the meeting it was suggested 
that basic fatigue data should be obtained on materials 
after soaking at maximum temperature for the appropriate 
time, the fatigue testing to be performed at temperatures 
appropriate to climb and descent manoeuvres. The same 
correspondent asked why continued use of alclad aluminium 
alloy sheet was allowed, stating that its fatigue strength 
could be about 40 per cent less than that of the bare forged 
material at any large endurance. It would surely be better 
to use anodic treatment and paint rather than a protective 
coating that seriously reduced the fatigue strength of the 
material. 

In his reply Mr. Heath said he was interested to hear 
the suggested compromise between the “short time” tests 
at elevated temperature and the “long time” fully pro- 
grammed tests, but he did not know if it would work 
only actual test results could show. Concerning the other 
point, he doubted that cladding could make so much 
difference to fatigue strength and thought that a compari- 
son of results using specimens containing typical joints 
would show that endurances were more alike than had 
been indicated. He wondered what paint the questioner 
would recommend for use at 300°C. 


DESIGN OF TRANSPARENCIES 

The discussion following Dr. Przemieniecki’s paper was 
opened by a speaker from the R.A.E. who said that 
although the lecturer had been dealing with the case where 
there was no heat flow on the inside of the transparencies, 
surely in any ordinary pressure cabin the heat flow on 
the inside would have an appreciable value? He pointed 
out also that transparencies had two rather unfortunate 
properties: (/) a very high resistance to heat flow which 
induced large temperature gradients, and (ij) low powers 
of radiation. He asked whether. in future aeroplanes 
designed to fly at very high speeds, dimensions of trans- 
parencies could be kept low enough to be cooled from the 
inside or, would it be possible to cool them from the 
outside? 

Dr. Przemieniecki said that for his paper he had con- 
sidered the simplest case only, i.e. the one where the heat 
flow on the inner surface was zero. This was perhaps a 
little unrealistic as there would obviously be some heat 
losses from this surface, but the principles remained the 
same in each case. The two unfortunate properties that 
the questioner had referred to were, of course, quite 
important, and in some practical calculations that had been 
made the temperature gradients had been found to be 
quite large. Some of the problems could be overcome by 
using a double panel construction and supplying cooling 
air between them. It was doubtful if it would be practical 
to supply cooling air to the external surface. 

Another speaker from the R.A.E. asked if the lecturer 
would agree that the choice of materials for transparencies 
was very limited because of all the conditions that had to 
be met. In the comparison between the fused silica and 
toughened soda-lime glass it had been shown that fused 
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silica had definite merits where the temperature change 
was high, but he felt that toughened soda-lime glass might 
still be considered favourably for use in civil air liners 
where any change of temperature would not be rapid. 

In replying to this speaker Dr. Przemieniecki said that 
he had been dealing with temperatures and temperature 
effects associated with Mach numbers greater than two, and 
therefore the use of plastic materials had been eliminated. 
He agreed that there might be rather more scope than he 
had indicated for use of toughened soda-lime glass, but 
pointed out the danger that at very high temperatures it 
would lose its toughening. 

Dr. Przemieniecki replied to a question on the use of 
transparencies and retractable blinkers in supersonic 
transport, by saying that this was not a practical solution 
and that he would prefer no passenger windows at all. 
Retractable shutters could perhaps be used on some special 
purpose windows, ¢.g. emergency hatches 

After commenting that he thought that the question 
of windows would be settled by the travelling public who 
obviously liked them, a speaker from the R.A.E. went 
on to say that he had been surprised to hear the lecturer 
assume zero heat transfer into the cabin. He stated that 
on fighters, a large proportion of the heat entering the 
cockpit entered through the transparencies. A usual cure 
was to increase the mass of air flow through the cockpit, 
but this usually involved increased air velocities over the 
transparencies and, hence, higher heat transfer coefficients, 
which was really no solution. He thought that the lee- 
turer's proposal of cooling the space between two panels 
would be difficult to achieve 

Dr. Przemieniecki replied that although the assumption 
of zero heat transfer on the inside surface was an unrepre- 
sentative case it could be used effectively for design pur- 
poses. This assumption led to conservative designs since 
cooling of the inside surface resulted generally in a 
reduction of the maximum thermal stress. He thought 
that no appreciable difficulty would be encountered in 
cooling the interspace, but due to a very high thermal 
resistance of transparency panels this method was not 
very effective in reducing temperature of the outer panel; 
consequently outer panels would have to be designed to 
withstand maximum temperatures 


SOME PRACTICAL ASPECTS OF KINETIC HEATING 
CALCULATIONS 

An aerodynamicist said that he was diffident on en- 
larging upon the complexities of the estimation of aero- 
dynamic heating rates. He congratulated Mr. Bore on 
his manipulation of the formulae and the very convenient 
form of presentation but requested him not to ask the 
aerodynamicists to produce charts showing the values of 
the parameters over a range of Mach numbers. 

A question was asked regarding the relation of enthalpy 
to Mach number. However, the complex nature of the 
problem meant that the answer could not be given briefly 
at this meeting. In later correspondence the questioner 
stated that he did not disagree with the form of equation 
(5), but felt that little error would arise in the estimation 
of recovery enthalpy (“ forcing enthalpy”) due to taking 
the coefficient of M* equal to 0-2. He agreed that the 
local temperature could be high in regions of much reduced 
flow speed, leading to changed values of the coefficient, 
but M, would be small in such cases, so that there would 
be little overall effect. He suggested that it was more 
convenient to use equation (1), thus avoiding any ambiguity 
over the value or definition of 

Mr. Bore replied that it was true that in many cases 
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little error would accrue from taking M? equal to 0-2, but 
it was quite easy to include the variations while preparing 
the charts, so there was little point in avoiding them. 

Larger discrepancies occurred between the strict local 
enthalpies and those obtained by taking the approximate 
coefficient when the local temperature was large and M, 
was also large. If M,=5 and M,=2:5, for example, the 
value of local enthalpy would be given about 2} per cent 
too large. 

The main reason for introducing the later forms of the 
energy equation was to facilitate calculations for regions 
where the local velocity was substantially different from 
free-stream velocity. In such regions he generally found 
it more convenient to estimate local Mach numbers than 
true local velocities which depended on the local tempera- 
tures, which were rather difficult to estimate when M, 
or U, were unknown. 


ADVANTAGES AND LIMITATIONS OF THE Use OF MODELS 

Concerning Mr. Sobey’s paper the same aerodynamicist 
said he thought that the idea of the Hoff study wing was 
am abomination. He had never liked the idea of infinite 
accelerations and he could only accept it as a purely 
academic problem. 

A speaker from Bristol Aircraft thought that Mr. 
Sobey had given rather a depressing picture of the use of 
models. It seemed that for a particular aircraft the use 
of models for determining the thermal effects did not 
have much future. Models had, however, been used with 
success in most fields of aerodynamics and structural re- 
search and they were valuable in that they provided a 
physical picture and a check on the mathematical theory. 

He thought that the Hoff wing could be used to obtain 
an experimental check on the theory, and that theoretical 
models could then be elaborated and developed and 
checked with physical models approaching nearer to real 
aircraft structures. 

For aeroelastic heating problems, simple models would 
be put in an air stream so that the heat transfer coeffici- 
ents could be easily deduced, and various characteristics 
could be measured and assessed and taken into account 
in theoretical studies. He felt there was a danger of having 
a wealth of theory with very little experimental evidence. 
Mr. Sobey replied that he agreed that more models should 
be tested to check with theory, but the difficulties were 
largely experimental. Unless the tests were directed at a 
specific simple structure it was very difficult to acquire 
any useful information. This was particularly so in the 
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case of rapidly transient phenomena, which should be 
avoided because of the limited test equipment available 

A speaker from the N.G.T.E. quoted the popular press 
as saying that future supersonic transports would travel 
at a Mach number of 2 or 3. As an engine man he said 
he would prefer a Mach number of 3 on the grounds of 
higher engine efficiency. He put forward the opinion 
that the future transport aircraft would be more important 
than the fighter aircraft. 

The Chairman here pointed out that the discussion 
should consider all aircraft in the very broad sense. 
Missiles must also be considered. 

A speaker from the R.A.E. supported the analogy with 
aerodynamics where any advances had come from small 
scale model tests rather than replica tests. At times, he 
said, the latter were necessary, but sometimes led to con- 
siderable difficulties. These difficulties were usually sorted 
out by the small scale model test when the problem was 
broken down into its component terms. 

He agreed that in order to deduce the most significant 
parameters, crude theories had to be constructed to isolate 
those parameters. 

Mr. Sobey replied that he did not really disagree but 
he felt there were times when replicas must be used. There 
was generally not the opportunity to do extensive tests 
from which to build up theory. A simple and single test 
was of great value. The small scale work piece was an 
instrument of scientific research, but did not obtain the 
useful structural information needed for preliminary 
design. 


CLOSING REMARKS 

In his closing remarks the Chairman thanked the four 
principal speakers and said they had given very useful 
information which would be of benefit to all those inter- 
ested in kinetic heating. He made two points in con- 
clusion. The first, a general point, concerned the vital 
importance of really understanding the test technique so 
that this could keep ahead of the actual aeroplane. 

The second point was that from the day’s discussion 
it might appear that they were on new and dangerous 
territory. This was not altogether true. There were a 
large number of complex aeronautical devices flying 
around with complete success and suffering from quite 
severe kinetic heating problems. Quite a lot was known 
about kinetic heating already, but he stressed that more 
knowledge was needed in order to design with extreme 
efficiency. 
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Icing Flight Development’ 


by 


B. T. CHEVERTON 


(Chief Flight Development Engineer, D. 


Introduction 

The aircraft designer must clearly ensure that his 
aircraft will fly in all climatic conditions normally to 
be expected. Over and above this, aircraft must be 
able to cope with a change in weather also. 

Since the last war, research into the problem of ice 
accretion on aircraft has done much to widen under- 
standing of this serious flying hazard. In the flight 
testing and development of ice protection systems for 
aircraft, however, one of the principal difficulties en- 
countered is that of finding natural icing conditions 
sufficiently uniform through which to fly the test air- 
craft. Because of this it has been found more profi- 
table not to seek natural conditions, but instead to fly 
in clear air at temperatures below 0°C, while spraying 
the test surface from a water spray rake fitted to the test 
aircraft 

It is nevertheless essential to prove that the simu- 
lated conditions are closely related to the natural 
conditions that would give the severity of ice accretion 
specified for test purposes, and it is desirable therefore 
to know the precise severity of the contributory natural 
factors to the icing risk, namely, ambient temperature, 
water concentration and droplet size. 


The Icing Risk 

Aircraft icing may occur in natural atmospheric con- 
ditions if flight is made in droplet cloud, the particle 
size of which is greater than, say, three microns, and 
when the forward facing surfaces are at a temperature 
of 0°C or below. When ice crystals are encountered 
in cloud a different situation may exist since, if the 
surface is not wetted by the thawing out of crystals, 
idhesion of ice will not take place 

It is possible to produce an icing condition by de- 
scending a cold aircraft into warmer but moist air, or 
to avoid an icing condition because of a surface kinetic 
heating effect associated with forward speed. An air- 
craft flying at 10,000 ft.. for example, on an I.C.A.N. 
day will be experiencing an ambient temperature of 

5°C, but at a true air speed of around 250 m.p.h., 
the wet air kinetic rise will be about 5°C and the icing 
risk will be associated with an ambient air temperature 
of only zero. The icing severity for a particular air- 
craft depends therefore on forward speed such that at 
a Mach number giving a wet kinetic rise of 40°C, the 
icing risk in droplet cloud may be nil. 


*Based on a lecture given to the Luton Branch of the Society 
on Sth December 1956. 


659 


apier and Son Ltd., Luton Airport) 


The Icing Severity Standard 

In flight, the icing severity experienced will vary from 
day to day and cloud to cloud, but it is necessary, for 
design purposes, to define a standard. This is given in 
A.P. 970 for military aircraft, and in B.C.A.R. papers 
for civil aircraft. 

Three degrees of icing conditions are promulgated. 
The most intense is Instantaneous, which may occur for 
distances of up to half a mile and with concentrations 
up to 5 gm./m*. “Intermittent ” conditions with con- 
centrations up to 2:5 gm./m*, may occur at intervals 
of three miles below 30,000 ft., and at intervals of 
twenty miles above 30,000 ft. over distances of up to 
three miles. The severity for Continuous Maximum 
conditions is shown in Fig. 1. The severity is tied to 
ambient temperature, water concentration and droplet 
size. It will be seen that the concentration at 30 
microns is half that for 20 microns, and that icing con- 
ditions may be expected at ambient temperatures from 
zero to —30°C 

It is by no means readily apparent which combi- 
nation of these severity conditions would prove most 
hazardous. A short period at high concentration might 
seriously impair engine performance, whereas a long 
period at relatively low concentration might destroy 
aerodynamic lift or increase weight and drag to un- 
acceptable limits 

It is partly to determine the relative seriousness of 
the several possible icing conditions that experimental 
flight testing is undertaken. Exploratory flight tests 
have been made at all conditions shown in Fig. 1 and 
also at Intermittent Maximum conditions. 


Formation of Icing Cloud 

The maximum amount of water vapour that can be 
present in an air mass without condensation occurring, 
depends mainly on the temperature of the air and be- 
comes less as the temperature decreases. If air is 
cooled, condensation therefore occurs when that maxi- 
mum is reached, the air then being saturated. 

In the atmosphere, cooling is most often due to the 
adiabatic expansion that occurs when air is moved 
upwards to a lower pressure level, so that the amount of 
water vapour present may be described in terms of the 
pressure level to which air must be carried up for it to 
become saturated and for condensation to water 
droplets to take place. This will be the condensation 
level or the altitude at which cloud forms. 

A cloud may form as the result of water vapour in 
the exhaust of an engine and this is an interesting study 
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Figure 1. Icing 
severity standard. 


but, from the icing consideration, an aircraft does not 
of itself condition the severity of the icing cloud 
through which it may subsequently fly. 

In clear air the condensation of water takes place 
into water droplets and not ice crystals, and if the 
temperature is below freezing point the droplets are said 
to be supercooled. Ice crystals do not form in clear air 
when saturation is reached, they are formed as the result 
of subsequent freezing of water droplets. 


Water Catch 

The icing risk has been considered briefly and it is 
seen that by flying faster the temperature range over 
which icing occurs may be reduced. It is not normally 
practicable to do this, although as a general maxim 
“the sooner out of cloud the better.” In flying faster 
however, the rate of catch is increased and it will be as 
weil now to consider severity in terms of rate of catch 
and ambient temperature. 

Clearly, catch rate will increase as water concen- 
tration and forward speed increase and it will be appre- 
ciated also, that large droplets have greater inertia and 
are more likely to be caught than small droplets. Catch 
increases also as frontal area is increased but that is 
inherent in design and may not be blamed on the 
weather. 


Figure 2. Water spray boom, showing 43 nozzle spray head 
and camera housing. 
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Rate of catch is equal to the product of the water 
concentration, the velocity of the aircraft, the frontal 
area of the section and, the catch efficiency. An esti- 
mation of catch efficiency is difficult and tedious to 
obtain, and this becomes a major difficulty in the 
estimation of heating requirements. It is nevertheless 
necessary for the designer to know the efficiency and 
hence, the rate of catch, and also the peripheral limits 
of catch in order to produce a protective system. 
Because of difficulty in assessing rate of catch or heat 
transfer, flight tests of a specimen are undertaken to 
confirm the suitablity of an anti-icing or de-icing design. 


The Test Aircraft 

At the Luton Flight Development Establishment 
some early tests were completed with an electrical mat 
mounted on the leading edge of a Tempest aircraft. 
The test section was flown for 100 hours in a variety of 
weather conditions at ambient temperatures down to 
—48°C. Surrounding paintwork was damaged by 
climatic conditions, but the weathering properties of the 
mat were good and further installations were made. 

A Lincoln aircraft was fitted with a water spray on a 
forward boom, more commonly known as “a proboscis.” 
Water was sprayed from a 43 nozzle spray head into the 
ducted spinner intake of a Napier Naiad propeller- 
turbine engine (see Fig. 2). These tests were undertaken 
partly to prove the effectiveness in icing conditions of 
Spraymats fitted to spinner and intake splitters. A 
camera fitted to the boom allowed a photographic’record 
of tests to be obtained—a scheme to fit a large mirror 
was rejected. 

Excellent results were obtained, the spinner remain- 
ing ice-free in moderate icing conditions with a heat 
input of less than 4 watt/in.* as is shown in Fig. 3. 
During initial tests an excessively high heat input 
resulted in the open-circuiting of one phase of the three- 
phase circuit, but complete protection continued to be 
provided, Fig. 4 demonstrates what may happen if 
protection is not provided; a typical cauliflower head 
appearance after 6 minutes in icing conditions is shown. 
It is also of interest that in conjunction with these 
tests an investigation was completed with a momentum 
separating type of intake designed to protect the engine 
in icing conditions. 

A Spearfish aircraft was used for spray nozzle 
development work. A tennis racket device was mounted 


Figure 3 (/eff). Napier Spraymat ice-free in moderate icing 
conditions with heat input approx. 4 watt/in.* 
Figure 4 (right). Spinner unprotected for 6 min. 
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in the spray of a single nozzle, and the area of spray 
was clearly established at several forward speeds and 
water flow rates. These tests, besides showing the 
clearly defined area of spray, produced the distinct types 
of ice normally encountered. 

Slush ice is dirty in appearance, appears to be wet, 


is easily dislodged and forms at relatively warm 
effective temperatures of, say, zero to —3°C. Glaze 
ice forms at temperatures around ~—7 to —8°C; it is 


relatively smooth, tends to be transparent but is dirty 
in appearance. Rime ice has a white milky appearance, 
and is opaque and forms at lower temperatures where, 
due to the rapid rate of freezing, the ice formation 
contains a good deal of entrained air. Quite obviously 
the appearance of ice varies with ambient temperature 
and does not readily fall into “ice-tight’”” compartments. 
For example, to describe ice formation that appears to 
be neither glaze nor rime we have coined the word 
“glime.” 

This slight digression from a discussion of test air- 
craft demonstrates the difficulty that will inevitably be 
experienced in classifying types of ice formation during 
flight tests, particularly when it is considered that the 
appearance of ice is considerably affected by lighting 
effect, as for example when the sun is shining on to the 
surface or, from behind the surface. I have seen ice 
apparently dispersed merely by a change in aircraft 
heading. 

The next aircraft used for icing research was a 
Vickers Viking; this aircraft was flown on icing tests 
for over 300 hours during the years 1951 to 1955. It 
had originally been intended to utilise the aircraft for a 
full scale de-icing investigation of mainplanes and 
empennage, and a Gypsy Queen auxiliary power unit 
capable of supplying sufficient electrical energy was 
installed. It was later agreed with the Bristol Aircraft 
Company to flight test a Britannia wing tip section in 
preference to other schemes. 

To cope with the size of section involved, an over- 
shoe to the aircraft fin was designed to house a variety 
of test sections. Cloud conditions were simulated by a 
water spray system with a spray rake mounted six feet 
in front of the test section. 

It is unfortunately not possible to describe at length 
the mosi interesting tests and results obtained with the 
Viking aircraft. It was, however, a splendid test 
vehicle. Flight testing is not always achieved in such 
comfort. 

Figures 5 to 8 show more recent installations on 
Lincoln aircraft, involving Beverley, Comet and 
N.A.C.A. 4414 type aerofoils. A number of heater mats 
have been added experimentally to the latter type 
section to provide easy clearance of run-back ice 


Water Spray System 

It has already been inferred that the water spray 
system is probably the weakest link in the test system. 
There is the obvious difficulty of housing a sufficient 
quantity of water for a test programme, and pipe lines 
have become frozen, even in well heated hangars. The 
introduction of anti-freeze was considered but it was 
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Ficure §. Beverley aerofoil. 
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N.A.C.A. 4414 type aerofoil mounted on Lincoln 
fuselage, with experimental heater mats. 


FIGURE 


Figure 8. Lincoln RF 402 showing instrument panel and 
section hot gas controls (numbers refer to Observer stations, 
a maximum of four Observers being employed). 
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Figure 9 (left), Poor water spray distribution during initial 
tests on Viking VL 229 test section. 


Figure 10 (right). Improved distribution after fitment of fair- 
ings to spray nozzles. 


thought that the Air Registration Board might not 
approve. 

There is then, the problem of passing a sufficient 
quantity of water to give high concentration but with 
atomisation adequate to give small droplet sizes. We 
have now reached the stage where we are able to 
achieve this, but it has only been possible as the result 
of lengthy testing in a cold wind tunnel and flight tests, 
as for example on the Spearfish aircraft. It is just not 
good enough to throw buckets of water at a test section 
or even to argue that conditions more severe than 
normal will adequately prove a section, particularly for 
de-icing tests. 

The next requirement is that the water spray shall 
be directed over a sufficiently large area, and from a 
distance sufficient to ensure that complete supercooling 
of droplets is achieved. Exploratory tests were necessary 
and Figs. 9 and 10 show the poor distribution obtained 


ee 


Ficure 11. Droplet sample from natural cloud. Volume 
median 11 microns. 
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at one time and the subsequent improvement resulting 
from fitment of nozzle fairings. 

In association with the water system there must be 
a nozzle air system, an additional air flow being 
necessary to achieve the required degree of atomisation. 
A large quantity of air is required, and it has been 
necessary to employ a cabin type supercharger for the 
job. To illustrate the importance, for another project, 
a Palouste turbo-jet engine is used for this purpose only. 

There is next the nozzle calibration necessary to 
ensure that correct droplet sizes are provided. This 
has entailed droplet sampling in flight and the pro- 
gramme has been of such magnitude that there was at 
one time a danger of mistaking a tail for the dog. 


Droplet Sampling 

If natural conditions are to be simulated, however, 
droplet size must be correct and it is most gratifying 
to be able to guarantee droplet size to within two 
microns. 

Figures 11 and 12 show firstly, a sample obtained 
in natural cloud and secondly, a sample with the artifi- 
cial spray system. Both are to the same scale. There 
is no suggestion that one is intended to simulate the 
other but the comparison demonstrates how far off 
from a true simulation we might be without precise 
calibration. 

The sampling apparatus, which may readily be in- 
stalled in other aircraft, consists of three independent 
components: an ice box in which oil slides prepared 
for exposure are stored, a sampling device for exposing 
the slides to a droplet stream for a controlled period, 
and an adapted Newman Sinclair camera fitted with a 
macro lens. Fig. 13 shows the observer station. 

The type of oil used is important, it being necessary 
to obtain an even film with good photographic pro- 
perties. Because the oil film takes and hold samples 
most efficiently at temperatures around 0°C, and be- 
cause it is essential to keep the slides clean, they are 
stored in a dust-proof ice box. An appreciable rise in 
temperature would cause the oil to run and destroy the 
even depth of film, whereas a decrease in temperature 
below 0°C may cause droplets to freeze on impact and 
form ice crystals. 

There would be a risk of condensation from 
observers’ breath if care were not taken. We have 


Figure 12. Droplet sample in simulated cloud. 
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FiGuRrE 13. Droplet sampling station in Viking VL 229. 


estimated the risk of contamination through perspiration 
as slight, but there is a real danger of critical evapora- 
tion of droplets from the oiled slides if too long a period 
elapses between the taking and photographing of a 
sample 

The over-exposure of a slide will cause too many 
droplets to be caught with the result that adjacent drop- 
lets coalesce, causing an exaggeration in the subsequent 
calculations. A correction factor determined experi- 
mentally may then have to be applied. 

To control the exposure time a simple pneumatic 
sampling device was used at first in the form of an 
aif-operated piston by which single slides are moved past 
an aperture in a tube that is raised by the observer to 
project into the droplet stream. The mechanism is 
designed to ensure that while the aperture is exposed 
the slide passes it once only. 

After sampling, areas of approximately one square 
millimetre are magnified and photographed on 35 mm. 
film; three or four shots are obtained to give representa- 
tive views, although, with a uniform spectrum, variation 
in mean droplet size calculated from each photograph 
should be negligible. The globular appearance by which 
droplet images are identified as droplets or otherwise 
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Ficure 14. Curve used to calculate volume median droplet size. 


FLIGHT DEVELOPMENT 


FIGURE 15 Catch 
limits associated with 
droplet size of (a) 20 
microns and (6) about 
55 microns 


is well illustrated in Fig. 12, together with the clear 
definition required for measurement. 

The film is next projected on to a suitable screen 
and by super-imposing a calibrated scale, the droplets 
are then counted into size bands of 5, 10, 15, 20 and 
so on, microns diameter. It is next required to determine 
the mean droplet size but before describing the method 
some explanation of its significance may be useful. 

The essential feature of droplet size that affects the 
amount and distribution of ice accretion is the droplet 
weight, which is proportional to the volume. However, 
since the droplet size is not uniform throughout any 
cloud it is necessary to establish an expression of drop- 
let volume in terms relevant to the icing severity pro- 
duced by the particular droplet size composition of the 
cloud. The relevant expression is not an average 
volume but a median volume which is stated in units 
of droplet diameter at which, for the example concerned, 
the total volume of droplets contained in the spectrum 
of diameter less than the median is equal to the total 
volume of droplets of larger diameter. 

This value, known as the volume median, is usually 
given when natural icing conditions are specified in 
meteorological data, and is determined from a sample 
by constructing a curve to data from the count. A 
curve, as shown in Fig. 14, is plotted to co-ordinates 
of nominal droplet size in microns and aggregate per- 
centage of the total weight of droplets counted. The 
droplet volume median is then determined by reading 
off the nominal droplet size at the position on the curve 
corresponding to 50 per cent by weight of the total. 

The critical nature of droplet size is shown clearly in 
Fig. 15. 


The Protective System 

he icing risk and the problem of providing artificial 
simulation having been considered, some reference 
must now be made to the protective system to be 
emplov~u, although a comprehensive list of the various 
systems that might be employed will not be attempted. 

One form of protection is to use a mechanical ice- 
breaking device on the ieading edge of an aerofoil 
surface. An inflatable rubber tube might be used so 
that when expanded at intervals it would serve to break 
off ice adhering to it. The disadvantages are plain in 
that it may not always clear satisfactorily and that in 
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high speed aircraft the presence of an inflatable bag 
is not compatible with aerodynamic efficiency. 

Another method of protection is to spread on the 
surface to be protected a film of some liquid, such as 
alcohol, which serves to lower the freezing point when 
mixed with water. Even mixing is essential but difficult 
to control and, apart from being messy, a large amount 
of fluid must be carried with consequent reduction in 
payload. 

Special pastes and even warm engine oil have been 
used for anti-icing and, although useful for limited 
application, as is for example anti-freezing grease on 
control actuators, these methods do not recommend 
themselves for use on aircraft lifting surfaces. 

There are two methods readily employable and both 
were tested on the Viking test aircraft—they are hot 
gas anti-icing and electrical de-icing methods. 

Anti-icing presumes that heating becomes effective 
before encountering an icing hazard. The entire surface 
on which water may catch or run back should be main- 
tained at a temperature above 0°C. De-icing presumes 
that ice may be allowed to form, being subsequently 
cleared in a selected sequence from separately heated 
mats. 

Anti-icing and de-icing systems may be compared 
but the overall suitability of one system as against the 
other is tied up with installation requirements, inter- 
changeability, liability to damage and so on, not exclud- 
ing cost. 

For flight test experimental work, however, it must 
be accepted that during an immediate flight programme 
we are mainly concerned with the mechanical reliability. 
endurance and functional ability of a component under 
test. 


Flight Testing an Anti-Icing System 


The hot gas anti-icing system as flight tested will 
now be considered. Hot air must be ducted to the test 
section which must be set at representative incidence 
to the effective air flow. The heat input in terms of air 
mass flow and temperature differential must be 
measured, further instrumentation will be required to 
allow for reproduction on full scale tests and may take 
the form of surface thermocouples and reference for 
internal structure temperature and air flow distribution. 
The hot gas may be supplied from engine exhaust or 
compressor tapping, or from a separate ram source 
through a combustion heater or other heat source. 

In an experimental installation the entire anti-icing 
system may not be standard. For example: on the 
Viking aircraft it was necessary to employ an electrically- 
heated element to replace a combustion heater. Correct 
simulation of the actual section was essential, particular 
care being taken neither to exclude nor introduce cold 
spots. 

In order to determine the lift coefficient at particular 
section incidence, pressure plots have sometimes been 
obtained in flight and, to determine section incidence, 
yawmeters have been used. 

Anti-icing systems may be classified into two kinds, 
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HEAT REQUIREMENT CHU/ HR/FT SPAN 


AMBIENT TEMPERATURE °C 
Ficure 16. Anti-icing heat requirements. 


either “wet” or “dry.” A wet system is designed to 
keep a protected surface ice-free while allowing water 
not evaporated to run back and form on unheated 
portions aft where it will not necessarily be an em- 
barrassment. A dry system presumes complete evapora- 
tion of all water caught and therefore requires a larger 
heat input. 

In determining the heat requirement, therefore, run- 
back is the main guide, although at very low ambient 
temperatures associated with low rates of catch ice may 
form directly on the heated portion. 

The testing technique is to set initially to a high 
heat input and then to reduce progressively, maintaining 
inlet temperature or air mass flow constant as pertinent 
to the system being tested, working always within the 
temperature limits for the type of aerofoil construction 
under test. At the heat input just insufficient to evapor- 
ate, or otherwise clear, all water running back to the 
unheated portion, a streak of ice will be detected. 

It has been observed that run-back forms initially 
behind a region of low heat input, as for example. 
multiple skinning on a wing. Run-back will also occur 
in a region of high water catch due to poor spray nozzle 
configuration or local roughness of the skin as, for ex- 
ample, when external thermo couples are fitted. 

A task for the Flight Observer is to detect this first 
instance of run-back and this calls for continuous view- 
ing and frequent change in aircraft loading to provide 
variation in solar lighting. In fact when working to 
close limits, we have observed that run-back ice clears 
more rapidly from a surface on which the sun is shining. 

For a particular icing severity condition tied to for- 
ward speed and section incidence, the maximum heat 
requirement may be determined. Theoretically, as 
shown in Fig. 16 at curve A, the heating intensity 
required at Continuous Maximum severity and zero 
speed is a maximum at 0°C. The heat input required 
is a compound of heat to raise the water temperature 
from ambient to surface temperature, heat to maintain 
a surface temperature and heat to evaporate the water 
caught. 

The maximum temperature at which icing may 
occur depends on forward speed and curve B shows 
the heat input actually found to be necessary. The 
dotted part of curve B is difficult to obtain, there being a 
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Figure 17. Ambient 5°C. Concentra- Ficure 18. Ambient 
tion 1-18 gm./m*. Encounter time 13 min. tion 1°14 gm./m*. 
“Cruise” configuration. “Cruise” 


very small ambient temperature difference between that 
giving zero heating requirement and that giving maxi- 
mum heating requirement. 

An interesting feature is the difference between A 
and B; briefly the reason is to be found in: 


(a) That not all of the water caught needs to be 
evaporated since some blows off the section and 
does not subsequently catch. 


(b) That the boundary layer becomes heated, thus 
allowing some run-back water to flow on to the 
unheated portion from which further evapora- 
tion takes place. 


In practice, it has also been found that water may 
run back in rivulets so that the entire heated surface area 
behind the catch area is not wetted even at the point of 
breakdown. This effect, and the practically impossible 
task of determining, firstly, transition points in wet air 
from laminar to turbulent flow, and secondly, the result- 
ing heat transfer coefficient, makes compilation of curve 
difficult. 


Figure 20. Ambient —15°C. Concentra- Figure 21. Ambient 


tion 0:98 gm./m*. Encounter time 20 min. tion gm./m*. 
“Climb” configuration. 


FLIGHT DEVELOPMENT 


“Climb” configuration. 


74°C. Concentra Figure 19. Ambient —10°C. Concentra- 


Encounter time 17 min. tion 1:10 gm./m*. Encounter time 15 min. 
configuration 


‘Climb” configuration. 


At — 30°C, however, there is no bonus since break- 
down occurs as the result of ice freezing on the heated 
portion and there is a small margin of heat input only 
between that necessary to keep the leading edge ice-free 
and that required to Keep the heated position ice- 
free. The first instance of run-back determines the ice- 
free heat input. This may, however, be an extravagant 
intensity, and it has been shown that a reduction in heat 
input value at 10°C may be tolerated without 
embarrassment 

As heating is further reduced the ice build-up in- 
creases and moves towards the leading edge, ultimately 
forming on the heated portion. It has been found that 
once this condition has been reached only a small further 
reduction in heating intensity is required to cap the 
leading edge. The difference in heating intensity 
between the ice-free and ice-capped state is therefore 
greater at —10°C than at —30°C. 

Figure 16 is of further interest since it shows that 
at constant indicated forward speed and section inci- 
dence, the higher concentration at —8°C ambient 
produces a more severe condition than the lower 


—30°C. Concentra- Figure 22. Ambient —10°C. Concentra- 
Encounter time 18 min. tion 1:10 gm./m*. (Simulated reduced heat 


condition with 10 min. at final heat input.) 
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concentration (and temperature) at —30°C. The design 
condition then for the particular section under review 
would be -8°C, and the heat to be supplied is the 
theoretical value multiplied by the ratio B to A. 

A number of photographs taken at various icing 
conditions are shown in Figs. 17 to 22. Tests were at 
Intermittent Maximum severity although heat input was 
that found adequate to cope with Continuous Maximum 
conditions; this accounts for the run-back ice ex- 
perienced. Attention should be given to the following 
features : — 

(a) The location of run-back ice with respect to the 

heated portion. 

(b) The chordwise extent of run-back ice. 

(c) The formation of rime ice behind run-back, due 

possibly to turbulence caused by the run-back 
deposition. 


Flight Testing a Cyclic De-icing System 

A cyclic de-icing system has essential 

features : — 

(a) Continuously-heated parting strips are provided 
along the leading edge and between primary 
heated areas. 

(6) A number of separately heated mats operate 
to a time sequence best calculated to ensure 
adequate shedding without danger of bridging 
the parting strips. 


PRIOR TO COMMENCEMENT AFTER FIRST APPLICATION OF 
OF CLEARANCE CYCLE TO ALL AREAS 


ABYER SECOND APPLICATION OF MEAT | AFTER THIRD APPLICATION OF MEAT! 


Figure 23. De-icing clearance cycle. 
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Figure 24 (left). Minimum acceptable leading edge gap width 
FiGure 25 (right). Clean shedding from primary areas. 


Clean shedding is helpful but not a criterion, since 
ultimate shedding with a cyclic system may be achieved 
by re-selecting mats on after flying out of cloud as 
shown in Fig. 23. 

A de-icing system must inevitably be a compromise, 
since it must cope with the most severe condition likely 
to be encountered without bridging of the leading edge. 
and it must operate in low severity conditions without 
producing run-back. Fig 24 shows what was taken as a 
minimum acceptable leading edge gap width at - 16°C 
ambient temperature. The condition was excessively 
severe, water concentration being 1-12 gm. /m*. Leading 
edge heating intensity was 14°5 watt/in. with 10 
watt/in.* to primary areas. The photograph shows the 
condition at the begining of the fourth primary cycle. 
that is after nearly 18 minutes at this severity condition 
Fig. 25 shows the subsequent clean shedding obtained 
from primary areas Bl and B2. 
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INTENSITY 


HEATING 


AIR °C 
Figure 26. Leading edge breaker strip heating requirements. 
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FiGuRE 27. Relationship between heating intensity and primary 
area shedding time 


During flight development of a de-icing system the 
following requirements must be established: 


(i) The heating intensity and breaker strip width 
particularly at the leading edge. Fig. 26 shows 
a typical requirement curve. 

(ii) Whether variation in this requirement justifies 
a multiple-way heating intensity switch, to 
cope for example with ambient temperature 
range of say 0 to ig oe 10 to 174°C 
and 174 to —30°C, requiring heating in- 
tensities of 8, 12 and 15 watt/in.° respectively. 

(iii) The optimum length of complete operation of 
cyclic sequence, this being allied to the time 
that each mat will be on, and equally impor- 
tant, the time that each mat will be off. Fig. 
27 shows a typical plot obtained of shedding 
time against ambient temperature at constant 
watts loading. The “heat on” time determines 
the “heat off” time and, hence indirectly the 
sizes of mats required and total wattage for 
cyclic area. 

(iv) The optimum heating intensity for primary 
areas. Here the observer comes into his own, 
it is not a matter one can prove subsequently 
with curves and figures. 

(v) The usefulness of graduated heat input from 
front to rear of the section must be determined 

(vi) Whether secondary areas are necessary 


Figures 29, 30 
and 3] Inade 
quate heat input 
Effectiveness of 
continuous heated 
areas on Spray- 
mat. Water con 
centration 0:24 
gm./m*. Ambient 
10 
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Figure 28. Control panel for electrical thermal de-icing 
experimental tests on Viking VL 229. 


There are, of course, other considerations. A 
considerable amount of rig-testing will have been 
necessary and it must be accepted that some testing 
may more usefully be done in closely controlled and 
observed laboratory conditions before flight testing 
begins, particularly when it is necessary to determine 
trends toward optimum configuration. 

With an electrical system, heating intensities are 
closely controlled and readily predicted. A selected 
watts loading remains constant at all forward speeds. 
altitude, ambient temperature and distances from the 
heat source. There is not the same need therefore for 
extra instrumentation to ensure simulation of heat out- 
put in full-scaled tests. This makes life easier for the 
observers, although it will be sobering to look at the 
special test cyclic control panel (Fig. 28). There was in 
this installation plenty for at least one observer to do. 

The observer crew requirements were largely as for 
the anti-icing investigation. The auxiliary power unit 
was operated by a Flight Engineer who had the least 
comfortable ride, being closeted in a sound-proof cabin 
with a 150 s.h.p. piston engine. His job was to maintain 
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A.P.U. r.p.m. constant thereby maintaining constant 
spray nozzle pressure and alternator output. 

No. | observer was able to view the section through 
a slit in the cabin roof (Fig. 24). Additional viewing 
was provided from a blister in the entrance door (Fig. 
25). Later in the programme a blister was fitted to the 
starboard side giving yet another view of the test section. 
No. 2 observer controlled the water flow rates and nozzle 
air pressures, selecting the concentration and droplet 
size required at each test level. No. 3 observer con- 
trolled watts loading to the heated section and was also 
chief permutator of the cyclic panel. No. 4 observer 
had the dual role of assistant photographer and droplet 
sampler. 

It will be instructive to observe the penalty to be paid 
at Continuous Maximum severity conditions of supply- 
ing insufficient heat as is shown in Fig. 29. Comparative 
views of satisfactory operations are shown in Figs. 30 
and 31. 


Recent Experience 

In conclusion some technical data relevant to a 
more recent clearance undertaken with the Lincoln 
aircraft, that for the Comet 4 is of interest. Twenty-four 
flights were made, involving a flying time of about 
46 hours. 

The flight conditions of climb and descent were 
demonstrated and the climb was shown to be the more 
severe; two combinations of climbing speed and 
incidence were used. 

It was also required to cover an air temperature and 
altitude operating range of from I1.C.A.N. +25°C to 
LC.A.N. —5°C. Demonstration was concentrated on the 
LC.A.N. + 25°C condition. 

The heat outputs and working temperatures to be 
expected on the Comet were set up on the Lincoln as 
shown in Table I. The water flow rates required to 
meet the maximum Comet rates of catch were employed. 

Climb on the Comet 4 is at 7,350 engine r.p.m. and 
a working air temperature of 200°C is maintained for 
anti-icing purposes by diluting compressor air by cold 
air bleed as necessary. The air flow available for anti- 
icing purposes becomes less as altitude is increased and 
varies from 0-088 Ib./sec. per foot run at —5°C on an 
LC.A.N. —5°C day, to 0:038 Ib. sec. per foot run at 
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TABLE 1 


Inlet Heat input 
Regime Ambient chu./hr./ ft 
run 


Climb ICAN + 25°C -—10 9980 
Descent ICAN + 25°C —10 7470 
Climb ICAN + 25°C —20 8850 
Climb ICAN + 15°C — 30 9100 
Descent ICAN + 15°C 30 6840 


Heat input is defined as =m Cp (tgi 
where m=mass flow of hot gas in Ib./sec 
where ftgi=inlet temperature °C. 
where to=ambient air °C. 


30°C on an LC.A.N. + 15°C day; the latter case im- 
plying a high altitude. 

Descent is at 6,000 r.p.m. and satisfactory perform- 
ance was demonstrated. At this lower r.p.m. a working 
temperature of 200°C is no longer obtainable, and both 
temperature and air flow are dependent on altitude, 
being 167°C and 0-077 Ib./sec. per foot run on an 
LCA.N. —S°C day at —S°C, and 149°C and 0-036 
Ib./sec. per foot run at — 30°C on an LC.A.N. + 15°C 
day. 

It was demonstrated that with the design heat inputs. 
anti-icing performance was satisfactory at B.C.A.R. Con- 
tinuous Maximum severity conditions. 

Based on the Lincoln tests, the A.R.B. agreed that 
airframe icing approval would be granted for the Comet 
4 in the event that the same dry air surface temperatures 
were achieved in flight as were associated with the satis- 
factory heat input conditions during Lincoln tests 
This was achieved on the Comet as evident from the 
clearance obtained. 


Footnote 

The Lincoln has since been converted to carry wing 
and tailplane sections of the Sud Aviation Caravelle, 
and tests to Continuous Maximum and Intermittent 
Maximum severity have been satisfactorily completed 

The aeroplane is now being converted to carry a 
wing tip section of the Avro 748 aircraft; this will be 
followed by testing of the A.W. 650 tailplane. 
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Iterative Solutions of the Equations for Plane Oblique Shock Waves 


A. R. COLLAR, M.A., D.Sc., F.1.A.S., F.R.Ae:S. 


(Department of Aeronautical Engineering, University of Bristol) 


1. THE PRINCIPAL EQUATION 
If a plane oblique shock wave, inclined to the free 
stream at the angle «, is produced in two-dimensional 
supersonic flow of Mach number M by (for example) a 
wedge which deflects the flow through an angle 4, the 
equation connecting these quantities may be written 
2 (M? sin? -—1) (1) 
an 6=cot ¢ = 
(y+ 1) M?—2 (M? 
In this form, 6 is given explicitly when M, « are fixed 
Similarly, we may obtain M explicitly when <, 6 are fixed; 
equation (1) may be written (see, for example, Liepmann 
and Puckett'”, Equation 4.27) 


l y+1) sin 4 
M- 2 cos 


Unfortunately, it is rare to meet a problem in which 
M, <= or ¢, 6 are given. It is nearly always the case that 
M, 6 are given and that it is required to find < and to derive 
pressure ratios, and so on. Equations (1) and (2), for 
given M, 4, yield a cubic equation for the determination 
of <: one of the roots is physically inadmissible, and of 
the other two, only the one corresponding to the weaker 
shock wave appears in practice. Because of the difficulty 
of solving a cubic and obtaining the relevant root, it is 
often the case that equation (1) is used to draw a chart or 
set of charts from which < may be read off when M, 6 
are given. However, such charts are often either too large 
to be convenient, or too small for accurate reading. The 
writer suspects that there may be many who, like himself, 
have neither the inclination nor the eyesight for reading 
such charts, even if one remembers where the chart is to 
be found. Accordingly, he feels that it is worth noting 
an arithmetical method of solution of equation (1) which 
automatically yields the correct root of the cubic equation 
determining <«. The method proposed is one of iteration; 
it will yield results to any desired degree of accuracy. 


2. THE REQUIRED SOLUTION 
In equation (1) write 
x=cot (3) 


and re-arrange the equation to read 


tan A | 


y+1 
+(4 M?+1)tan | 
We then adopt the iterative process 
B 
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(ag? 


x= 


where 
A=M?—-1, 


— y+! 4 ¢- 
B= M* tan 06, 


~ 


(6) 


| 
c=(% M?+1)tané, 


Note that A, B, C are positive for all M, 6. We begin the 
iteration with x, large. If, in fact, the initial choice is 
infinitely large, the next value is evidently A'/*, and this 
is thus a suitable choice for x It is easy to see that this 
is greater than the root sought, since 


A'/?=(M?—1)'/?=cot My 


where u is the Mach angle. Since « is necessarily greater 
than yn, it follows that A'/? is greater than cot ¢. However, 
it is, of course, not necessary to choose x,=A'/?, if a closer 
approximation is known or can be inferred. It is, in fact, 
not necessary to begin with x, large: the iteration (5) 
converges to the greatest positive root of the cubic implied 
by (4) (if a positive root exists; that is, if an attached 
shock is possible) provided only that the initial value 
x, is greater than the second (positive) root. If an attached 
shock cannot exist, then after a number of steps x,,, 
becomes imaginary. 


3. THE CUBIC EQUATION 
If we rationalise equation (4) we obtain 


f (x)= + Cx*—Ax+(B—AC)=0, (7) 


where A, B, C are given by equation (6), so that 


B—ac= (7 M?+1)tan , 


It is to be noted that the coefficient of x* is positive, that 
of x negative, while the constant term is positive, for all 
values of M, 6. This implies certain properties: in par- 
ticular, that the cubic has one, and only one, real negative 
root. To prove this, we observe that f (x) changes sign in 
the range —oO< x < 0; there will therefore be either one 
or three real negative roots. But if there are three, then 
f(x) must have its maximum and minimum values for 
x negative. This, however, is not so, for the constant term 
in f (x) is negative. There are thus not three negative 
roots: there must be one and one only. If the remaining 
two roots are real, they must be positive. 

This, of course, accords with the well-known fact 
(demonstrated, for example, by the hodograph of the flow) 
that for given M, 6 there are two shock inclinations ¢ 
which are physically possible; of these, it is always the 
weaker shock (smaller <, larger x) which appears to be 
stable and which is observed in practice. 
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4. EXAMPLES OF THE ITERATIVE PROCESS 
As a first, simple example, let M?=10, tan 6=0-1. 
Then, with y= 1-4, we have 


A=9, B=12, C=1°3. 
These numbers give, for the roots of the cubic (7), the 
values —3-73267, 0-03350, 2-39917. The iteration pro- 


cedure is, by (5) 
12 


= 
{ x, +13 


If we begin with x,=3, the successive values of x,, as 
computed with a slide rule, are 


(9) 


x, x, xX; x, 
3 2-49 2-42 2-40 2:40 


Similarly, if we begin with x,=0-5, we obtain 


Xo x, x, xs x, x; 
153 218 235 23 24 24 


This is an example of an initial choice lying between the 
positive roots. However, if we begin with x,=0, the 
process will not converge; one step gives x, imaginary. 
This example shows that, in the absence of any know- 
ledge of the required value, it is safest to begin with 
x,=(M?—1)'/?, since this is known to be greater than the 
root sought. 

If a desk calculating machine and table of square roots 
is available, a calculation can be carried through very 
rapidly without recording any of the steps in the iteration. 
For example, let 


M=2:-37, 6=15° (tan 6=0-26795). 
Then the iteration formula is 


10-1445 


We begin by looking up 4-6169'/? in the Tables, setting it 
on the machine and adding 2-0740; use the sum to build 
up 10-1445, set up the multiplier found, subtract from 
zero, and add 4-6169. The square root of the answer is 
again found from the Tables and the procedure repeated. 
In this way we find the required root to be 1-2516. 

Finally, as an example of a flow in which an attached 
shock wave cannot exist, let 


*=10, tan 6=1-0. 


120 
49— 
Fatt { x, +13 


Then 


This yields the sequence 


x, x, Xs x5 x, 
1-225 0-751 0-523 0-355 0-121 imaginary 


5. CONVERGENCE 

It is usually easier to use iterative processes of the kind 
proposed than to demonstrate rigorously their convergence 
or to find the limits to which they converge. The follow- 
ing discussions are not claimed to be rigorous, but they 
may help in giving some insight into the convergence or 
otherwise of the process. 
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5.1. Analytical Treatment 
Let the system be such that an attached shock can 
exist. We know from the previous discussion (Section 3) 
that f(x)=0 has three real roots, one negative. Let ¥ 
be the greatest (;,ositive) root, X, the second positive root. 
Now the iterative process is 


B 1/2 


and the desired solution is ¥,, where 


B 1/2 
x,= 
X,+C 


If now a step in the iteration begins with x,, where 
x, >X,, 


it follows, since A, B, C are positive, that 


or, in view of (12) and (13), 
Bass 


Also, in view of (12), 


Now since x, > X, and X, is the greatest root of f(x), 
it follows that f (x, is positive, whence 


From (14) and (15) we see that 
Xn 


and each iteration therefore yields a closer approximation 
to X, than its predecessor. 
If we begin with x, < X,, we obtain similarly the final 
result 


provided f (x,) is negative; that is, provided x, lies between 
X, and X,. If x, < X,, then f(x,) is again positive, and 
the iterative process clearly diverges. 

Finally, if no positive roots exist, equation (14) is no 
longer applicable, but equation (15) is still true, since 
f (x,) is now positive for all positive x. Hence the iteration 
procedure eventually leads to a negative value for x,., ,” 


5.2 Graphical Treatment 
We can regard the iteration process (5) as an attempt 
to find the intersection of the two lines 


y=x ‘ . (16) 


The curve (16) passes through the first and third quadrants. 
As regards (17), it is (in view of the radical) symmetric 
about the x-axis. In the third quadrant, it has a branch 
asymptotic to y=—A'/?; g(x) is increasingly negative as 
x increases and the curve tends to the asymptote x= —C. 
Clearly (16) and (17) will have one intersection in this 
quadrant. 


—— 
ie 
B B 
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For—C < x <(B—AC)/A, g(x) is imaginary. It becomes 
real again for x =>(B—AC)/A, and in the first quadrant 
rises from zero with initial slope infinite, curves over and 
tends ultimately to the asymptote y=+A!'’*. Differentia- 
tion of equation (17) shows that the slope decreases 
monotonically from infinity when g=0 to zero for x 
infinite. 

Figure 1 shows the curves (16) and (17) for the first 
example of Section 4, that is, M?=10, tan 6=0-1. The 
iteration procedure is equivalent to following the paths 
PORST ...«@ to the root X,. 
For if we begin with x infinitely positive and substitute 
in (17), y=g(co)=A!'/*,. Substitution of this value of 
y in (16) gives the point P with x also equal to A'’*. This 
value of x in (17) gives Q, which in (16) gives R, and so 
forth. The two paths POR ....P,Q,R,... . shown 
correspond to the two iteration tables of Section 4. 

It is also clear that if the curve of g(x) does not inter- 
sect y=x in the first quadrant, an iteration will give a 
zigzag path alternating between the curves until a vertical 
step misses the curve g(x): that is, an imaginary value of 
g is obtained. 

Finally, it is also clear that if g(x) and y=x approach 
tangency (whether they just intersect or just miss) a large 
number of small steps in the iteration will occur. 


6. OTHER ITERATIVE PROCESSES 

The iteration process (5) is, of course, not unique; it is 
however, the most satisfactory the author has found. 
However, a few remarks on other formulae may be of 
interest. 

Equation (1) may be transformed by the substitution 


z=M? «—1 
to the form 


M? “1/2 
€ M?—z,) = 1) tan 0. 


If this process is begun from z,=0 (which, like x,=A'!? 
earlier, corresponds to ¢,=,) it converges to the required 
result. However, each step requires two substitutions of 
z,, so that the process takes longer, even though the 
number of steps required may be less. 

An iteration of interest is derived from (5), but with 
the negative sign, namely, 


x,+C 


If we begin with x infinitely negative, we find the next 
value to be —A'/*. Fig. | then demonstrates the quasi- 
spiral path P,Q.R,S,T, . .. . which converges in this case 
on the negative root: 


x, Xs x x 


=—401 -366 —3-73 —3-79 


When the negative root has been found, it can be 
extracted from the cubic, and the positive roots, if any, 
then found. However, this iteration does not always 
converge; it does so only if the negative root lies at a 
point for which 


If the (negative) slope exceeds unity, the process will 
clearly diverge. For example, if we attempt to solve (11) 
with the negative sign, beginning with x,=—3, x, is 
imaginary. However, we can proceed in reverse; (11) 
gives x, in terms of x,,, as 


120 


x,= — 13, 


and if we take n=0, this gives the reverse sequence 


—13-75 -—13-67 -~13-67 


| | y 
ASYMPTOTES 


R- —— 


PART OF FIRST QUADRANT, 
ENLARGED 


FiGure |. 
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We can also use this reverse iteration to obtain con- 
vergence toward the second positive root. If we revert 
to equation (9), for example, and write it as 


12 
z,= 1 3, 
we find, if we begin with x, =2, 
1-1 0-24 0-042 0-0336 


0-03350 


while the same result is achieved in two steps from an 
initial choice of zero. 

However, these iterations are principally of academic 
interest, as an underlining of the conclusions of Section 5 
on convergence. The important requirement for an itera- 
tive process is, in this case, that it can be used with 
confidence in all circumstances and that it should lead 
to the root sought, if such a root exists. The iteration (5) 
satisfies these requirements. 


7. COMPLETION OF THE SOLUTION 

The iteration process determines x=cot ¢; when this 
has been found, the Mach number normal to the shock, 
namely, M sin ¢, may be calculated and the pressure and 
density ratios, and so on, then determined. However, with 
a little manipulation, formulae involving x directly may 
be established; we have 
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yM? tan 6 
x+tand 


tané 


(19) 


(20) 
M? +1) tan 6 


and the temperature ratio can be found from these. For 
example, using the illustration of Section 4, we have, with 
y=1-4, 
M?=10, tan 6=0-1, x=2-39917. 

Then 

Ps T, 1-5602 

— =1-5602; —* =1-3705; — = — 

Pr T,  1-3705 

Lastly, if M, is the Mach number after the shock, it is 
given explicitly by the energy equation in the form 


= 1-1384. 


or with y=1-4 
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The Effect of Initial Lateral Tension or 


Slackness on the Deflection of Long, 


Thin Plates under Normal Pressure 


H. L. COX, M.A., F.R.Ae.S. and D. DAWSON, G.I.Mech.E. 
(The National Physical Laboratory) 


HE BASIC THEORY of the deflection of long, thin 

plates under normal pressure is extended to cover 
the effects of initial lateral tension or slackness in the plate 
and the effect of deformation of the supports at the edges 
of the plate under the lateral tension developed. The 
procedure is illustrated by application to a case of a plate 
under initial tension with its edges held apart elastically, 
and the computed values of the strains are compared with 
values measured experimentally. The application to an 
initially slack plate is also illustrated by a diagram. 


1. BASIC THEORY: THE DEFLECTED FORM 
The deflection (w) of a long plate under normal 
pressure (p) is governed by the von Karman equation, 


4 


where T is the tension per unit length of the plate across 
its width (in the direction Ox), and El is the modulus of 
section per unit length of the plate for bending in this 
direction. For moderate values of w (say up to one-tenth 
of the width 6) T may be regarded as independent of x, 
but its constant value may depend upon initial slack or 
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tension in the undeflected plate and upon the stiffness of 
the side supports which hold the edges of the plate apart. 

The von Karman equation (1) may be rewritten in 
the form 


1 d‘*w d*w pb (1a) 
(aby? dx* dx? T 
where w’=w/b, x‘=x/b and 2?=T/El, and the solution 
may be written in the form 


w’=(pb/T) (x/b, ab) 


The general form of v is 


Y= —4(x/bP +A (x/b)+B+C cosh sinh zx (3) 


and the values of the arbitrary constants A, B, C and D 
are determined by the four conditions of deflection and 
rotational support, two at each edge of the plate. 

The form of the function y is therefore completely 
determined by these edge conditions and by the value of 
ab, and the effect of elastic restraint against movement of 
the edges of the plate in the plane of the plate is entirely 
comprised in the value of ab. 


% P. 
Pa =] 
Pi 
4 
4 
T, 
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2. BASIC THEORY: DETERMINATION OF THE MEMBRANE 

STRESS 

If the edges of the plate were held at a fixed distance 
apart, the centre line of the plate in the deflected form 
would be extended by a fraction proportional to (w/b)?’. 
Therefore this extension takes the form (pb/T)* 4 (2b), 
where the function » of (zb) depends otherwise only on 
the condition of support of the plate edges against deflec- 
tion normal to the plate and against rotation. The 
condition of support against movement of the edges of 
the plate in the plane of the plate is a separate issue yet 
to be considered. 

The actual extension of the centre line of the plate is 
{(T—T,)/E’t} where T,, is the initial tension in the plate, 
\’ = E/(1—o) for the plate material, and f¢ is the thickness 
of the plate; but the effective extension is further increased 
by the elastic deflection of the supporting members at 
the sides of the plate under the loading T—T,. Thus the 
total effective extension of ‘the centre line is 
{(T—T,)/E’t} +2{(T—T,)/kb}, woere k is the stiffness 
(load per unit width per unit defleci«.) of the supports 
at each edge of the plate to movement in **:* vlane of the 
plate*. The value of 7 is therefore to be determined from 
the equation 


*If the two stiffmesses are k, and k,, (1/k,)+(1/k,) should be 
substituted for (2/k). 
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T~f, ( 


The mean stress f, in the plate 


T _ El 
= 


73 (t/b)* (aby? 


so that {f,/E(t/b)} is a function of ab only. On the 
other hand the bending stresses which are functions of 
zb are also proportional to w and therefore to pb/T; and 
so, of course, is the deflection also. 


3. EXTENSION OF APPLICATION OF ROYAL AERONAUTICAL 
SOCIETY STRUCTURES DATA SHEETS IN THE 02.09 SERIES 
The data presented in R.Ae.S. Structures Data Sheets of 

the 02.09 series represent solutions of equation (4) for the 

case 7,=0 and (1/k)=0. If we define a nominal pressure 


where f,’ is the initial lateral stress, equation (4) reduces 
to 


T 


and the Data Sheets are once more directly applicable. 


° 
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Membrane Stresses and Defiections of a Long Plate 
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The procedure is to assume a value of p, and to read 
the value of f, from 02.09.01 or 02. Then, the values of 
f,, and k being specified, the value of p is found by multi- 
plying p, by the factor 


) (1+ 


The values of 6 and f,—f, read from the Data Sheets 
must also be multiplied by the same factor, because they 
are proportional to pb/T that is to (p/p,)(p,b/T), 
whereas the values read from the data sheets are propor- 
tional to (p,b/T). 

The procedure is demonstrated in Table II and com- 
parison with experimental measurements is described in 
Section 5 and Table L. 


4. APPLICATION TO PLATES INITIALLY SLACK 

The general method is applicable if the plate were 
initially under transverse compression, that is when f/f,’ 
is negative. However, if —f,’ > (x*/3)(t/b)* the plate 
would buckle and therefore a lower value of f,’ would 
represent initial slack in the plate. If the slack consisted 
in waviness of half wavelength A and half amplitude A, 
the corresponding strain would be (<A /2A) and the value 
of would be roughly — {3-6+(nzA 
where n is the number of half waves across the width b. 
Using this formula to assess the value of f,’ from measure- 
ments of (A/t) and n, the general method may be used 
to compute the behaviour of an initially slack plate. 

It is true that this application must be limited ‘to 
positive values of f, (final); but this restriction is of no 
consequence, because until f, > the initial waviness 
would persist. 

The application of the method to a slack plate is illus- 
trated in Fig. 1. 


5. EXPERIMENTAL CHECK 
The plate tested was 24 in. long, 4 in. wide and th in. 
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TABLE I 


Change of Total Lateral 

Lateral Stress (lb./in.*) Movement at 
Top Surface Bottom Surface Edges 
(in.) 


— 2,500 13,600 00049 


Pressure 
(1b./in.*) 


Note: Compressive stress is indicated by a negative sign. 
Since the measured change of average stress, 5,600 Ib./in.-, 
corresponds to a loading of 350 Ib./in. on the side supports, 
the value of Ib./in.2. Taking 
E’ =33x10® |b./in.?, 

1+ = 2°5 

The computed values of the stresses and deflection are 
listed in Table II. It will be seen that the computed value of 
f,—f,’ is in good agreement with the measured value, although 
the computed value of f,—f,’ is rather too low. 


thick. The ends of the plate were clamped to supports 
20 in. apart; these clamps held the edges of the plate 
firmly against rotation, but afforded only elastic restraint 
against movement of the edges in the plane of the plate. 
One support was adjusted relative to the other so as to 
apply an initial lateral tension of 1000 Ib./in.* to the plate. 
This lateral tension was measured by means of resistance 
wire strain gauges glued to both top and bottom surfaces 
of the plate at the centre of its length, and these gauges 
were subsequently used to measure the stresses developed 
when normal pressure was applied. The normal pressure 
was applied by dead weights distributed over a sheet of 
sponge rubber one inch thick laid on top of the plate. 

Measurements were also made of the lateral movement 
of the side supports by means of dial gauges in contact 
with each support just below the plate fixing. The results 
of the test are listed in Table I. 


TABLE Il 
(All stresses in lb./in.2: deflections in inches) 
p nom. | 3 3 40 


f, act. (read from 02.09.02) 4000 6900 11400 
f,’ 3000 5900 | 10400 


f, 6900 7306 9300 11400 
p nom. 


p act. 


. nom. (read from 02.09.02) 


act. 

(f,—f,) nom. (read from 
02.09.02) 

(f, —f,) act. 9950 
0°47 
21350 

1450 
20350 
450 
10400 


i 
3 
| 
1:36 2-15 | 2:94 3-08 | 4°47 | 6-04 
| 2:21 2-59 | 285 2:97 3-29 | 5-05 
| 
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Graduates’ and Students’ Section 


Annual Dance: Friday 20th November 


The dance this year will be held in the Library of the 
Society at No. 4 Hamilton Place, W.1. Dancing will be, 
as usual, to the music of Hugh McCamley. Double tickets 
cost only 14s. and include free drink and refreshments. 
No single tickets will be available this year. 

Get your tickets from your nearest Committee member 
or representative, or write to Mr. J. D. Duncan, 8 Carlisle 
Avenue, St. Albans, Herts. Cheques and Postal Orders 
should be made payable to “ Graduates’ and Students’ 
Section, Royal Aeronautical Society.” 


Visit to Hanover 

We would be interested to know if there would be 
support for a visit to the Hanover Fair (24th April-3rd 
May 1960) on the same lines as the Paris Air Show visit 
Part of the Fair is devoted to Air Transport and should 
prove of great interest. The cost is bound to be greater 
than the Paris visit and will depend on the numbers inter- 
ested, so write as soon as you can to the Hon. Visits 
Secretary: Mr. A. R. M. Pickering, “* Riverside,” Vicarage 
Walk, Bray, Berks. 


FUTURE LECTURES AND VISITS 


The Flight Testing of Rotorcraft 


Just to remind you of Sqdn. Ldr. Gellatly’s lecture on 
Wednesday 11th November; full details were given in the 
October JOURNAL. 


Atomic Energy Research Establishment, Harwell 

There may still be a few places left for this very popular 
visit on Saturday 21st November, when these words appear 
in print. Write quickly to the Visits Secretary, at the 
address above. 


Aerodynamics of Racing Cars 

Mr. F. A. Costin will be speaking to us on Wednesday 
2nd December on a very specialised subject of which he 
is an acknowledged expert. Shapes of record breaking 
Lotus, Vanwall, Maserati and Lister cars have been 
evolved from his careful work, and those who have met 
him will know that his fund of experience will ensure 
an interesting and entertaining evening. Frank Costin 
will be responsible for future shapes of projects from the 
Speedwell Company. 


College of Aeronautics, Cranwell 
A visit has been arranged for Saturday Sth December. 


Letter to the Editor 


Sir, 

I wish to comment on Mr. Wilby’s letter published in 
the Graduates’ and Students’ Section of the July 1959 
JOURNAL regarding the professional status of members of 
the Society. 

I wish to concur fully with Mr. Wilby’s sentiments. In 
my experience as an Aeronautical Engineer in both the 
U.K. and Canada, I have noticed this apparent lack of 
knowledge concerning the Royal Aeronautical Society and 
its qualifications and also the professional status of its 
members. 

In Britain (16 years ago) I found that preferential 
treatment in job status and salary was given to holders of 
B.Sc. degrees, who had no previous experience in the 
Aeronautical or any other Industry and were not members 
of the Society in any grade, over holders of a Grad. 
R.Ae.S. qualification. 


In Canada, when prospective employers question you 
regarding Professional Qualifications, an answer involving 
the Royal Aeronautical Society brings, in a surprising 
number of cases, the query * What does that mean?” or 
“Who are they?” The same confusion does not exist 
regarding the status of Mechanical, Electrical and Civil 
Engineers and membership in their respective Associations. 

I personally feel that today, professional recognition 
is becoming increasingly important for the Engineer as 
his work is becoming increasingly important to world 
progress, and it is essential that the public recognise and 
support this fact. 

I submit that the R.Ae.S. has fallen somewhat behind 
other institutions in recognising (or acting on this recogni- 
tion of) this fact of their raison d'étre. 

I have always been proud to be a member of the 
R.Ae.S. and continue to attach great importance to my 
association with this body. | certainly feel however, that 
more importance should in future be attached to the points 
raised in Mr. Wilby’s letter 

Yours faithfully, 
(signed) J. MACNAUGHTON, 
Grad.R.Ae.S., P.Eng., M.C.A.L. 


(A check on current advertising shows that only 
Blackburn and General Aircraft Ltd. specifically mention 
A.F.R.Ae.S. as a desirable qualification. The R.A.F. 
brochure on commissions in the Technical Branch does 
not mention it The Royal Navy has at present no 
vacancies for careers as i Technical Officer for anyone 
except Cadet entrants (aud Electrical Engineers) and we 
must therefore reserve judgment—Section Editor). 


Lunar Probes 


In his recent lecture to the Section about moon probes 
Mr. J. E. Allen outlined the flight paths and orbits which 
could be achieved and went on to examine the accuracy 
requirements for launching. He also discussed the design 
of the moon-probe vehicles and the launchings which 
have already been made by the U.S. and Russia. 

Mr. Allen showed that very small variations in the 
final velocity of a probe had a big influence on the journey 
time taken to reach the moon. He also gave figures for 
the accuracy needed if a probe were guided only during its 
launching phase. For a vehicle intended to hit the moon 
the acceptable velocity error was only 75 ft./sec—in a 
total speed of about 36,500 ft./sec—if the injection angle 
was exactly correct. If the correct injection speed were 
achieved, the acceptable error in angle was only 0-01°. 

For more complicated missions even greater accuracies 
were essential. If the probe was to pass round the moon 
and return to make a braking ellipse around the Earth, 
the initial accuracy required was 1 ft./sec. in velocity 
if the direction were correct; if the velocity were correct 
the direction error should be no more than 0-001°. Such 
accuracies were impossible with current guidance systems 
and mid-course or terminal corrections to the probe’s 
trajectory would be needed. 

The relative success of Russian and U.S. moon-probe 
launchings was mentioned in the discussion. Although 
the Russians may well have made abortive launchings, 
Mr. Allen thought it significant that recent Lunik shots 
had been timed to coincide with Kruschev’s visit to the 
U.S. and the second anniversary of the Sputnik I launch. 
This indicated that the Russians could make successful 
launchings on planned dates, which had not always been 
the case with U.S. efforts. 

Mr. Allen’s complete lecture will be published in a 
future issue of the JOURNAL.—J. R. C. 
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Reviews 


SOUNDING ROCKETS. Homer E. Newell, Jr. McGraw- 
Hill. 334 pp. Illustrated. 97s. 


In Dr. Newell's words “ The essential feature of a 
sounding rocket is its ability to carry scientific instruments 
into the upper atmosphere.” In his book he has gathered 
together descriptions of a very wide range of rockets with 
this ability and has outlined both the nature and the pur- 
pose of the instruments which are carried. The fact that 
Dr. Newell has nineteen collaborators, as authors of 
various chapters, stems from the fact that the book is 
based on papers presented in November, 1956 to a Sym- 
posium on High Altitude Sounding Rockets organised by 
the American Rocket Society. In the course of compila- 
tion the material has been kept up to date so that it, in 
effect, gives the “state of the art” at the start of the 
International Geophysical Year in July 1957. 

To the separate descriptions of some fifteen basic 
sounding rockets and their variants, Dr. Newell has added 
two introductory chapters and two concluding chapters. 
The first introductory chapter outlines the essential 
features of sounding rockets and goes on to give the basic 
mathematics of rocket propulsion, the thermodynamics of 
rocket jets and the theory of the convergent-divergent 
nozzle. Why these particular topics are given so much 
weight is puzzling for they do not seem any more peculiar 
to sounding rockets than other important aspects of rocket 
flight, such as aerodynamics, which are not discussed. 
The second introductory chapter outlines the application 
of rockets to geophysical and solar research. In six pages 
the physical properties of the atmosphere (density, pres- 
sure, temperature, composition, clouds, winds and turbu- 
lence, the ionosphere, magnetic field, aurora, airglow, 
cosmic rays, meteors) are summarised. In the following 
seven pages the principal methods and the results obtained 
in measuring the above properties by rocket techniques 
are given. The high degree of condensation required to 
cover such a wide field in thirteen pages results in an 
account which should be of interest to those on the fringe 
of this field but is unlikely to appeal to either the novice 
or the expert. 

The bulk of the book consists of fourteen chapters, 
each of which is devoted to a particular rocket or rocket 
family. In addition to American rockets (Aerobee, Viking, 
and so on) the sounding rockets developed by the British 
(Skylark), the French (Veronique and Monica) and the 
Japanese (Kappa and Sigma) are described in detail. For 
each rocket the historical background is given; the pro- 
pulsion, launching and constructional features described 
and the performance characteristics stated. In some in- 
stances rocket costs are quoted. The techniques involved 
in launching from balloons or aircraft as a means of 
achieving greater altitude are described. 

In his penultimate chapter Dr. Newell describes briefly 
the rocket firing ranges available in the U.S.A., Canada, 
Australia, French N. Africa and Japan. In his final 
chapter he discusses that complement to the sounding 
rocket—the artificial satellite. In this he derives the basic 
orbits and the relationships which have to be established 
in order to achieve orbital motion. The Vanguard launch- 
ing vehicle is very briefly described. 

I would have welcomed an extra chapter which 
attempted to bring together and compare the main charac- 
teristics of the rockets so that the reader was left not with 


a confused impression of many rockets, but an integrated 
picture of the trends in sounding rocket development. 

The main value of the book undoubtedly lies in its 
being a detailed, in places too detailed, history of sounding 
rockets up to mid-1957.—J. J. GaIT. 


GUIDED MISSILE ENGINEERING. Allen E. Puckett and 
Simon Ramo (Editors). McGraw-Hill, London, 1959. 497 pp. 
Illustrated. 77s. 6d. 


Despite its title, this book is not really concerned with 
engineering: it might more suitably have been called “The 
fundamental theory of guided missiles.” In seventeen 
independent chapters, by separate authors, the theory 
underlying nearly all aspects of guided missile design and 
operation is developed ab initio up to a point where a 
reader previously unfamiliar with the theory can see how 
it applies to actual missiles. The editors chose most of 
the authors from their own organisations, the Hughes 
Aircraft Company and the Ramo-Wooldridge Corporation, 
but there need be no regrets over this touch of clannishness, 
for the book is an excellent one. 

The first section is a deft outline of aerodynamics by 
A. E. Puckett, and this is followed by chapters on aero- 
dynamic stability and control, rockets, and ram-jets. Then 
come surveys of information theory and missile guidance, 
navigation, random phenomena, servos and feedback 
control, and gyroscopy. Next there are three chapters on 
electronics—entitled Radio Propagation, Microwaves, and 
Radar. Thus far, the book is consistently good: the 
general plan seems to be that each author should start at 
the most elementary level and develop the essentials of the 
theory, sometimes to quite an advanced stage, without 
digressing into matters of detail. The authors succeed in 
this aim and are often refreshingly direct in cutting through 
the tangle of jargon which shrouds their subject, most 
notably perhaps in F. P. Adler’s presentation of infor- 
mation theory and its application to missile guidance, 
topics notoriously difficult to explain. The overlapping, 
omissions and sudden changes of style which often mar 
multiple-author books are avoided. The last four chapters 
on analogue computers, digital computers, trajectory 
analysis and flight simulators, and the integration of design, 
are however less satisfactory, chiefly because of overlapping 
and a drift into generalities. 

The editors, like the compilers of a verse anthology, 
can scarcely hope to escape criticism over their choice of 
subjects, on which opinions are bound to vary. To the 
reviewer it is astonishing, especially in view of the book's 
title, that there is no chapter on missile structural design: 
the intercontinental ballistic missile, which is, at least to 
the general public, the guided missile par excellence, 
depends more upon efficient structural design than upon 
any other factor. The balance of the book is questionable 
in a few other respects, too. For example, it is rather 
surprising to find a 40-page chapter, backed up by an 
Appendix, on such a specialised subject as instrumental 
gyroscopy, when the other chapters are much wider in 
scope. 

With these few reservations, the book can be strongly 
recommended as an introduction to guided missile theory; 
and it conforms to the high standard of production which 
one expects of McGraw-Hill books.—b. G. KING-HELE. 
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BRITISH CIVIL AIRCRAFT 1919-59. Vol. 1. A. J. Jackson. 
Putnam, London, 1959. 571 pp. Illustrated. 63s 

In his preface, the author of this latest book in Putnam’s 
historical series points out that the progress of civil aviation 
in Britain falls into two well-defined 20-year cycles, the 
first distinguished by relatively unfettered operation of 
privately-owned light aeroplanes, the second dominated by 
fiercely competitive transport fleets and strict legislation. 
Somewhere about halfway through the first cycle, the 
author was inspired to begin collecting photographs and 
data of aeroplanes bearing British civil registration marks. 
G-EBZZ had just been followed by G-AAAA and both had 
been in the news, the former with a solo flight to West 
Africa and the latter with a new altitude record. Even at 
that date the possibility of compiling anything like a 
complete record must have seemed remote, but by fortu- 
nate coincidence, he, living in Essex, made contact with 
another of like mind living in the North of England. Their 
collections, being regional in origin, were found to be 
complementary and so began the fruitful partnership 
between A. J. Jackson and the late E. J. Riding, to whom 
this book is dedicated. 

Their complete records, compiled for personal satis- 
faction, form an unparalleled archive of one aspect of 
civil aviation over a period of forty years, but are obviously 
far too voluminous and detailed for publication at length. 
Mr. Jackson has had a difficult task in epitomising this 
record in a book of reasonable compass even though his 
publishers have allowed him two pint pots into which to 
distil his imperial gallon. No room could be spared for 
excursions into anecdote or design philosophy, and some 
readers may regret this, but Mr. Jackson is to be congratu- 
lated on providing a very readable text containing all the 
essential facts, dimensions, fates and dates on his subject. 
This volume ranges from A.B.C. to Dudley Watt, the 
remainder being reserved for Volume 2. 

The production and quality of the photographs is fully 
up to Putnam’s highest standards, and what a store of 
hidden treasure they reveal! For the first time this 
reviewer sees in print the Airspeed Ferry with all three 
engines inverted, the Avro Cabin Cadet, the Bristol 83E 
with geared Titan, the Parnall-built Airdisco-engined C.11 
gyroplane, the Avro 538 sporting variant of the Spider 
scout, the Beardmore W.B.X., the Dart Flittermouse, the 
only civil Bristol Scout C and the D.H. Hound with its 
original direct-drive Lion. The other 500 photographs, 
although less rare, are all excellently chosen and in many 
cases hitherto unpublished. There are over 100 three- 
view general arrangement drawings, which exhibit Mr. 
L. E. Bradford’s customary skill and accuracy, although 
unavoidably constrained to a small scale in certain cases. 

The main body of the book deals with 210 types and 
variants of production civil aircraft, giving each main type 
or family group a concise and informative account of 
development and operation averaging some 1500 words. 
This section is followed by four Appendices, of which the 
first briefly describes and illustrates a further 141 civil 
aeroplanes, mainly prototypes; the second deals similarly 
with 45 civil-registered military aeroplanes; the third 
describes the mere dozen civil-registered types of which no 
illustrations have so far come to light; while the last gives 
in tabular form the registration particulars and individual 
histories of aircraft described in the main text, so far as 
they are on record. 

No book of this size and scope can entirely avoid errors, 
but very few are apparent and those are of minor impor- 
tance. For instance, the Genet-Avian G-AABR (a most 
interesting close-cowled variant) is captioned as “wooden” 
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although stringers visible through the fuselage fabric are 
typical of metal construction; surely, too, the Cirus III 
Avian AACV and Tommy Rose’s Hermes Avian AAHJ 
had metal fuselages and antedated AATL which is quoted 
as being the first metal Avian? There is some confusion 
over the type number of the Avian IVM—on page 91 it is 
given as Avro 616 but elsewhere as Avro 594. The Beard- 
more W.B.X is identified in the text as G-EAQI whereas 
the photo clearly shows it as EAQJ; on page 133 Mr. J. E. 
Carberry is wrongly styled “ Sir John,” which may annoy 
him, and there are a few mis-spelled names (e.g. Wooley 
Dod), which should have been avoided. These will 
not mislead any reader but ought to be eliminated from 
the next edition, in which, too, the opportunity should be 
taken to bring last-minute data up to date, e.g. Britannia 
AOVT should be added to the B.O.A.C. fleet list in 
Appendix D. An author’s difficulties in gleaning the last 
grain of information are not confined to current events, 
for, after this volume had gone to press, news came from 
Australia indicating almost certainly that Harry Butler’s 
preserved Bristol monoplane VH-UQI was bought by him 
from the Disposal Board at Waddon as C-5001 at an earlier 
date than G-EAER was registered; consequently the sur- 
mise given in good faith on pages 186 and 534 must now be 
discounted. 

An interesting sidelight on the small productive capacity 
of the British aircraft industry in the late twenties is 
implied on page 164, where it is mentioned that Bluebird 
IV production was subcontracted to Saunders-Roe because 
Blackburn’s had just received a military order (presumably 
for Ripons). What Mr. Jackson does not say is that Saro 
themselves were so overwhelmed by the Bluebird com- 
mitment that they in turn sub-contracted the wings to 
Boulton and Paul at Norwich. Mr. Jackson writes as a 
lifelong flying club enthusiast and is an active and talented 
amateur holder of a P.P.L., who has flown a representative 
sample of the vintage light aeroplanes he describes. He 
seems, moreover, to have captured much of the zest which 
was so characteristic of Eddie Riding’s own reporting. 

This is a book in the high tradition set by J. M. Bruce’s 
“British Aeroplanes 1914-18,” and a thoroughly worthy 
companion to Thetford’s two precursors in the same 
format. One has only one fear, that Volume 2 may not be 
out in time for Christmas—but that is no reason for not 
ordering both volumes now.—SENRAB. 


NINE LIVES. Group Captain A.C. Deere, D.S.O., O.B.E., D.F.C. 
Hodder and Stoughton, London, 1959. 262 pp. Illustrated, 15s. 


Alan Deere, the most famous of the New Zealand 
fighter pilots, though a forceful character, as his autobio- 
graphy makes clear, is a very modest man. His personal 
courage in conditions of desperate war against the enemies 
Luftwaffe and Fatigue was a great inspiration to all 
who worked, flew and fought with him. His story of a 
very unpleasant yet thrilling experience as a pilot in the 
Battle of Britain is probably the best yet published. As 
a word picture it has everything to commend it. The book 
is written in the first person and the author takes the 
reader with him, as a guide shows a visitor a military 
museum. He shows from personal experience how the 
battle ebbed and flowed, as he saw it rather than as it 
affected him. Of course, he was tremendously affected, 
as several incidents make clear, but his personal enthusiasm 
for the job of making war on a vicious enemy submerges 
his personal sense to a remarkable degree. This book is 
well worth reading.—A.S.C.L. 
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Additions to the Library 


Air Aces of the 1914-1918 War. Bruce Robertson 
(Editor). WHarleyford Publications, Letchworth. 1959. 
211 pp. 45s. To be reviewed. 

Asbestos: its Industrial Applications. D. V. Rosato. 
Reinhold (C. and H.). 1959. 214 pp. Illustrated. 
46s. A small book on uses of asbestos in industry. 
Chapters include insulation, textiles and plastics. A 
bibliography is appended. 

British Civil Aircraft 1919-59, Volume I. A. J. Jackson. 
Putnam, London. 1959. 571 pp. Illustrated. 63s. 
Reviewed on p. 677. 

Construction of Research Films, The. D. H. Densham. 
Pergamon. 1959. 104 pp. 30s. No. 33 in the Agardo- 
graph series, this has been compiled by the Adviser on 
Films at the Ministry of Supply. Thirteen chapters deal 
with all aspects and a useful glossary is appended. 

Famous Airships of the World. J. A. Sinclair. Muller, 
London. 1959. 144 pp. Illustrated. 9s. 6d. There 
are 18 other titles in this series—airports, air routes, 
battleships and so on through the alphabet to water- 
ways. The late Captain Sinclair was one of the few 
people still left with an interest in airships. Such a 
small book cannot be exhaustive and most of it deals 
with the bigger craft of the later years. 

Handbook of Supersonic Aerodynamics. Volume 3, 
Section 7. Three Dimensional Aerofoils. U.S. Bureau 
of Ordnance and Johns Hopkins University. U.S.G.P.O. 
1957. Sectional pagination. $1.50. 

Numerical Control of Machine Tools, The. R. C. Brewer. 
Engineer’s Digest. 1959. 69 pp. Illustrated. 6s. A 
second edition of Engineer's Digest Survey No. 5, 
originally published in that journal in September 1958 
and January 1959. A long list of acknowledgments 
indicates co-operation from many firms. 

Outer Space Propulsion by Nuclear Energy. U.S. Con- 
gress. U.S.G.P.0. (H.M.S.O.). 1958. 232 pp. 5s. Sd. 
Minutes of hearings before subcommittees of the Joint 
Committee on Atomic Energy. This Committee contri- 
buted to the Nautilus project and now has Project 
Rover underway. “It was the consensus among wit- 
nesses who appeared before the Committee that of all 
propulsion systems currently under consideration, 
nuclear propulsion appears to offer the greatest advan- 
tage for carrying large payloads great distances, such 
as interplanetary travel.” 

Proceedings of the Society for Experimental Stress 
Analysis. Volume XVI, No. 2. C. V. Mahlmann and 
W. M. Murray (Editors). S.E.S.A. 1959. 165 pp. 


Illustrated. No price. 17 papers on a range of 
aspects of stress analysis, including photoelasticity, 
triaxial stress testing, acoustical testing, impact loading 
and the use of digital computers. 

Raumfahrt—technische Uberwindung des Krieges. Eugen 
Sanger. Rowohlt, Hamburg. 1958. 142 pp. D.M. 1.90. 
To be reviewed. 

Rockets and Satellites Work Like This. John W. R. 
Taylor. Phoenix House London. Roy Publishers, 
New York. 72 pp. Illustrated. 9s. 6d. One of the 
“Science Works Like This” series, intended for boys 
but suitable for anyone who wants a simple, logical 
presentation of the history and development of rockets 
and satellites. It is well co-ordinated and illustrated, 
and could be used as a reference book for a short, 
clear description of nine basic categories of guided 
weapon and eleven types of guidance system. 

S.A.E. Transactions 1959. Society of Automotive Engin- 
eers. New York, S.A.E. 1959. 672 pp. Figs. $8 
to libraries. $11 to non-members. Contains some of 
the more notable papers read in 1958 including Van 
Every on design for very high speed flight, Stepniewski 
on wings and vertical thrust in VTOL, Ekvall on strength 
and limitations of materials for airframes (Titanium) 
and Stasinger on inertial navigation systems. 

Seven Skies, The: a Study of B.O.A.C, and its Forerunners 
Since 1919. John Pudney. Putnam, London. 1959. 
320 pp. Illustrated. 30s. To he reviewed. 

Sixth Anglo-American Aeronautical Conference. Folke- 
stone, 1957. Royal Aeronautical Society, London. 1959. 
524 pp. 6 guineas. The complete record of the Technical 
Sessions, including the 16 papers, all formal discussions 
and the list of delegates. The subjects covered include: 
N.A.C.A. Research on V.T.O.L. and S.T.O.L., Nuclear 
Propulsion for Aircraft, Control of Supersonic Pro- 
pulsion Systems, Electrical Control of Power Plants, 
Kinetic Heating, Titanium, Non-metallic Materials at 
High Temperatures, Blast Loading of Aircraft Struc- 
tures, Safety in Relation to Structural Damage, Aircraft 
Instruments, Aeroelasticity, Boundary Layer Control 
Experimental Hyperballistics, Civil Turbine Engines 
and The Human Pilot. 

Two Residents of Bexley, Part I—Sir Hiram Maxim. 
Patricia H. Hazael. Bexley. 1959. 55 pp. Illustrated. 
No price. A teacher's thesis which, with its local 
interest and modern approach is more readable and 
interesting than the conventional works. 
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AERODYNAMICS 
BOUNDARY LAYER see SCIENCE—-GENERAL 


COMPRESSIBLE FLOW 


An experimental investigation of the interference between bodies 
of revolution at transonic speeds with special reference to the 
sonic and supersonic area rules. G. Drougge. Appendix: 
Application of the transonic principle of equivalence to a pair 
of bodies of revolution. S. B. Berndt and L. O. A. Hilding. 
F.F.A. Report 83. 1959. 

Certain aspects of the sonic and supersonic area rules have 
been investigated experimentally, mainly by interference tests 
on a pair of slender bodies of revolution at transonic and low 
supersonic speeds. The investigations make it possible to 
test the transonic similarity rules for bodies of revolution, and 
a convincing confirmation is obtained.—(1.2.2 x 1.2.3). 


Fium Dynamics see also STABILITY AND CONTROL 


The effect of density on jet flow at subsonic speeds. A. H. 
Craven. C.0.A. Report 120. July 1959. 


On the assumption that the velocity and density distributions 
across a jet of one gas issuing into a stream of a second gas 
are of exponential form, the momentum integral approach 
has been used to find the variation of centre-line velocity down- 
stream of the potential core. The results for the particular 
case of a jet issuing into a gas at rest are equivalent to those 
found theoretically by Blottner and experimentally by Keagy 
and Weller. The results for the jet issuing into a moving 
stream agree with the author's limited experiments.—(1.4). 


A_ fluid-dynamic mechanism of meteorite pitting. D. T 
Williams. Smith: Cont. to Astro. Vol. 3, No. 6. 1959. 
Flow lines on meteorites were photographed to determine the 
nature of air flow inside the “thumb pits” in the surface. 
(1.4.3 x 8.1). 


Loaps see WINGS AND AEROFOILS 


STABILITY AND CONTROL see also TESTING AND INSTRUMENTS 


Effect of body incidence on two afterbodies with a rearward 
facing jet. A. H. Craven. C.0.A. Note 92. Dec. 1958. 


= = 
| 
lee 
: 
af 
4 
: 
2h 
2 


THE LIBRARY—REPORTS 


The results of an experimental investigation into the effect 
of a jet issuing from two afterbodies (one a right cylinder and 
the second conical) at incidence to a uniform subsonic flow 
are given. The tests were performed at a Reynolds number 
of 0:3x10® based on a body diameter and maximum tunnel 
velocity.—(1.8.0.2 x 1.4) 


The influence of non-stationary stability derivatives on the 
snaking motion at high subsonic speed. J. Yff. N.L.L. Report 
TR.F.208. Nov. 1957. 

Some lateral stability calculations with two degrees of freedom 
for flight at high subsonic speed have been performed to 
study the influence of the non-stationary stability derivatives 
on the damping of the snaking motion.—(1.8.1.1) 


WINGS AND AEROFOILS 


Efforts d'un fluide incompressible sur un profil oscillant autour 
d'une incidence non nulle. G. Couchet. O.N.E.R.A. Pub. 96. 
1959. (In French).—(1.10.1.1 X 1.6.3 X 2). 


Experimental investigation of transonic flow past two-dimen- 
sional biconvex circular-arc airfoils at small angles of attack 
R. Kawamura and K. Karashima. Aero. Res. Inst, Tokyo 
Report 342. April 1959 

The results of an experimental investigation of transonic flows 
past biconvex circular-arc aerofoils with  thickness-chord 
ratios of 10°6, 8-08 and 5°44 per cent are presented. Aero- 
dynamic characteristics of the aerofoils are analysed from the 
surface pressure distributions obtained by interferometry 
(1.10.2.1) 


TESTING AND INSTRUMENTS see also FLIGHT TESTING 


Measurement of aerodynamic damping and stiffness derivatives 


in’ free oscillation with automatically recycled feedback 
excitation. K. J. Orlik-Riickemann. N.R.C. Report LR-246. 
June 1959.—(1.12.6 x 2 x 1.8.0.2) 


Interpretation and applicability of results of wind-tunnel flutter 
and control surface buzz investigations, A. G. Rainey. AGARD 
Report 219. Oct. 1958 

Several aspects of the réle of the wind tunnel in the overall 
picture of flutter investigations, model simulation, tunnel 
suitability and application to flight are considered. The types 
of flutter amenable to wind tunnel study, the degree of model 
simulation required, the value of the small scale aeroelastic 
component model and, the application of wind tunnel results 
in the determination of critical flutter ranges or margins of 
safety are described.—(1.12.1 x 2). 


A summary of design techniques for axisymmetric hypersonic 
wind tunnels. Y-N Yu. AGARDograph 35. Nov. 1958 

A family of perfect fluid co-ordinates for axisymmetric effusers 
with the exit Mach number range of 5 to 20 has been com- 
puted by the method of characteristics. A conical source flow 
front and a centre line Mach number distribution are used as 
initial conditions. Design techniques and other related problems 
of interest in the development of hypersonic facilities are 
summarised.—(1.12.1.3 « 1.12.6.1). 


AEROELASTICITY 


See AERODYNAMICS—WINGS AND AEROFOILS 
TESTING AND INSTRUMENTS 


AIRCRAFT 
See DESIGN AND CONSTRUCTION 
DESIGN AND CONSTRUCTION 


Tire-to-surface friction-coefficient measurements with a C-123B 
airplane on various runway surfaces. R. H. Sawyer and J.J. 
Kolnick. N.A.S.A. Report 20. June 1959. 

An investigation was conducted to obtain information on the 
tyre-to-surface friction coefficients available in aircraft braking 
during the landing run. The tests were made with a C-123B 
aeroplane on both wet and dry concrete and bituminous pave- 
ments and on snow covered and ice surfaces at speeds from 
12 to 115 knots. Measurements were made of the maximum 
(incipient skidding) friction coefficient, the full-skidding 
(locked wheel) friction coefficient, and the wheel slip ratio 
during braking.—(4.2.2.3 x 5.3 x 3.6). 


United States Air Force development of aircraft arresting gear. 
V.V.Vary. AGARD Report 227. Oct. 1958. 

For many years Navy aeroplanes have required arresting gear 
and catapults for carrier operation. This has meant a consider- 
able weight penalty to the aeroplane structure, which has had 
to be designed to withstand horizontal arresting accelerations 
up to 6g and also large vertical accelerations due to the rela- 
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tively high vertical speed of contact with the deck. A 
relatively low stalling speed—to keep within required launching 
and arresting energies—has been another penalty. Land based 
Air Force aeroplanes have, until recently, gone in for longer 
runways and, consequently, can operate at greater loads and 
speeds. This makes the design of arresting gear for Air Force 
aeroplanes much more difficult than for the short-stroke, 
high-g Navy aeroplane. Various devices by which engage- 
ment can be effected and the aeroplane successfully arrested 
are discussed.—(4.2.3) 


AIRCRAFT OPERATION 


See also DESIGN AND CONSTRUCTION 

FLIGHT TESTING 

AIRPORTS 
Origin and prevention of crash fires in turbojet aircraft. 1. 1. 
Pinkel et al. N.A.C.A. Report 1363. (Supersedes N.A.C.A. 
T.N. 3973). 1958 
The tendency for the jet engine rotor to continue to rotate 
after crash presents the probability that crash-spilled com- 
bustibles suspended in the air or puddled on the ground at 
the engine inlet may be sucked into the engine. Studies showed 
that combustibles drawn into the engine in this way ignite 
explosively within the engine Flames that could set the 
major fire appear at the engine tailpipe and often at the 
inlet when this explosive ignition occurs. This ignition may 
occur on the hot metal of the engine interior even after the 
combustor flame is extinguished and the engine rotor is coast- 
ing to rest.—{5.3 x 27.4) 
An evaluation of the factors which influence the selection of 
landing approach speeds. F. J. Drinkwater et al. AGARD 
Report 230. Oct. 1958 
The aim was to increase understanding of factors influencing 
pilots’ choice of approach speeds by an evaluation of landing 
techniques of a wide variety of aeroplanes, using the constant 
speed, constant angle, type of approach.—(5 x 13). 


AIRPORTS 


Measurements and power spectra of runway roughness of an 
airport. (Airport 1). AGARD R.M. 1. March 1959. 
Measurements of the surface roughness obtained by a profile 
survey method are presented in the form of tables giving the 
survey data and in the form of power spectra.—(6.1), 


Measurements and power spectra of runway roughness of an 
airport. (Airport 2). AGARD R.M. 2. March 1959. 
Measurements of the surface roughness obtained by a profile 
survey method are presented in the form of tables giving the 
survey data and in the form of power spectra.—(6.1). 


EXTRA-ATMOSPHERIC TECHNOLOGY 


See also AERODYNAMICS—FLUID DYNAMICS 
MATHEMATICS 
MISSILES 


Minimum propulsion for soft moon landing of instruments. 
D.S. Carton. C.0.A. Note 94. July 1959.—(8 x 27:3). 


FLIGHT TESTING 


See also AIRCRAFT OPERATION 


Review of aircraft altitude errors due to static-pressure source 
and description of nose-boom installations for aerodynamic 
compensation of error. W. Gracey and V. S. Ritchie. N.A.S.A. 
Memo 5-10-59L. TJ.L. 6491. June 1959. 

Aircraft altitude errors due to static pressure source are 
briefly reviewed. Wind tunnel test information concerning 
nose-boom pressure installations tor aerodynamically com- 
pensating static-pressure errors due to position and angle of 
attack at subsonic and low supersonic speeds are also presented. 

{13 x 1.12.2) 


Several methods for aerodynamic reduction of static-pressure 
sensing errors for aircraft at subsonic, near-sonic, and low 
supersonic speeds. V. S. Ritchie. N.A.S.A. Report 18. June 
1959. 

Tests were conducted in transonic wind tunnels to investigate 
and verify experimentally methods for aerodynamically reduc- 
ing errors due to sensor position, bow-wave passage, and 
angle of attack.—{13 x 1.12.2). 


The performance testing of military aircraft-carrier aircraft. 
R. P. Dickinson. AGARD Report 232. Oct. 1958. 

The operation of carrier-borne aircraft and the ways in which 
it differs from the operation of other military aircraft from the 
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viewpoint of performance assessment are described. The factors 
affecting the determination of recommended approach speed 
and of wind speed limitations for catapult take-off are discussed 
and an indication is given of the extensive instrumentation 
needed during trials.—(13.3 x 5). 


FUELS AND LUBRICANTS 


Aviation fuel problems at high altitudes and high aircraft 
speeds. E. M. Goodger. C.0.A. Report 119. Jan. 1959. 

Much useful data has appeared concerning the problems 
incurred by continued increases in operational altitudes and 
aircraft speeds. A representative amount of these data is 
ge in a form suitable both for general information and 
or project design reference.—(14), 


INSTRUMENTS AND EQUIPMENT 


R.P.M. indicator and overspeed trip. R. J. Harron. N.RC. 
Report MI-819. May 1959. 

The function of this dual-purpose instrument is (/) a coarse 
r.p.m. indicator and (ij) an adjustable overspeed trip relay 
having two overlapping ranges of 11,000 to 16,000 r.p.m. and 
15,000 to 20,000 r.p.m.—(18.1). 


MATERIALS 
See STRUCTURES—-THEORY AND ANALYSIS 


MATHEMATICS 
See also MIsSSsILes 


Periodic orbits of a planetoid passing close to two gravitating 
gy R. R. Newton. Smith: Cont. to Astro, Vol. 3, No. 7, 
A study of periodic orbits in the restricted three-body problem 
in which the small mass passes close to each of the large 
masses is presented. The orbits are classified according to the 
types of motion involved and methods are given for estimating 
initial conditions for various classes of orbit. Sample orbits 
determined by a digital computer are presented.—(22.3 x 8). 


Actes du coiloque de calcul numérique. Perigueux 1957. Pubs 
Se et Tech N.T. 80. (In French.) 1959.—422.1). 


MECHANICAL ENGINEERING 


On helical-spring suspension. L. G. Lidin. Inst. of Machine 
Elements. Chalmers University of Tech. Report 6. 1959. 

The problem of determining the natural frequencies for a device 
suspended by helical springs is treated. The method is applic- 
able also to other spring types, provided their elastic properties 
are known. There are two different types of springs that are 
used for this purpose, namely springs with hinged ends (hooks) 
and compression springs with ends clamped. An account of 
the elastic properties of both types is given. Then for each 
type the general frequency determinant for a spring-suspended 
system is derived from the movement equation. The significance 
of the natural frequencies for forced vibrations is demonstrated. 
General schedules are given for obtaining the determinants for 
any case of helical-spring suspension. Experiments have been 
made of both types of springs.—(23). 


MISSILES 


1958 NASA/USAF space probes (Able-1). Final Report Vol. 1. 
Summary. Space Technology Labs. N.A.S.A. Memo 5-25-59W. 
TIL 6490. June 1959. 

The three NASA/USAF lunar probes of 17 Aug., 13 Oct. and 
8 Nov. 1958 are described. Details of the programme, the 
vehicles, the payloads, the firings, the tracking, and the results 
are presented. Principal result was the first experimental veri- 
fication of a confined radiation zone of the type postulated by 
Van Allen and others.—(25.3 x 27.3 x 8). 


Lieu des racines dune équation algébriaue dépendant d'un 
paramétre: application a la stabilité et au guidage des fusees. 
F. Salles. Pubs Sc et Tech. 351. 1959. (in French).—(25.1 x 


22.1). 
POWER PLANTS 


See also AIRCRAFT OPERATION 
ExTRA-ATMOSPHERIC TECHNOLOGY 
MISSILES 


The effect of combustor outlet temperature profiles on ramiet 
thrust. M. R. Williams and S. W. Greenwood. C.0.A. Note 59. 
July 1957. 

The desirability of obtaining an even distribution of temperature 
at the outlet of a ram-jet combustor is indicated theoretically. 
Results of experiments on a particular combustion system over 
a limited range of conditions using a series of mixers are 
presented.—{ 27.4). 
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Theoretical performance of liquid hydrogen with liauid oxygen 
as a rocket propellant. S. Gordon and B. J. McBride. N.A.S.A. 
Memo 5§-21-59E. TIL 6492. June 1959. 

Theoretical rocket performance for both equilibrium and 
frozen composition during expansion was calculated at four 
chamber pressures (60, 150, 300 and 600 Ib./in.* abs.) and 
a wide range of pressure ratios (1 to 4000) and oxidant-fuel 
ratios (1-190 to 39°683). Data are given to estimate performance 
parameters at chamber pressures other than those for which 
data are tabulated.—(27.3 x 34.1.1). 


SCIENCE—GENERAL 


La resistance thermique des joints. J. J. Bernard. AGARD 
Report 212. Oct. 1958. (In French.)}—(32.2.2 x 34.3.3). 


Recherches expérimentales sur les ondes de choc produites par 
les décharges condensées dans les gaz rares. M. Cloupeau. Pubs 
Se et Tech. 352. 1959. (In French.)}—32.2.1 X 1.2.3.2). 


STRUCTURES 
THEORY AND ANALYSIS 


An improved first-approximation theory for thin shells. J. L. 
Sanders. N.A.S.A. Report 24. June 1959. 

A theory is derived in which all strains vanish for any rigid- 
body motion in contrast to the results of Love's theory. 
Expressions for the stress resultants and couples which satisfy 
the homogeneous equilibrium equations are given in terms of 
three stress functions. The special forms of the equations of the 
new theory in the case of a circular cylinder are given.— 
(33.2.4.11). 


Sandwich structure for high temperature vehicles. H. R. 
Ashley. AGARD Report 216. Oct. 1958. 

A review of high temperature sandwich structure from the view- 
points of theoretical structural efficiency and of practical design 
and manufacturing is presented. The physical properties of 
five materials (R.R.58 light alloy, F.V.520 stainless steel, 6 Al 
4 Va titanium alloy, Nimonic 90 anc “5 titanium molybdenum) 
are summarised and, in conjunction with core data, used to 
present a series of graphs showing the efficiency of both, honey- 
comb and corrugated sandwich at high temperatures for 
specific structural index values. Thermal stress factors for these 
materials at various temperatures are also presented. The prac- 
tical design and production aspects of sandwich construction is 
discussed.—(33.2.4.14.9 x 21.2.3). 


Transient thermal stresses in the elastic slab. C. H. J. Johnson. 
ARL Note ME 232. March 1959. 

As part of an investigation of thermal stresses in turbine blades 
a solution has been obtained to the problem of determining the 
transient thermal stress distribution in an infinite elastic slab 
due to the sudden simultaneous application of a temperature 
T: to one face and a temperature T2 to the other face, the initial 
conditions being zero and with zero surface tractions. The 
solution is given in series form, the successive terms of which, 
in effect, give the stress at a particular depth in a semi-infinite 
solid with a constant temperature boundary condition. This 
result is analogous to the solution of problems in transmission 
line theory in terms of multiple reflected waves. The solution 
in this form converges quite rapidly except for large values of 
the time.—(33.2.4.5.9). 

TESTING 

A calibrator for strain gauge bridges. I. G. Scott. ARL Note 


SM 251. March 1959. 
A very simple calibrator for strain gauge bridges is described. 


(33.3.1). 
THERMODYNAMICS 


See also POWER PLANTS 
ScIENCE—~GENERAL 


The freezing and melting times of slabs and cylinders. D. C. 
Baxter. N.R.C. Report MK-1. Feb. 1959. 

General results for the time required to freeze or melt slabs 
or circular cylinders of a pure substance initially at the fusion 
temperature and exposed to a constant ambient temperature 
have been obtained on the PACE analogue computer.—(34.3.4). 


Heat conduction in solids at hypersonic speeds. L. Broglio. 
AGARD Report 209. Oct. 1958. 

Part I is an extension of a general approach to the heat conduc- 
tion problems presented in a previous paper to a body whose 
thermal properties are dependent on temperature. Part II is 
an application of Part I to the heat conduction problems in 
hypersonic flight. Part III is a numerical application of Part | 
reting II to a hollow steel hemisphere at hypersonic speed.— 
(34.3.3). 
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APPOINTMENTS 


This section of THE JouRNAL is available for advertisements of appointments in the Industry, 
the Ministries, Research Establishments, Universities and Colleges. 


Remittances—-Cheques and posta! orders should be made payable to the 
Journal of the Royal Aeronautical Society, One Arundel Street, Strand, 
London, W.C.2. 

The Society reserves the right to decline any copy or advertisement at its 
discretion and accepts no responsibility for delay in publication or for 
clerical or printer's errors, although every care is taken to avoid mistakes. 


Press Day—20th of the month preceding publication. 
Rates—8/- a line. Each paragraph is charged separately and name and address 
must be counted Semi-displayed setting £4 Os. Od. per column inch 
Box Numbers—1/- extra Replies should be addressed to: Box 000, care of 
THE Journat of the Royal Aeronautical Society, One Arundel! Street, 
Surand, London W.C.2 


UNIVERSITY OF BRISTOL 


1) Applications are invited for a LECTURESHIP (Grade II) in 
AERONAUTICAL ENGINEERING. Applicants should have 


Bla ckb url) a good Honours Degree and some relevant research or teaching 
experience, preferably in the field of High Speed Flow. 


offer The salary will be in the range £900—£1,300 (bar}—£1,650. 

interesting and progressive Careers to The appointment will be made as soon as possible and the 

Designers and Technicians in the —— candidate will be required to take up duty not later 
‘ than Ist January, 1960. 
Aircraft and Gas Turbine industries. 

If you are an A.F.R.Ae.S. or possess Applications, including the names of two persons to whom 

: , reference may be made, should be sent to the Registrar, The 

an equivalent qualification, University, Bristol, 8, from whom further particulars may be 


obtained. 


please write to: The Technical Staff Manager 


BLACKBURN & GENERAL AIRCRAFT LIMITED 
BROUGH, EAST YORKSHIRE 


C106/a 


FATIGUE 
DATA SHEETS 


AN INTRODUCTION TO 
AERONAUTICAL ENGINEERING 


Vol. 1: Mechanics of Flight 


A supplementary issue to the 


By A. C. Kermode, O.B.E., M.A. 6th Edition 

The well-known standard work on the subject Fatigue Data Sheets will be issued 

It forms an essential and extremely valuable | shortly. The new sheets will cover 

book for those engaged in all branches of , 

aeronautical work. 25/- net. | Basic Statistical Concepts and 
| Methods; the Effect of Inter- 


THE HELICOPTER AND 


HOW IT FLIES | ference in Pin Joints; Fatigue 
By John S. Fay. Meets the need for a simple 

book giving basic information about helicopter Testing Machines Data and 
theory without resorting to mathematical . 

formulae. Containing over 100 illustrations, this . Geometric Stress Concentration 


is a most useful book for the pilot, apprentice, 
and the technician transferring to helicopters. 
12/6 net. 


Factors 


FROM ALL BOOKSELLERS 
further particulars are available from 


PITMAN AERONAUTICAL BOOKS =| = THE ROYAL AERONAUTICAL SOCIETY 
Parker St., Kingsway, London, WC2 } 4 Hamilton Place, London W.1 
| 


Grosvenor 3515-9 


{ADVERTISEMENTS NOVEMBER 1959 


JOURNAL OF THE ROYAL AERONAUTICAL SOCIETY|\ 15 


: 
< 
Rep 
~ 
a 
Pe 
‘ 
; 
| 
9 
4° 
ee 
ay 
ae 
| 
| | 
in 
14 


TRADE MARKS SECTION 


AIRCRAFT MATERIALS LTD. 


STRUCTURAL MATERIALS 
and COMPONENTS 


FIRTH-VICKERS STAINLESS STEELS LTD. 


STAINLESS STEEL 


ASSOCIATED ELECTRICAL INDUSTRIES LTD 


ELECTRICAL EQUIPMENT 
for Aircraft 


Associated Electrical Industries Ltd. 


Aircraft Equipment Group 


HUNTING AIRCRAFT LTD. 


BOULTON PAUL AIRCRAFT LTD. 


KELVIN G HUGHES (AVIATION) LTD. 


(HENRY HUGHES & SON LTD. 
KELVIN, BOTTOMLEY & BAIRD LTD.) 


THE BRITISH REFRASIL CO. LTD. 


LIGHTWEICHT, HIGH TEMPERATURE INSULATION 


K.L.C. SPARKING PLUGS LTD 


KLG 


SPARKING PLUCS AND 
IGNITION EQUIPMENT 


COWTY GROUP LTD. 


DOWTY 


UNDERCARRIACES 
HYDRAULIC G ELECTRICAL EQUIPMENT 
FUEL SYSTEMS FOR CAS TURBINES 
RUBBER SEALS 
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LOCKHEED PRECISION PRODUCTS LTD 


hydraulics 


SAUNDERS-ROE LTD 


-SAUNDERS-ROE 


LIMITED 

OSBORNE EAST COWES - ISLE OF WIGHT 

Phone Cowes 2211 and at Trafalgar $448 
CFH/SKO/15 


JOSEPH LUCAS (CAS TURBINE EQUIPMENT) LTD 


SMITHS AIRCRAFT INSTRUMENTS LTD 


SMITHS 


SMITHS AIRCRAFT INSTRUMENTS LTD. 


THE UNITED STEEL COMPANIES LTD. 


REDIFON LTD. 
” 
Kediton RED FOX 
HEAT RESISTING 
FLIGHT SIMULATOR STEELS 
DIVISION S. FOX & CO. LTD. SHEFFIELD 

CRAWLEY - SUSSEX - ENCLAND mide 

ROTAX LTD WESTLAND AIRCRAFT LTD 


ROTAX 


AIRCRAFT ELECTRICAL ENGINEERS 


Ss- WESTLAND 
The Hallmark of British Helicopters 
Westland Aircraft Limited, Yooit, England 
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Aircraft Division Sales Department), Poyle 
é 
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Yeovil, Somerset. 
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The non-technical monthly for the enthusiast 


* 


Exciting exclusive articles on past and present 
aircraft 


* 


Superb photographs — scores of them in each 
issue — keep you up to the minute on new 
developments 


* 


Large scale drawings provide information on 
world’s aircraft 


At your bookseller or 


newsagent | /6d 
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Or send £1.2.6d for a year’s subscription to Dept. A.S., AIR 
PICTORIAL, ROLLS HOUSE, BREAMS BUILDINGS, LONDON, EC4 


if you would like a FREE COPY fill in the form 


AIR 
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Passenger Comfort Assured 
in the 

NORMALAIR vapour cycle cooling systems will provide full refrigeration even 
when an aircraft is stationary with engines stopped, or taxying with engines 
turning at low r.p.m. They can also serve as a means of reducing the relative 
humidity of the cabin air supply, thereby ensuring an acceptable effective temper- 
ature in the cabin. The separate ground conditioning trucks used with air cycle 
systems can, therefore, be dispensed with. 

As well as offering these important advantages, NORMALAIR vapour cycle 
systems make much smaller demands on power supplies than comparable air cycle 
systems. In general, they compare favourably with air cycle systems in the matter 
of weight, and where heat loads are high, are significantly lighter. 

In addition to these aircraft systems, which are available for aircraft of all classes, 
and which can be based on either reciprocating or centrifugal compressors, 
NORMALAIR are producing systems for guided missiles and other applications. 


enciano 


The only designers and manufacturers in 
Great Britain of complete cabin air conditioning 
systems for all classes of military and civil aircraft 


NORMALAIR (CAMADA) LTD TORONTO NORMALAIR (AUSTRALIA) PTY LTD MELBOURNE 


PRINTED BY THE LEWES PRESS. WIGHTMAN & CO. LTD.. LEWES. SUSSEX. ENGLAND. AND PUBLISHED 
BY THE KUYAL AERONAUTICAL SOCIETY. 4 HAMILTON PLACE. LONDON, W.1, ENGLAND. 


; 
| 
4 
4 
j \ 
4 
Ve 
\\ 
Ven, 
4 
\ 
art 
Ve: 
. 


: 
: 

{8 

; 4 

4 
|. 

| 
5 

\ BE: 
we 

=i 

; 
i 

= 
: 

os 

ae 

| 

: 


